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SUMMARV 

The  development  of  a mathematical  model  of  the  nitration  section  of 
a continuous  process  for  the  manufacture  of  TNT  Is  described.  Imple- 
mentation of  the  model  on  a high-speed  digital  computer  is  also  discussed , 
as  are  the  steady  state  and  dynamic  studies  that  were  conducted  on  the 
resulting  computer  simulation.  The  way  in  which  optimized  process  oper- 
ating conditions  and  improved  process  control  strategies  were  developed 
using  the  simulation  is  Illustrated. 
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INTRODUCTION 


Computer  simulation  of  a chemical  process , simply  defined , is  a 
procedure  by  which  a mathematical  model  of  the  process'  relevant  chemi- 
cal and  physical  characteristics  is  first  constructed,  after  which  the  re- 
sultant system  of  equations  is  programmed  for  solution  on  a high  speed 
computer.  Such  a model  is  usually  based  upon  a combination  of  prevailing 
theory  and  empirical  relationships,  and.  when  properly  formulated,  will 
predict  all  process  input/output  relationships . The  applicability  of  such 
a model  dapenda  heavily  upon  the  ratio  of  theory  to  empiricism  employed 
in  model  development.  For  example,  a totally  empirical  model  such  as  a 
least  squares  polynomial  approximation  Is  useful  only  within  the  bound- 
aries of  experimental  observation,  whereas  a totally  theoretical  model, 
for  example  Newton's  Laws  of  Motion  as  a description  of  projectile  trajec- 
tories. will  apply  under  virtually  all  conditions.  For  moat  complex  systems 
(such  as  a chemical  process)  an  optimum  balance  between  theory  and 
observation  must  be  sought. 

Upon  completion  of  mathematical  modeling  and  computer  program- 
ming, the  resulting  computer  simulation  can  then,  ofter  verification  against 
real  process  data,  be  used  to  study  the  process  without  incurring  the  cost, 
time  loss,  and  (especially  in  the  case  of  explosives  manufacturing  processes) 
the  risk  involved  in  pilot  plant  or  full-scale  plant  experimentation.  Most 
often,  studies  of  both  the  static  and  dynamic  states  of  the  process  are  carried 
out,  the  objective  of  the  former  being  optimisation  of  operating  and  design 
conditions,  and  the  objective  of  the  latter  being  development  and  testing  of 
control  strategies . 

It  has  always  been  possible  to  develop  mathematical  models  of  com- 
plex physical  systems.  However,  the  value  of  such  models  remained  ques- 
tionable until  the  application  of  digital,  analog,  and  hybrid  computers  made 
possible  the  rapid  solution  of  the  complex  equations  which  usually  resulted. 
It  has  only  been  in  the  lust  10  to  IS  years  that  mathematical  modeling  of  com- 
plex prooeesea  haa  yielded  elgnlficant  accomplishments.  Moat  notable  among 
these  have  been  the  simulations  used  In  the  aerospace  industry,  especially 
in  the  manned  space  flight  program. 

Although  employed  luccesefully  as  a process  study  tool  for  a number 
of  years  in  private  industry , the  application  of  chemical  process  simulation 
to  the  U,  S,  Army  Armament  Command's  explosive  and  propellant  manufac- 
turing proceseee  is  Just  now  beginning  to  be  appreciated.  For  the  past 
three  and  a hall  years,  the  Chemical  Process  Technology  Division  of  the 
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Manufacturing  Technology  Directorate  at  Plcatlnny  Arsenal  has  been  en- 
gaged In  an  extensive  effort  to  Improve,  via  the  application  of  computer 
simulation,  a continuous  process  for  TNT  manufacturing.  This  work  has 
been  carried  out  Jointly  with  Imperial  Chemical  Industries,  U.  S.  (ICIUS) 
the  operating  contractor  at  Volunteer  Army  Munition  Plant  in  Chattanooga, 
Tennessee,  and  has  to  date  gone  through  three  distinct  phases;  Phase  I, 
initial  development  of  a kinetics  and  mass  transfer  model  of  the  nitration  sec- 
tion with  subsequent  implementation  on  a hybrid  computer;  Phase  II,  con- 
version of  the  hybrid  nitration  model  to  a pure  digital  model , followed  by 
improvement  of  the  nitration  model  based  on  new  plant  data  and  additional 
theoretical  oonalderations , simulation  of  nitration  section  hydraulics,  and 
development  of  a simulation  of  the  purification  section  of  the  process;  and 
Phase  III,  extensive  exploitation  of  the  nitration  and  purification  computer 
models . 

Project  summary  reports  (lief  1 , 2 , 3)  have  previously  been  prepared 
covering  these  efforts,  However,  no  unified  text  has  been  generated  which 
ties  together  all  the  Interrelated  tasks  and  provides  a much  needed  explan- 
ation of  the  present  structure,  capabilities,  and  limitations  of  various  models 
developed.  The  purpose  of  this  report  la  to  satisfy  these  needs  as  well  as 
to  clarify  material  in  the  existing  reports  which  has  either  become  out- 
dated or  which  the  author  feels  is  unclear  as  presently  stated . 

Because  of  the  extensive  amount  of  information  which  is  prosented , 
the  report  has  been  prepared  in  three  volumes.  Volume  I describes  the 
modeling  and  almulation  effort  that  hao  been  carried  out  on  the  nitration 
section  of  the  process,  Includiiig  a description  of  the  process  chemistry 
and  equipment.  Volume  II  deals  exclusively  with  the  purification  section 
model,  and  Volume  III  involves  the  almulation  of  the  hydraulic  phenomena 
occurring  in  a nitration  stage. 


THE  CHEMISTRY  OF  TNT  NITRATION  - A BRIEF  DESCRIPTION 

The  production  of  TNT  is  essentially  carried  out  via  a thrne-etep 
reaction  in  which  toluene  and  its  mono-  and  dlnltro-substltuted  derivatives 
are  contacted  with  mixtures  of  nitric  and  sulfuric  uold  called  mixed  acid . 
Bach  step  in  the  reaction  corresponds  to  the  successive  ^lectrophillc  sub- 
stitution of  a hydrogen  ion  (H  ) by  a nltronlum  ion  (NO|  ) on  the  aromatl^ 
ring  of  the  toluene  molecule.  The  nitric  acid  serves  as  the  source  of  NO| 
ions,  while  the  sulfuric  acid  acts  as  both  a catalyst  and  dehydrating  agent, 
first  reacting  with  the  nitric  acid  to  generate  NO|  ions  and  then  picking  up 
the  water  which  forms  as  a product  of  the  nitration  reaction. 


A simplified  scheme  for  the  mononitration  of  toluene  in  mixed  acid 
can  be  written  as  follows  : 


HNO,  + H,S04  — NO,  + HSO,  + H,0 


^ HSO4 


(Toluene) 


(o-MNT) 


+ H,S04  (2) 


Adding  reactions  1 and  2 gives  the  overall  reaction: 


* * NO 

+HNO,  +H,0 


Similarly,  when  toluene  is  replaced  by  mononitrotoluene  (MNT) 
or  dinitrotoluene  (DNT) , the  stoichiometric  reactions  are  written  as 


+ HNOj 


H,S04 


CH, 

^NO, 


(2,4-DNT) 


CH, 

h so  o,n^^  no, 

*HNO.  .H, 

Kim 


no, 

(a -TNT) 
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The  temperature  as  well  as  the  composition  of  mixed  acid  required 
to  efficiently  carry  out  each  of  the  above  reactions  is  different,  reflecting 
the  relative  ease  with  which  nitration  Is  accomplished.  Thus,  Heaction  3 
is  usually  conducted  at  relatively  low  temperatures  (40°  to  50°C)  with  weak 
acid  solutions,  while  reactions  4 and  5 require  progresBlvely  highui  temper 
Btures  (60°  to  100°C)  and  stronger  acid  mixtures  to  be  driven  to  the  desired 
degree  of  completion. 

The  strength  of  mixed  acid  is  usually  determined  by  the  amount  of 
water  present,  In  aqueous  media  Reaction  1 oont  ols  the  generation  of 
nltronium  ions.  However,  when  nitration  is  conducted  in  an  anhydrous 
medium,  the  acid  mixture  consists  of  nitric  acid  and  oleum  (100%  lltSO^ 
in  which  SO|  dissolved)  and  the  controlling  nltronium  ion  generation 
equilibria  are: 

HNO,  + 2H,S,0|  ^ NOJ  + HS,0^  + 2H,S04  (B) 


HNOa  + HS,Ot 


NOji  4 2HSO| 


It  is  also  worth  noting  that  Reactions  3,  4,  and  5 are  highly  exothermic. 

Although  not  shown  In  Reactions  3 through  6,  a significant  amount 
of  isomerism  occurs , resulting  in  the  production  of  a variety  of  nitrated 
components  in  addition  to  the  desired  o-TNT . A complete  description  of 
the  Isomerism  in  toluene  nitration  with  relative  amounts  of  each  component 
formed  is  given  in  Figure  1.  The  relative  amounts  of  isomers  formed  de- 
pend upon  nitrating  temperature  and  add  composition.  Thus  Figure  1 must 
be  viewed  as  a typical  distribution  only.  In  the  continuous  TNT  process  at 
Radford  Army  Ammunition  Plant,  only  the  a,  p,  y,  and  e isomers  have  been 
positively  Identified,  although  by  the  analytical  technique  used  on  process 
samples,  it  Is  not  possible  to  separate  the  a and  q isomers  (Ref  4)  . 

In  addition  to  its  function  ss  a nitrating  agent,  nltrlo  acid  also  plays 
the  role  of  an  oxidant  under  the  usunl  conditions  found  In  toluene  nitra- 
tion reactions.  All  nitrated  components  sre  capable  of  undergoing  some 
form  of  oxidation.  Two  types  of  oxidation  reactions  are  observed;  those 
which  result  in  destruction  of  the  aromatic  ring,  yielding  gaseous  pro- 
ducts, and  those  which  result  in  oxidation  of  the  methyl  group  or  hydroxyla- 
tlon  of  the  aromatic  ring,  yielding  stable  organic  compounds,  Typically, 
MNT  undergoes  hydroxylation  to  the  cresol  which  is  later  destructively 
oxidized  to  gaseous  products;  DNT  undergoes  complete  ring  destruction 
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to  yield  COj , CO,  NO,  NOj  and  tetranltromethane  (TNM);  and  TNT  under-  O 

goes  methyl  group  oxidation  to  yield  a variety  of  organic  products  as  shown 
in  Figure  2. 

In  a typical  nitration  of  toluene  to  TNT , isomerism  and  oxidation 
account  for  a yield  loss  of  about  8 to  10%  based  on  molar  feed  of  toluene . 


PROCESS  DESCRIPTION 

The  continuous  TNT  process,  which  is  the  subject  of  the  simulation 
study,  is  an  adaptation  of  technology  developed  by  Canadian  Industries 
Limited  and  is  currently  the  only  continuous  process  for  manufacture  of 
TNT  operating  in  the  United  States. 

Nitration  is  carried  out  in  a series  of  six  reaction  stages , each  con- 
sisting of  a 425-gallon  nitrator  and  a ISO-gallon  separator,  except  for  Stages 
1 and  3 which  have  two  nltrators , Because  of  the  ratios  of  mixed  acid  to 
toluene  employed  in  the  process,  a two-phase  system  results,  with  nitra- 
tion occurring  in  the  acid  phase  as  diffusion  and  solubility  phenomena 
establish  equilibrium  concentrations  of  the  organic  species  in  that  phase , 
Separation  of  organic  and  acid  phases  is  carried  out  after  each  nitration 
stage  and  enables  countercurrent  flow  of  the  reactants  to  be  maintained. 

A flow  diagram  of  the  nitration  section  of  the  process  is  shown  in  Figure  3. 
and  a detailed  diagram  of  a nitrator/separator  nitration  stage  is  shown  in 
Figure  4. 

In  the  countercurrent  process,  all  of  the  toluene  is  introduced  into 
the  first  stage,  where  mononitration  and  a small  amount  of  dinltratlon  oc- 
curs. After  separation,  this  so-called  "nitrobody"  stream  moves  forward 
through  the  reactor  train , becoming  more  highly  nitrated  as  it  progresses 
from  vessel  to  vessel.  Oleum  is  fed  to  Stage  6 and  by  moans  of  the  mixed 
acid  recycle  streams,  moves  in  the  opposite  direction,  decreasing  in  strength 
as  it  flows  toward  Stage  1 . Nitric  acid  is  added  to  each  stage  in  order  to 
maintain  the  correct  mixed  acid  composition.  Mixed  acid  recycle  streams 
also  contain  some  dissolved  nitrobody  as  well  as  nltrosylsulfurlc  acid 
(HNOSO4 ) which  forms  as  a result  of  the  oxidation  side  reactions . 

As  pointed  out  previously,  nitration  becomes  progressively  more 
difficult  as  the  number  of  nitre  groups  on  the  aromatic  ring  of  the  sub- 
strate increases.  Thus,  the  temperature  as  well  as  the  mixed  acid  strength 
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must  Increase  from  Stage  1 to  Stage  6 . Nominal  temperatures  for  the  eight 
nitratorsare  55,  50,  70,  80,  85,  90.  85,  and  100°C  respectively. 

From  Figure  4 it  is  seen  that  a nltrator  is  a cylindrically  shaped 
vessel  with  an  agitator  fitted  inside  a draft  tube  at  the  center . The  agi- 
tator provides  the  required  interphase  mixing  in  the  nitrator  and  also 
creates  a head  difference  between  material  on  either  side  of  the  draft  tube 
which  in  turn  provides  the  motive  force  for  fluid  flow . There  are  no  pumps 
in  the  system.  Rapid  Internal  circulation  of  material  over  the  nitrator 's 
cooling  coils  is  maintained . There  is  also  a recycle  of  acid  between  the 
separator  and  its  own  nitrator  (internal  recycle) . Plow  of  this  stream  can 
be  adjusted  to  maintain  proper  cooling  in  the  separator  or  to  change  the 
acid*‘tO“nitrobody  ratio  and  residence  time  in  the  nitrator . 

The  acid  stream  flowing  from  a separator  to  the  previous  nitration 
stage  (external  recycle]  passes  through  a weir-‘llke  device  called  a decanter. 
By  adjusting  this  apparatus,  the  nitrobody/acid  interface  in  the  separator 
can  be  regulated  to  insure  adequate  phase  separation . However , under 
normal  operating  conditions  of  55  tons  per  day  TNT  production , a signifi- 
cant amount  of  entrainment  still  occurs,  especially  in  Separator  3, 

The  crude  TNT  leaving  Separator  C flows  to  the  purification  section 
where  residual  nitrating  acid  and  asymmetrical  isomers  are  removed  prior 
to  final  packaging  of  the  product. 


DEVELOPMENT  OF  THE  MATHEMATICAL  MODEL 


The  Single  Vessel  Model 


The  majority  of  the  work  conducted  under  the  TNT  process  computer 
simulation  effort  has  been  devoted  to  the  construction  of  a mathematical 
model  which  describes  the  nitration  section  of  the  process  in  terms  of  the 
kinetic  and  mass  transfer  phenomena  which  are  believed  to  occur  there 
and  account  for  most  of  the  process'  behavior.  As  previously  Indicated, 
models  of  the  purification  section  and  of  the  hydraulic  phenomena  which 
govern  fluid  flow  in  a given  nitration  stage  have  also  been  developed  but 
will  not  be  discussed  in  this  volume. 


The  kinetic  and  mass  transfer  model  of  the  nitration  section  of  the 
process  (hereafter  called  the  nitration  section  model)  is  the  end  product 
of  numerous  refinements  and  Improvements  which  have  been  made  to  a 
first- generation  model  developed  for  Picatinny  Arsenal  in  May  1970  by 
Kinotrol,  Inc.  of  Houston.  Texas,  functioning  as  a consultant  to  ICIUS 
at  Volunteer  Army  Ammunition  Plant  in  Chattanooga,  Tennessee.  The 
structure  and  contents  of  the  original  nitration  section  model  can  be 
found  in  Reference  1 . 

As  in  the  orieinal,  a fundamental  approach  was  taken  in  developing 
the  present  nitration  section  model.  Wherever  possible,  theoretical  equa- 
tions are  employed  to  represent  the  chemical  and  physical  phenomena 
which  govern  the  nitration  process.  The  use  of  empirically  derived  equa- 
tions and  correlations  is  kept  to  a minimum  and  has  been  resorted  to  only 
after  exhaustive  research  for  theoretical  alternatives  proved  fruitless. 
Phenomena  described  by  the  model  include  kinetics  of  nitration  and  oxi- 
dation of  organic  species,  diffusion  of  the  organic  species  between  the 
acid  and  organic  phases,  solubility  of  organic  species  in  tho  acid  phase, 
and  solubility  of  nitric  acid  in  the  organic  phase,  The  way  in  which  real 
plant  data  was  used  in  the  development  of  the  nitration  section  model  will 
be  described  when  parameter  fitting  is  explained  in  a later  section  of  this 
report. 

The  easiest  way  to  explain  the  nitration  section  model  is  to  take  a 
single  vessel  in  the  nitration  train  (there  are  14  vessels  in  series)  and 
develop  a complete  material  balance  for  each  component  in  both  acid  and 
organic  phases . The  generalized  nitrator  in  Figure  5 showing  actual 
stream  flows  could  be  used  for  this  purpose.  However,  the  analysis  is 
far  simpler  if  each  stream  is  broken  down  into  its  acid  and  organic  phases 
with  the  subsequent  addition  of  similar  phases , so  that  the  vessel  input 
and  output  consist  of  a single  acid  phase.  A,  a single  organic  phase  P, 
and  a gas  phase,  G.  Figure  6 illustrates  this  simplifying  concept.  Thus, 
for  component  i in  the  acid  phase,  a mole  balance  can  be  written  as; 

(Rate  of  moles  i in  with  A ] + 

0 

(Rate  of  generation  of  moles  1 (Rate  of  accumulation  of 

in  the  acid  phase  within  moles  1 within  the  acid 

nitrator  n)  - of  nitrator  n) 

(Rate  of  moles  i out  with  A ) (1) 
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The  balance  for  component  i In  the  organic  phase  follows  the  same  form  as 
Equation  1.  For  the  gas  phase,  the  balance  is  simplified  by  the  fact  that 
products  of  the  gaseous  reaction  are  insoluble  in  both  the  acid  and  organic 
phases . Thus , 

(Rate  of  generation  of  (Rate  of  gaseous  i out 

gaseous  species  1 in  b with  0 ) 

nitrator  n)  (2) 


For  a specific  component,  say  a-TNT,  the  molar  balances  for  vessel  n 
are  written  as: 

“A  « ,(x~TNT 

n, a-TNT  n,o-TNT  n,o-TNT  dt  ’ 


1* 

-p 


for  the  acid  phase,  where  S is  the  acid  phase  generation  term,  and 

dP 


-f  - P B iP-TNT 

n, a-TNT  n,o-TNT  n, a-TNT  dt  ' ' 


for  the  organic  phase , where  S is  an  organic  phase  generation  term , 

The  assumption  of  a perfectly  backmixed  reactor  allows  for  the  equivalency 
of  the  reactor  contents  and  the  reactor  outflow . 


Equations  similar  to  3 and  4 can  be  written  for  all  the  components 
in  each  phase  for  every  vessel  In  the  nitration  train . In  the  current  ver- 
sion of  the  nitration  section  model , 10  chemical  components  are  considered 
in  the  organic  phase,  14  in  the  acid  phase,  and  3 In  the  gaseous  effluent 
(see  Figure  8) . Before  proceeding  further,  a few  comments  on  these  com- 
ponents are  in  order.  First,  all  nitrobody  components  originating  from 
m-MNT  have  been  lumped  together  in  them-DNT  and  m-TNT  terms  in  order 
to  simplify  the  kinetics  of  the  asymmetrical  isomer  reactions.  Secondly, 
a single  term,  TNBX,  is  used  to  represent  all  components  other  than  TNB 
formed  from  methyl  group  oxidation  of  TNT.  This  Includes  trinitrobenzyl 
alcohol,  trlnitrobenzaldehyde,  and  trinltrobenzolc  acid.  Lumping  of  oxi- 
dation products  into  one  term  was  necessitated  by  a lack  of  process  data 
for  these  components , and  is  justified  considering  the  small  concentration 
of  each  component  which  is  present. 

Thus,  from  Equations  3 and  4,  It  can  be  seen  that  for  a given  vessel 
the  model  consists  essentially  of  a series  of  first  order , non-linear , coupled , 
ordinary  differential  equations  which  make  up  the  vessel's  material  balance. 


:t 
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The  way  in  which  the  individual  vessels  are  connected  and  the  procedure 
for  generating  the  flow  and  composition  of  the  actual  process  streamu 
from  the  ideal  phases  will  be  discussed  at  length  in  later  sections  of  this 
report , 


The  heart  of  the  vessel  material  balance  equations  and  therefore 
the  heart  of  the  nitration  section  model  is  the  generation  terms . These 
terms  are  actually  equations  which  represent  the  kinetic , mass  transfer , 
and  solubility  processes  occurring  within  a given  vessel . For  each 
component  balance  equation,  three  generation  terms  exist;  reaction  rate 
R .) , diffusion  (D  ) , and  bulk  mass  transport  (M  .) . 


Thus  Equation  3 can  be  rewritten  as: 

A®  + 

n.oTNT  “n.oTNT  n.oTNT  "n.oTNT 


\,aTNT 


dA 

, n.oTNT 
dt 

(5) 


Since  it  is  widely  accepted  that  all  reactions  occur  in  the  acid  phase,  it 
follows  that  ^ B 0,  Also  mass  transfer  into  the  acid  phase  must  be  at 
the  expense  or  Aiaterlal  in  the  organic  phase  so  that  oTNT  * “I^^.oTNT 
and  m|  ' “M^.oTNT-  Thus,  Equation  4 can  be  written  as: 

v°  - n - M - P ^ j^n  .cTNT 

*^n,oTNT  n.aTNT  n.oTNT  *^n,aTNT  dt 

(6) 

Note  that  there  is  no  longer  a need  for  acid  and  organic  phase  superscripts 

on  the  mass  transfer  terms . Development  of  the  generation  terms  R^  | 

0 , and  M , will  now  be  discussed  in  detail.  ' ’ 

n 1 1 , n 1 1 


Reaction  Kinetics 

Reaction  rate  expressions  have  been  developed  to  represent 
the  nitration  and  oxidation  of  toluene  and  its  nitro- substituted  derivatives 
at  each  stage  in  the  reactor  train . The  stoichiometric  reations  which  are 
considered  in  the  nitration  section  model  Include: 
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Nitration 


Toluene  + HNO*  — ► uMNT  ♦ H,0 
Toluene  + HNO|  — ► mMNT  + H,0 

aMNT  + HNO,  oDNT  + H,0 
mMNT  + HNO,  — > mDNT  + H,0 

ODNT  + HNO,  — ► oTNT  + H,0 
mDNT  + HNO,  — ► mTNT  + H,0 

Oxidation 

DNT  + 13HNO,  + 11H,80«  - TNM  + 6COu 

* IIHNOSO4  1SH,0 

DNT  + 12iHNO,  + 144H,S04  — ► 7COu 

* 14IHN0S04  ^ 16iH,0 

HNOSO4  H,0  + r [01  HjSO, 


(Reaction  1) 
(Reaction  2) 

(Reaction  3) 
(Reaction  4) 

(Reaction  S] 
(Reaction  6) 


(Reaction  7} 


(Reaction  8) 


+ NO  -0  HNO,  (Reaction  0) 

4TNT  + 9HNO,  + SHjSO,  TNBOH  + TNBAL 

+ TNBA  + TNB  + CO,  + 12H,0  + BHNOSO,  (Reaction  10) 

It  Is  pointed  out  that  Reactions  1 through  10  are  not  mechanistic  equations 
but  rather  stoichiometric  equations  which  account  for  the  appearance  and 
disappearance  of  the  various  chemical  components  during  the  nitration 
process . 

As  mentioned  briefly  In  the  section  on  process  chemistry,  aromatic 
nitration  in  mixed  acid  occurs  via  electrophlllic  substitution,  where  In 
most  instances  the  electrophile  has  been  proven  to  be  the  nitrruium  Ion. 

However,  It  has  also  been  well  establlahed  (Ref  6)  that  as  the  molar 
ration  of  water  to  sulfuric  acid  in  the  mixed  acid  Increases  to  unity , the 
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concentration  of  NO|  falls  off  rapidly  to  nearly  apectroacoplcally  undetect- 
able  levels.  Thus,  in  the  case  where  [H|0]/[H|S04]  weak  acid 

region)  and  nitration  la  still  seen  to  occur  at  a rapid  rate , either  the  aro- 
matic substrate  is  reactive  enough  to  require  only  trace  amounts  of  nitron- 
iuro  ion  for  reaction  to  proceed , or  another  nitration  mechanism , say  one 
dependent  on  molecular  nitric  acid  is  in  ehfect.  The  former  is  more  likely 
but  the  latter  cannot  be  ruled  out . 

In  the  continuous  TNT  process,  all  of  the  monitratlon  and  a signifi- 
cant amount  of  dlnltratlon  occurs  in  mixed  acid  where  [M|0]/[H|S04]  > 1. 
so  that  in  developing  a kinetics  model  for  this  weak  acid  region  one  of  the 
two  mechanism  hypothesis  must  be  chosen  . It  was  first  decided  to  test 
both  hypotheses  against  actual  process  performance  data . While  the  nltron- 
ium  ion  mechanism  was  eventually  shown  to  be  preferable  and  was  subse- 
quently included  in  the  overall  nitration  section  model,  a detailed  descrip- 
tion at  the  weak  aoid  model  baaed  on  molecular  nitric  acid  ia  nevertheless 
also  presented  here  for  reasons  of  completeness . 

The  molecular  nitric  dependent  weak  acid  kinetics  model  was  based 
on  the  work  of  Kobe  and  Lakemeyer  (Ref  7)  and  McKinley  and  White  (Ref  8) 
who  investigated  the  factors  affecting  the  mononltratlon  of  toluene  in  mixed 
aoid  where  the  molar  ratio  of  water  to  sulfuric  acid  was  greater  than  two, 

A mechanism  is  proposed  whereby  the  nitrosonlum  ion , NO  , is  first  formed , 
then  attacks  the  benzene  ring  and  is  subsequently  oxidized  by  nitric  acid 
to  yield  the  mononitrated  product.  This  mechanism  can  be  represented  as 
follows: 


Although  not  definitively  stated  in  the  literature , Equation  B appears  to 
be  the  rate  limiting  step  since  data  taken  by  McKinley  and  White  show  a 
strong  dependence  on  the  concentration  of  molecular  nitric  acid  in  the 
acid  phase.  McKinley  and  White's  nitration  studies  also  indicate  that  the 
sulfuric  acid  concentration  has  a significant  influence  on  the  rate  of  toluene 
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nitration  In  weak  acid.  Finally,  Kobe  and  takemeyer  demonstrated 
that  the  presence  of  nitroaylaulfurlc  acid  in  weak  nitrating  acid  mixtures 
causes  a maximum  in  the  rate  at  4 mole  % HNOSO4 . 


All  of  the  above  observed  behavior  was  combined  to  yield  an  expression 
for  the  bimolecular  rate  constant  for  nitration  in  weak  acid.  This  expres- 
sion should  be  applicable  to  toluene,  MNT,  or  DNT  as  long  as  the  condi- 
tions of  a weak  acid  nitrating  medium  are  met.  The  form  of  this  expresaion 
is: 


RatOj^- 


k,  • IH,S04]®  i0"“*^HNOSO4  ' (10) 

where  [M]  la  the  concentration  of  the  organic  reactant  in  the  acid  phase; 

[N]  is  the  concentration  of  the  organic  product  in  the  acid  phaae; 
^HNOSOi  fraction  nitrosylsulfuric  acid  in  the  acid  phase; 

and  a and  0 are  adjustable  parameters. 


The  dependence  of  weak  acid  nitration  reactions  on  temperature  was  ex- 
pressed through  an  Arrhenius-type  relationship  and,  reaction  was  limited 
to  the  acid  phase.  Thus,  Equations  8 and  10  can  be  restated  as; 


- fs 

.«M  .K  , H '•T  Tb' 

_1  A.  •'k*  ^ 


[M]  8 


8*V,  [HNO,]  [HjSO*]®  lO'**  ^^HNOSO,  ‘ 


where 


is  the  volume  of  the  acid  phase; 
is  the  Arrhenius  frequency  factor; 
is  the  activation  energy;  and 
Tj^  is  a reference  temperature . 

Plant  dita  fc  r the  continuous  process  at  Radford  indicates  that  toluene 
will  react  to  .itrotoluene  juat  aa  fast  as  it  can  diffuse  into  the  acid 
phase.  Thus  for  Reactions  1 and  2,  mass  transfer  is  the  limiting  factor  and 
the  weak  acid  kinetics  developed  above  are  of  academic  Interest  only . How- 
ever, for  conversion  of  MNT  to  DNT  in  weak  acid,  the  weak  sold  mechanism 
would  apply . 

The  above  model  is  obviously  based  on  rather  unstable  theoretically 
grounds  and  even  though  a reasonably  good  fit  to  plant  data  was  achieved 
with  it,  it  was  scraped  in  favor  of  a more  theoretically  consistent  nitron- 
ium  ion  dependent  system  of  equations. 

In  this  approach  the  reaction  rate  la  still  expressed  as  in  Equation  9; 
however  the  expression  for  k|  is  simply: 


k| 


''.■'nQ 


(12b) 


where: 


E 


K ■ k 

N 


.::n  (A. . JL) 

e R NT  tJ 


(12b) 


tNoj'  1 
Q “ (HNO,  1 


(12c) 


Utilisation  of  this  model  depends  on  the  ability  to  compute  the  nitronlum 
concentration  from  the  equilibrium  reactions  which  are  in  effect.  The 
controlling  equilibria  are  chosen  to  be  (Ref  18) ; 


HNO,  + H,S04 
H,0  + H,S04  • 


NO,  -t-HSO,  +H,0 


H|0  +HSO4 


(12d) 

(12e) 


14 


o 


The  prooedure  by  which  the  NO|  concentration  Is  calculated  follows  below. 

In  the  weak  acid  region . the  Initial  and  equilibrium  concentrations 
are  given  as: 


Comoonent 

Initial 

At  Equilibrium 

HNO| 

[HNO,]^ 

[HNO,] 

HtSO« 

[H,S04l^ 

[H18O4] 

H|0 

NOJ 

0 

[H,q 

[NOjl 

H,0* 

0 

(H,0*) 

H8O4 

0 

mol] 

By  writing  material  balances  foi;  each  chemical  element,  a ayatem  of  equa- 
tlona  can  be  developed  which,  In  conjunction  with  the  equilibrium  Baua- 
tlona  12d  and  12a,  can  be  uaed  to  solve  for  the  desired  NO|  oonoentratlon 
equilibrium . The  elemental  material  balances  are  baaed  on  the  necessary 
equivalence  of  the  moles  of  chemical  elements  prasem  Initially  and  the 
moles  of  chemical  elements  present  at  equilibrium . 

Per  Nitrogen; 

(HNO,)^ 

» (HNO,J  + iNoJ) 

(12f) 

For  Sulfur: 

[H»804l^ 

° [H,S04)  ^ [HSO4] 

(12g) 

For  Oxygen; 

[H»0]q  + 

3 IHNO,]^  + 4[H,804lp  ■ 

[H,0] 

+ 

3[T  4(H,SO«]  ^ 2[N0^1 

+ 

4[HSOj  + 

(12h) 
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Making  uaa  of  Bquatlono  12f  and  12g  in  Equation  12h  glvea; 


+ {Nojl  ■ [H,Ol  + [H,0*1 


For  Hydrogan: 

2(H,Oljj  + [HNO,]^  + 2lH,804ljj  • 2[H,Ol 

* (HNO|l  -f  2[H|SO«]  4 (H8O4] 

♦ 3tH,0^1  (12J) 

Multiplying  Equation  121  by  2 and  aubtraotlng  from  Equation  12J  glvaa: 

[HNO,l^  - 2[NoJ]  ♦ 2[H,804)j^  ■ tHNO»l 


4 2[H|804)  ^ [H8O4)  [H|0^] 


(12k) 


And  aubtraotlng  Equatlona  12f  and  g from  Equation  12k  and  glvai; 
(H8O4I  • [H,0^1  (NO^ 


Now,  lat  tha  initial  oonoantrationa  ba  givan  aa; 


IH.OI. 


lH,804lg  ■ 8 


|HNO,]p  ■ F 


and  tha  aquillbrlum  oonoantrationa  of  nltronlum  and  blaulfata  Iona  be 
given  aax  and  y raapaotlvaly . Than  at  aquillbriuro,  tha  oonoantration 
of  tha  remaining  oomponanta  oan  be  axpraaaad  aa; 


[HNOil  >P>x  (fromBql2f) 

[H1SO4]  B s-y  (from  Eq  12g) 
lH|0^1  *•  y-x  (from  Bq  122) 
lH|0]  ■ W X - (y-x)  (from  Bq  121) 

Tha  aquillbrlum  oonatant  for  Equation  12d  la  given  byi 

(NOil  (H80;i  [H|Ol  (X)  (y)  (W42x-y) 

■ IH18O4I  (HNO,l  ■ (a-y)(P-x) 

and  that  for  Equation  12a  by: 

[Bio'll  (HSU)  (y-x)  (y) 

K|  ■ [h,01(H|S04]  ■ (W+2x-y)  (S-y) 


(14) 


Baaad  on  tha  ImpoaelbUlty  of  nagative  oonoantrationa . tha  aaaumed  axla- 
tanca  of  all  poaalbla  oomponenta  at  equilibrium , tha  following  conatralnta 


rouat  apply: 

0 < X < y 

If  y <F 

0 < X < F 

If  F<  y 

0 < y < (W*2x) 

if  (W^2x)  < S 

0 < y < S 

If  S < (W+2x) 

Knowing  K| , K| , F,  W,  and  S, 

and  utllialng  the  above  constraints , Bqua- 

tlona  13  and  14  can  be  aolved  for  x,  tha  nltronlum  Ion  concentration  at 
aqullibrluin , From  Equation  14, 

y»  - K,  (8-y)  (W-y) 

X-  2K,  (8-y)+y  ' ' 


Equation  IS  can  be  aolved  for  x and  y by  an  Interval  halving  procedure 
after  which  the  solution  values  can  be  eubatltutad  In  Equation  13  tu  give 


a value  of  . If  convergence  on  is  not  achieved , the  procedure  la 
repeated  using  the  current  value  of  y to  adjust  thr  upper  or  lower  boundary 
of  the  solution  interval  for  x.  In  actual  execution  of  the  model,  between 
20  and  40  iterations  are  required  to  converge  on  Ki  for  each  pass  th)  ough 
the  vessel  equations , 

Before  moving  on  to  the  strong  acid  region  it  should  be  pointed  out 
that  a general  mechanism  for  aromatic  nitration  in  aqueous  mixed  acid  has 
been  proposed  (Ref  8)  based  on  activity  coefficients,  but  lack  of  data  pre- 
vented using  it  in  the  nitration  section  model.  An  attempt  was  made,  how- 
ever, to  do  without  any  theoreMcally  based  weak  acid  mechanism  and  Just 
use  a curve  fit  of  laboratory  data  developed  by  Vinnik  (Ref  10)  for  nitra- 
tion of  p-MNT  in  mixed  acid.  The  lower  and  of  the  curve  had  to  be  extra- 
polated since  no  data  was  provided  for  H1SO4  conoentratlons  lass  than  80%. 
In  all  cases,  the  weak  acid  mechanism  of  Equations  8,  12a,  12b  and  12c 
gave  far  better  agreement  with  observed  plant  behavior,  and  It  la  these 
equations  that  constitute  the  currant  model . 

For  strong  acid  nitration  (l.e.  [HaOl/fHiSOil  < 1)  It  is  firmly 
established  that  the  nitronlum  ion  is  the  active  nitrating  species.  Bennett 
(Ref  11)  developed  an  expression  for  the  blmolecular  nitration  rate  con- 
stant, k| , In  the  reaction 


- k,  IHN0,](DNT]  (16) 

whsre  the  dependence  on  the  NO|  Ion  concentration  as  well  as  the  concen- 
tration of  proton  acceptor  oomponenta  (HSO« , HfSOi , and  HS1O7  } are 
considered.  Bennett's  expression  is  given  as  follows: 

k|  ■ Q (k'  [HSOj  +k"  (H,S04]k"  (HS,0,‘l^  (17) 

(NO^ 

where  Q - [hno,] 


In  order  to  use  Bennett's  equation  to  represent  strong  acid  nitration 
it  la  necessary  to  determine  the  nitronlum  and  proton  acceptor  Ion  concen- 
trations throughout  the  nitration  section  of  the  process . In  the  model, 
four  strong  acid  regions  sre  defined  in  which  these  ions  should  exist,  and 
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U than  equilibrium  equations  are  establlahed  which  allow  calculation  of 

the  ionic  component  concentrations.  The  first  two  regions  are  exactly 
those  dealt  with  by  Bennett.  However  the  third  and  fourth  regions  have 
apparently  never  been  examined  experimentally  so  that  the  hypothesis 
presented  below  constitutes  an  extension  of  existing  theory  . 

The  first  region  (aqueous  region)  Is  ohsraoterlzed  by  the  presence 
of  water  subject  to  the  constraint  [H|0]/[H|804  ] < 1 , The  equilibria 
established  In  this  situation  are  exactly  those  of  Equations  12d  and  12e 
for  which  the  procedure  to  determine  the  equilibrium  concentrations 
has  been  given. 

The  remaining  three  regions  are  water  free , with  nitric , sulfuric , 
and  pyroaulfurlo  (H|S|Ot)  acids  making  up  the  nitrating  acid.  These 
will  be  referred  to  as  oleum  regions. 

In  the  first  oleum  region  (nitric  limited) . pyroaulfurlo  acid  Is  In 
large  excess  of  the  nitric  aold  to  the  extent  thst 


[HNO}]  <.5[H,S|Ot) 

This  Is  the  strong  oleum  region  described  by  Bennett.  The  equilibria 
governing  concentrations  in  this  region  are; 


t> 
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HNO,  + 2H,S,0,  ^ NOJ  + HS,0?  + SHiSO,  (IB) 

H,SO,  + H,S,Ot  HjSO,  + H8,07  (19) 

Equation  18  Is  believed  to  go  completely  to  the  right,  thus  leaving 
Equation  IB  as  the  limiting  equilibrium.  A two-step  reaction  occurs, 
the  nitric  acid  first  completely  dissociating  to  nltronium  ions  via  Equation 
IB,  followed  by  the  establishment  of  equilibrium  among  sulfuric,  pyro- 
sulfuric  acids  and  their  Ionic  counterparts  via  Equation  IB,  The  concen- 
trations are  represented  as  follows; 
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Initial 

After  18 

After  19  (Equilibrium) 

HNO| 

[HNO,]^ 

0 

0 

HtS|Oy 

[H,S,Ot1q-2[HNO,1^ 

[H,S,07]^-2[HN0,1^ 

HIS04 

[H,S04l^ 

[H,S04l^+2lHN0,J^ 

[H,S04l^  +2[HNO,l^- Y 

HS,Of 

0 

(HNO,l^ 

[HNO,]^  + Y 

NO^ 

0 

IHNO.J^ 

[HNO,l^ 

H,S04 

0 

0 

Y 

where  Y ie  the  concentration  of  HS|Ot  formed  via  Equation  19.  The  mase 
action  expreealon  for  Equation  19  it: 

[H,S0t]  [HS,o;i 
^eq=  [H,S04][H,S|0t] 

Subetltution  of  the  tabulated  values  yields; 

Y([HNO,]q  + Y) 

' (tHiSO*]^  +2tHNO,l^  - Y)  (tHiSiOy]^  - 2tHMO,)^  - Y) 

(21) 


If  Kgo  and  initial  concentrations  are  known , Y can  be  found  by  successive 
IteraUons  of  Equation  21 . Bennett  has  suggested  a value  of  1 .0  for  Kgq . 

In  the  second  oleum  region  (intermediate) , the  pyrosulfuric  acid 
is  in  moderate  excess  of  the  nitric  such  that; 

.5(H,S,07]  < [Hi  .]  <[H,S,Ot1 

The  dominant  equilibria  in  effect  here  are  (Ref  19); 

HNOj  + 2H,S07  ^ NO^  + HSjOt  + 2H,S04 

HNOs  + 2H,S04  ^ NOJ  + HjO"^  + 2HSO4 

H,0‘^  + HS|0;'  2H,S04 


(22) 

(23) 

(24) 
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Again  it  is  assumed  that  Equation  22  goes  completely  to  the  right.  Be- 
cause  of  the  limits  imposed  on  the  nitric  acid  concentration , it  must  be 
true  that  the  nitric  acid  remaining  undissociated  after  reacting  with 
H|S|07  (Eq  22)  must  always  be  less  than  the  HS«07_  formed  by  the  same 
reaction.  Thi^s,  there  will  always  be  enough  HS1O7  present  to  soak  up 
all  of  the  H|0  formed  via  equilibrium  Equation  23.  This  continuously 
drives  Equation  23  to  the  right  until  all  the  remaining  nitric  acid  is  com- 
pletely converted  to  nltronium  ions.  Equations  23  and  24  can  be  added 
to  give 


HNOt  HS1O7 


NOt  2HSO4 


In  effect  a two-step  reaction  occurs.  Nitric  acid  reacts  with  the  pyro- 
sulfuric  acid  via  Equation  22 . and  then  the  remaining  undissociated  nitric 
acid  reacts  completely  with  the  HS1O7  produced  from  the  first  reaction. 
Both  reactions  go  to  completion  and  no  equilibria  are  in  effect.  Concen- 
trations are  as  follows: 


Initial 


After  22 


After  25 


H,S,07 

H1SO4 

NoJ 

HS1O7 


(H,S,07l, 

IH,S04]„ 


IHNO,], 


[H»S04l^MH,S,07lp 

.5[HjS,07ljj 

.5[H,S,07l„ 


[HNO,]^  - .5IH,8,07Jj 


[H,S04)^j+  [H»S,07)^ 

[H*S|07lp  - [HNOjljj 
2[HNO,]q  - [H,S|07l|j 


While  Gillespie  maintains  that  Equation  25  is  completely  to  the  right, 
there  is  some  basis  for  the  occurrance  of  the  rever-.e  reaction  (Ref  20) . 
For  this  reason,  the  model  has  been  structured  so  Lhat,  as  an  option. 
Equation  25  can  be  treated  as  an  equilibrium  with 

(HNOjIIHSiOtI 

K„o  = TT <26a) 

^ [N0,](HS07l* 


u 


Initial  concentratlont:  are  aaautned  to  be  those^exlstlng  after  Equation  25 
has  gone  completely  to  right  at  which  time  NO|  and  HSO4  begin  to  rocom- 
blne  via  the  reverse  reaction . Equilibrium  concentrations  are  found  by 
assuming  a value  for  Keq  and  solving  Equation  25a  by  interval  halving . 

In  the  third  oleum  region  (oleum  limited) , nitric  acid  exists  in  ex- 
cess of  pyrosulfuric  acid,  i.o.  [HNO|]  >[H|S|07]  . The  same  equilibria 
are  in  effect  here  as  in  the  intermediate  region , but  now  Equation  25  can 
apply  only  until  the  HS|Of  produced  by  Equation  22  is  completely  dissi- 
pated. There  will  still  be  some  undissociated  nitric  remaining  which  must 
then  undergo  equilibrium  dissociation  via  Equation  23 . In  effect , a three- 
step  reaction  takes  place:  1)  Nitric  acid  -*■  NO|  in  the  presence  of  H1S1O7 , 
2)  Nitric  acid  ■*  NO|  via  Equation  25  until  all  the  HSfOf  is  dissipated, 
and  3)  Equilibrium  dissociation  of  the  remaining  nitric  acid  via  Equation  23 . 
Concentrations  are  as  follows; 


InlUttl 

After  22 

After  25 

After  23  (Bgulhbrium) 

UNO, 

[HNO,]^ 

IHN0,|^-  .5|H,S,0,1^ 

(HNO,)jj-  (H,8,0,]^ 

IHNO,]^  - tH|S,0,l^-  1 

H.S.O, 

0 

0 

0 

H,80, 

4, 

tH,SO,l^+  2k 

NO, 

0 

|H,8,0,lj^  + I 

HS,0i 

0 

8lH,S,0,ly 

0 

0 

HSO,' 

4, 

0 

0 

U'|8,0,1^ 

[H, 8,0,1^  + 2i 

H,0 

0 

n 

0 

t 

22 


((HNO,]^  - *)([H,S04l^  + [H,8,0,]^  - 2*)* 

(27) 

If  Kgq  and  the  initial  concentrations  are  known  > z can  be  found  by  auo- 
ceeelve  Iteration  of  Equation  27 . 

Aa  was  done  for  the  secon^  oleum  region,  provision  has  been  made 
to  consider  recombination  of  NO|  and  HSO4  via  Equation  25.  If  this  option 
la  ohoeen  In  the  modal,  Equations  23  and  2Sa  must  be  eolvad  slmultaneouelv 
to  give  the  equilibrium  concentrations  of  the  various  ionic  species , A pro- 
cedure analogous  to  that  used  In  the  weak  acid  region  to  solve  Equations 
12d  and  12e  le  employed. 

Now,  knowing  Initial  concentrations , and  having  aolved  for  the 
appropriate  concentrations  at  equilibrium.  Bennett's  proposed  rate  ex- 
pression can  be  applied  to  strong  add  nitration  as  follows: 

» (k'tHSOj  + k"[H,S04l  +k"'[H8|OT])  [M][HNO,]pQ 

(26) 

where  Q is  the  fraction  of  nitric  acid  existing  as  nltronlum  ions; 

[M]  Is  the  conoentratlcn  of  nitrotoluene  reactant; 

[N]  Is  the  concentration  of  nitrotoluene  product;  and 
k » k I and  k are  rate  conetants. 
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It  Is  emphRslsed  that  Equation  28  is  continuous  over  the  four  strong 
acid  regions  previously  described . 

The  inHuence  of  temperature  upon  reaction  rate  is  taken  into  account 

I II  m 

by  asauming , as  Bennett  did , that  the  constants  k , k , and  k exhibit 
an  Arrhenius  dependence. 

Thus,  Equation  26  can  be  rewritten  as: 

e X(*T  "T;i)(kj^  (HSOj  + k^  tH.SO*] 

+ kp  |H8,0;])[M1(HN0,1^Q  (29) 


I 

1 

(30)  i 


(k^  (HSOll  +k^  (H,S04l  +kg  [HS,Ot"])[HNO,1^jQ  (31) 

where  k^^  and  are  the  velocity  coefficient  and  activation  energy 
associated  with  the  nitration  reaction: 

M + HNOj  N + H,0 

and  is  the  volume  of  the  acid  phase . 

Subsequent  to  Bennett's  work,  strong  experimental  evidence  was 
uncovered  which  contradicted  a portion  of  his  basic  hypothesis , Melander 
demonstrated  that  proton  departure  from  the  transition  compound  during 
nitration  la  klnetlcally  insignificant  (Ref  12) . Therefore, Bennett's  con- 
tention that  DNT  nitration  in  mixed  acid  is  a termolecular  process  depen- 
dent on  both  the  nitronium  ion  and  the  proton  acceptor  ions,  is  most  likely 
incorrect.  In  addition,  for  two-phase  nitration  systems,  no  maximum  in 


and  Equation  29  can  be  further  simplified, 

~°N  A 1\ 

R«eN“kj,a  R 'T  "Tr/  (m)  y 
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u the  rate  at  92%  H1SO4  haa  been  noted. 

Returning  to  Equations  le  and  17,  it  appears  that  there  is  sufficient 
justification  to  eliminate  the  proton  acceptor  component  concentration  terms 


from  the  expression  for  k|  so  that 

, [NOJ] 

‘'I  Til.s’0,1 

and 

(31a) 

Y «V^k.  [HNO,l 

(31b) 

with  the  dependence  of  nitration  on  the  concentrations  of  SO| , H1SO4 , and 
H|0  considered  only  In  so  far  as  those  concentrations  affect  nitronlum  Ion 
equilibria . 

Although  the  nitration  section  simulation  has  been  modified  to  util- 
ise this  simpler  strong  acid  nitration  model,  the  option  to  use  Bennett's 
mechanism  has  been  retained . 

Far  less  la  known  about  the  kinetic  mechanisms  that  govern  the  oxi- 
dation reactions  which  always  occur  to  some  extent  during  TNT  manufac- 
ture along  with  the  desired  nitration  reactions . With  regard  to  oxidation , 
three  key  assumptions  are  made  in  the  model;  (1)  oxidation  reactions 
occur  exclusively  In  the  acid  phase  and  therefore  compete  with  the  nitra- 
tion reactions , (2)  the  controllinjg  component  In  oxidation  Is  the  same  as 
that  which  governs  nitration  (NO|  ion) , and  (3)  oxidation  reactions  ex- 
hibit biinolecular  kinetics  and  are  thus  subject  to  the  same  reaction  mechan- 
isms as  nitration . Also,  because  of  a lack  of  plant  data  on  oresol  forma- 
tion, oxidation  of  MNT  s not  represented  In  the  model.  This  however,  is 
not  a critical  shortcoi  .g  since  cresole  are  eventually  deotructlvely  oxi- 
dized to  gaseous  products  and  their  contribution  to  overall  yield  loss  can 
be  compensated  for  by  adjustment  of  the  DNT  oxidation  rates.  Equations 
A-42  through  A-95  are  employed  in  the  model  to  represent  oxidation. 

II  Diffusion 

I Diffusion,  as  employed  In  the  nitration  section  model.  Is  the  process 

by  which  the  nitrobody  compcnents  move  bi  tween  the  two  phases  after 
their  acid  phase  saturation  solubility  has  been  reached.  This  phenomenon 
i depends  solely  upon  the  concentration  gradients  of  the  nitrobody  species 


* P 

4 f* 
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aoroBi  each  phaae  and  on  the  fluid  propertlea  of  the  medium  aa  reflected 
by  the  maae  transfer  coefficient . 


1 


u 


At  steady  state , diffusion  of  components  toward  the  acid/organic 
interface  must  be  equal  to  that  away  from  It . 


Thus, 

where 


*1*  at 


*%b 


(32) 

(33) 


(See  Figure  8 for  an  explanation  of  the  nomenclature . ) 

Now,  if  we  assume  that  the  acid  phase  la  always  saturated  with  nitrobody , 


^ X * 


(34) 

(38) 


so  that  by  substitution  of  Equation  36  into  Equation  33  we  get 


X 


(36) 


At  the  interface , the  two  phases  must  bo  in  equilibrium . The  equilibrium 
relationship  assumed  to  be  in  effect  is  that  the  mole  fraction  nitrobody 
component  in  the  organic  phase  must  be  equal  to  the  acid-free  mole  frac- 
tion of  that  component  in  the  acid  phase . 
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Thus , at  the  Interface , 


^nb* 


A * ■ X A* 
^nb  ^eq 


■0  that 


lil  V 

eq  T 


y,*  ■ x,./x,» 


Substitution  of  the  expressions  for  and  into  Equation  32  gives 


but  also  ftrom  Equation  32 , 


so  that 


X • + X. 

i Ha  1 


D.  ■ n “ 

i '(p  <p 


^ (■ 


Now  let 


(40) 


eo  thit 


(47) 


Equation  47  la  true  for  all  organic  oomponenta  except  toluene . For  toluene 
the  reaction  to  MNT  la  ao  faat  that  it  la  aaaumed  to  occur  at  the  interface , 
which  ellmlnatea  the  acid  phaae  dlffualon  of  toluene  and  foroea  AjqI  fo 
lero.  Aleo  the  extremely  rapid  movement  of  toluene  to  the  Inter  face  war- 
ranta  a aeparate  maee  tranafer  ooafflolent . Thue , 


•^Tol 


^Tol 


***Tol 

^nb 


(48) 


Aa  Indicated  prevlouely , the  rate  of  oonverelon  of  toluene  to  MNT 
via  Reactlone  1 and  2 la  eolaly  dependent  upon  the  rate  of  toluene  dlffualon 


Ri  Ri  * D. 


Tol 


R.  R, 


1 


(48) 


Howtver , th«  ratio  of  oMNT  to  mMNT  produced  ia  dependent  on  the  kinetic 
expreaalona: 


Substitution  of  Equation  49  into  Equation  S3  gives 


/_L  -L^ 

- 1 - k,,  •“  R W ■ Tr/ 


(B|  - B|)  / 1 
R \ T 


1 ^ kii  e 


■ Tr^ 


Now,  alnca 


/J_  JLN 

kii  0 ■ R \ T TrV>>>  1, 


we  can  re- 


write Bquation  SB  as 


-B,, 


. C..  . V ( J L\ 

kn  e R \ T Tr  / 


(B.  - B,)/r(^  -5^)  D . (t  -tr) 

n ■ e 


At  this  point,  it  should  also  be  noted  that  while  ohsmioal  reactions 
occur  In  both  the  nitrators  and  separators  due  to  the  aver  present  solubll'* 
lied  organic  components  in  the  acid  phase,  diffusion  is  assumed  to  be 
nonexistent  in  the  separators  because  of  the  drastic  reduction  in  inter- 
face area  which  occurs  there. 
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Bulk  M«m  Tranitftr  and  Solubility 

Tha  remaining  generation  terms  in  the  material  balance  equations 
(Equatlone  5 and  6)  are  those  for  bulk  mass  transfer.  Bulk  mass  trans- 
fer, Is  in  effect  .the  maane  by  which  equilibrium  nitrobody  solubility  is 
enforced  on  tha  acid  phase . At  any  point  in  the  nitration  section , the 
equilibrium  amount  of  nitrobody  in  the  acid  phase  is  given  by 


(SB) 


A 

whara  haa  a value  & pendant  on  tha  local  temperature  a a well  as 

the  local  acid  phase  composition.  For  a given  vassal , nitrobody  is  present 
in  the  acid  phase  by  virtue  of  that  which  entered  with  the  incoming  acid 
phase,  that  which  diffused  in  from  tha  organic  phase,  d^.  and 

that  which  resultad  from  a change  in  solubility  created  by  combination  of 
various  vessel  input  streams  which  are  obviously  not  in  equilibrium , The 
last  of  these  is  considered  the  bulk  mass  transfer,  Thus,  to  enforce  equili- 
brium on  the  acid  phase  leaving  the  vessel,  it  must  be  true  that 


and 


M . « X ^ A_ 
nb  eq  T 


iihb  °i  " ^nb 


nb 


l5nb  °1 


(SB) 


(80) 


I »■ 
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where  represents  the  bulk  movement  of  nitrobody  into  the  acid  phase 
which  must  occur  in  order  to  maintain  equilibrium  solubility . Negative 
values  of  are  reasonable  and  indicate  nitrobody  dropping  out  of  solu- 
tion at  the  new  vessel  conditions . The  individual  component  bulk  maos 
transfer  rates  are  computed  from  the  overall  bulk  nitrobody  transfer  rate 
as  follows; 
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wh«ri  the  index  1 repreeente  etch  nitrobody  component,  The  change  in 
acid  phaae  nitrebody  ccmpoaitlon  due  to  reaction  ia  not  Included  in  the 
bulk  maee  tranefer  term  eince  reaction  only  converta  one  nitrobody  com- 
ponent to  another  and  doea  not  raault  in  a net  molar  change. 

An  expreaelon  for  the  equilibrium  eolubllUy  of  bulk  nitrobody  in 
the  acid  phaae  hai  been  developed  by  eubjeotlng  a limited  quantity  of 
publiahed  aolubility  data  to  a multiple  regreaaion  procedure.'*'  Solubility 
ia  correlated  agalnat  temperature  and  acid  phaae  compoaition  to  give: 

■ 10**  (fli) 


where 


p ■ 3.62  + .01023T  - 3.806  x„  _ + 7.7.18  x 

n|U 


MNT 


- 8.4S3  - 6.262  x 


HNO, 


(62) 


^Laboratory  aolubility  data  for  TNT,  DNT,  and  MNT  In  mixed  acid  at  vari- 
oua  temperaturea  ia  currently  being  generated.  Thle  data,  being  more  ax 
tanalve  than  what  ia  uead  for  the  exleting  aolubility  correlation , will  be 
ragreaaed  to  generate  a new  nitrobody  aolubility  expreeaion  which  ie  ex- 
pected to  significantly  improve  the  nitration  eaction  modal . 
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where 


T le  temperature  in 

X Is  mole  fraction  water  In  acid  phase 
H|  O 

Xmnt  fraction  MNT  in  acid  phase 

la  mole  fraatinn  TNT  in  sold  phase 
TNT  ^ 

mole  fraction  nitric  acid  in  acid  phaae 

nNQ| 

Finally,  an  expreaalon  haa  been  developed  to  repieaent  the  aolu- 
bllity  of  nitric  acid  in  the  organic  phaae.  It  la  the  only  acid  component 
that  shows  any  appieclable  solubility  in  tho  organic  phaae,  The  exprea 
Sion  used  in  the  model  is; 


Bxprasalons  for  all  of  the  terms  in  the  vessel  material  balance  equa> 
tlona  (Equations  6 and  6)  have  now  been  developed  so  that  if  all  vessel 
inputs  are  known  and  can  be  divided  into  acid  and  organic  phase  compon* 
anta,  simultaneous  solution  of  the  vessel  equations  will  give  tho  output 
of  that  vessel  in  terms  of  moles  per  hour  uach  component  in  each  phase, 

The  vessel  equations  can  be  solved  either  for  the  trai  slant  state, 
in  which  esse  Integration  of  the  flrat  order  ordinary  dlfforei  Inl  equa  - 
tions is  required,  or  fur  the  steady  state,  in  which  case  the  derlvutlvu 
terms  are  sol  equal  to  sero,  resulting  In  a system  of  non-llneor,  simul- 
taneous, coupled,  algebraic  equatlone,  Doth  of  thuio  cases  have  been 
solved  and  the  solution  procedures  including  the  computer  programs 
gensrsted  will  bn  described  In  dotull  In  subsequent  iiectlons  of  this  report, 

Formulation  of  the  Ovarell  Nitration  Section  Model 
from  the  Single  Vnaiiel  Modei 

Up  to  this  point,  discussion  has  centered  around  Ihu  details  uf  the 
single  vessel  modal  without  any  regard  to  the  way  In  which  the  compUite 
14-vessel  nitration  train  (B  nltrators  and  6 leparntora)  is  actually  simu- 
lated. "Connootlon"  of  the  individual  vessels  In  the  nitration  sectlun  is 
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made  possible  by  use  of  a physical  separation  model  which  takes  the 
ideally  separated  organic  and  acid  phses  computed  by  the  separator 
vessel  model  and  constructs  the  three  reel  streams  which  actually  flow 
to  other  vessels.  These  Include  the  organic  stream  forward  to  the  next 
nitrator , the  internal  acid  recycle  stream , and  the  external  acid  recycle 
stream  to  the  previous  nitration  stage.  Because  the  separation  process 
is  not  ideal,  a portion  of  the  acid  phase  is  entrained  in  the  oiganic  stream 
leaving  the  separator  and  similarly . a portion  of  the  organic  phase  is  en- 
trained in  the  two  exiting  acid  streams . 

In  developing  the  equations  for  the  physical  separation  model , the 
entrainments  are  defined  as  follows: 

e = molar  fraction  of  the  acid  recycle  stream , both  internal 
and  external,  which  is  entrained  organic  phase. 

e = molar  fraction  of  the  organic  stream  which  is  entrained 
acid  phase . 

The  following  assumptions  are  also  made;  (1)  the  composition  of  ^ 

entrained  material  is  the  same  as  that  of  the  bulk  phase  from  which  it  ori- 
ginated, (2)  the  composition  of  the  internal  and  external  recycle  streams 

are  identical,  and  (3)  the  values  for  e and  e are  constant  for  a given 

OB  ao 

separator . * 

The  separation  process  is  described  schematically  in  Figure  9 . 

It  Is  again  stressed  that  this  physical  separation  model  operates  on  the 
contents  of  a separator  after  the  kinetic  and  mass  transfer  processes  in 
I that  separator  have  transpired . 


^An  option  exists  in  the  model  whereby  e and  e can  be  computed 
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based  on  correlations  that  were  developed  using  linear  multiple  regres- 
sion techniques.  However,  because  of  the  relatively  high  standard  error 

\ of  estimate  associated  with  the  correlation  coefficients,  e and  e are 

i oa  ao 

usually  assigned  their  constant  values . 


o 


Now,  let 


= moles/ hr  of  aclrl  phase  in  Stream  i 

= moles/hr  of  organic  phase  in  Stream  i 

Note  that  the  flow  rates  A and  cp  of  the  acid  and  organic  phases  in  the 
input  stream  ere  known . 

First  consider  the  separation  of  Stream  0 Into  Streams  1 and  2 . From  the 

definition  of  e and  e 
oa  ao 


<Pt 

®oa  ” Ax  + <pi  (64) 

A| 

®oa  Aj  + <P|  (65) 


These  can  be  converted  Into  molar  ratios  as  follows; 


oa 


molar  ratio  of  organic  phase 
to  acid  phase  in  recycle 
streams.  (66) 


1-e  = molar  ratio  of  acid  phase  to 

organic  phase  in  organic 
stream.  (67) 


■1 

'I 

j 

1 


i 


3 


* ► 


< IT 


Material  balances  around  node  A yield: 
<Po  = 9i  + <Pi 
Aq  = A,  + A, 
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Substituting  Equation  66  for  (pj  and  Equation  67  for  At  gives; 


u 


<p  = A.  e'  + <p- 
^o  ‘ oa  ^ 2 


A = Ai  + <p,  e' 
o ‘ oa 


Ellinlnatlng  <pt  and  rearranging , gives 


A.  1-e'  <p  /A 
f - — e ao  ^0^  o 

ar  “ A,  1-a'^^  e'oa 


fraction  of  the  entering  acid 
phase  that  appears  In  the 
acid  recycle  stream , 


Solving  Equation  70  for  <p| , substituting  Equation  72,  and  rearranging 
gives; 


t ^ ‘ — * 1 - f e'  A /o  (73) 

po  <p^  ar  oa  o'  ^o  ' ' 

From  fp^  and  f^^  the  composition  of  each  component  in  the  organic  stream 
can  readily  be  calculated.  Let  A.,  be  the  moles/hr  of  Component  j in  the 
acid  phase  of  Stream  i.  Define  analogously  to  be  the  tnoles/hr  of 

Component  J In  the  organic  phase  of  Stream  1.  For  Stream  2,  these  can 
be  computed  as  follows; 


*0) 


<Pli  “ f <P  J 

po^oj 


stream  1 must  now  be  separated  Into  its  internal  and  external 
recycle  portions,  The  volumes  of  the  entering  phases  (Stream  0) 
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can  be  computed  from  their  oompoBltlons: 

= ft  /hr  of  organic  phase  In  Stream  0 

3 

B ft  /hr  of  acid  phase  in  Stream  0 

Q-p«  ■ Qa«  “ total  ft* /hr  for  Stream  0. 

TO  <pO  AO 

The  total  flow  of  Stream  1 can  be  computed  as; 


Q*  • f Q*  (1-f  ) Q 
^T|  ar^Ao  po'  ^90 


The  fraction  of  Stream  1 that  appears  In  Stream  3 (Internal  acid  recycle) 


'ar  Qao  " %o 


where  Q ■ maximum  ft  /hr  of  internal  recycle  ( a known  constant) 

max 

N • notch  setting  on  a gate  valve  In  the  Internal  recycle 
line  which  Is  modeled  as  an  equal  percentage  valve, 

f(N)  ■ fraction  of  maximum  flow  corresponding  to  the  notch 
setting . 

The  flow  rates  of  the  individual  uomponents  in  each  phase  of  Streams  3 and 
4 can  be  computed  as  follows; 


^r  ^r  ^oj 
^r  (l-'po’  ‘Poj 
^^■'ir’  ^ar\j 
<l-'ir)  *Po1 


u 


Thun  it  is  seen  that  the  composition  of  the  separator  outlet  atreams 
can  be  calculated  by  the  physical  separation  model  from  the  results  of  the 
separator  vessol  model  and  that  thia  in  turn  forma  the  basis  for  connection 
of  the  six  nitration  stages . The  way  in  which  the  stages  are  connected 
will  become  clearer  in  the  next  section  of  the  report,  which  deals  with  the 
implementation  of  the  vessel  and  physical  separation  models  on  a digital 
computer  to  form  the  complete  nitration  section  simulation . 


COMPUTER  IMPLEMENTATION  OF  THE  MATHEMATICAL  MODEL  - 

THE  SIMULATION 

Separate  computer  programs  have  been  prepared  to  aimulate  both 
the  dynamic  and  steady  states  of  the  nitration  process . While  differing 
in  solution  procedure . titere  is  no  difference  in  their  basic  structures; 
both  are  ways  to  solve  the  vessel  equations  to  give  vessel  outputs  for 
known  Inputs,  In  fact,  the  solution  obtained  by  running  the  dynamic 
simulation  long  enough  to  reach  steady  state  is  exactly  that  obtained  by 
executing  the  steady  state  simulation,  within  the  limits  of  numerical  accur- 
acy , A complete  list  of  the  vessel  equations  (the  mathematical  model  of 
any  nitration  section  vessel)  ordered  according  to  a logical  solution  se- 
quence la  given  In  Appendix  A . 

The  Dynamic  Simulation 

The  dynamic  simulation  is  a computer  program  which,  for  a given 
set  of  nitration  section  Independent  variables  Including  (1)  feed  rates 
(Ib/hr) , (2)  composition  of  raw  materials  (wt%) , (3)  nitrator  tempera- 
tures ( c) , and  (4)  recycle  notch  settings,  will  operate  upon  the  vessel 
and  physical  separation  equations  to  give  the  composition  and  flow  rate 
(molea/hr)  of  the  exiting  streams  from  each  of  the  14  process  vessels  as 
functions  of  time.  A numertcal  integration  technique  known  as  Euler's 
Method  is  employed  to  solve  the  differential  equations . Although  this  is 
the  simplest  of  the  many  numerical  integration  algorithms  available,  its 
efficiency  In  solving  the  vessel  equations  is  satisfactory . 

Formulation  of  the  Equations 

In  order  to  solve  the  unsteady  state  vessel  equations  via 
Euler's  Method,  it  is  first  necessary  to  reformulate  the  continuous  differ- 
ential equations  (Equations  6 and  6)  into  their  discrete  numerical  equiva- 
lents . 


as 


For  Component  1 (aMNT)*,  Equation  5 is 


R 


+ D 


n.i 


+ M - A 
,1  n.i 


~dt <“2) 


and  Equation  6 is 


d P. 


n,i 


dt 


(83) 


Equations  82  and  83  are  not  exactly  correct  since  the  units  of  all  terms  on 
the  left  side  are  molea/hr  while  those  on  the  right  are  moles/hr* . To  cor 
rect  for  this,  the  derivatives  are  rewritten  as: 


and 


where 


dt 


d A 


n.i 


dt 


dt 


(V 


d P. 


n.i 


dt 


3 

V if  the  vei eel  volume  (ft  ) 
n 

f.  is  the  volume  fraction  acid  phase  in  the  vessel 

* 3 

p . la  the  density  of  the  acid  phase  (moles/ ft  ) 

” 3 

Pp  is  the  density  of  the  organic  phase  (moles/ft  ) 

X]  la  the  mole  fraction  of  a MNT  in  the  acid  phase 

is  the  mole  fraction  of  a MNT  in  the  organic  phase. 


Therefore  Equations  B2  and  63  can  be  rewritten  as; 


V f p 

A^^R,  - R,  +D,  +Mi  -A,  - ^ 

See  Appendix  A for  component  number  assignments 
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Note  thet  In  Equation!  84  and  85  the  veaael  aubacripta  have  been  dropped 
and  the  aaaumptlon  that  fraction  acid , denalty , and  total  molea  change 
little  with  time  hae  been  made . 


The  numerical  approxlmatlone  for  the  oontlnuoua  derivatlvea  are: 


d At 
“dt“ 

d Pi 
dt 


A t+At  . t 
At  - At 


At 
t+At  _ 


'i 


At 


(86a) 


(86b) 


where  the  auperecript  t repreaenta  a quantity  evaluated  at  the  current 
time  and  t-«-At  Indicated  computation  at  the  end  of  a time  atep . Substitution 
Into  Bquatlona  84  and  85  yields: 


(A? 


VfA  Pi 


+ R,  -R,  +D,  +M,  - IA,‘*^^-At*l 


(87) 


lPi°-Dt  -Mt  -Pil*» 


VU-tA>  Pp 
P^  At 


lPt‘^^^-Pi‘) 


(88) 


Solving  for  and  Pi**^^  wo  get; 


Ai^'^^'»A,^+  (At°+  (Ri  -R,)  +D,  ■•■M,  - Ajl^At^  (89) 

“ Pt*  + (Pi°  - Dt  - Ml  - Pt  Atp  (90) 
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where 


(ditnenslonleee) 


At 

At.  = _JE 

* '"a  Pa 


P«  At 


(dimenelonlaee) 


Bquetlona  enalogoua  to  89  and  80  may  now  be  written  for  all  com-  ^ 

ponenta  In  the  organic  phaae  (except  P|  which  la  given  by  Equation  63) 
and  for  all  oomponenta  In  the  acid  phaae.  and  then  aolved  via  the  Euler  ^ 

algorithm . Thla  procedure  waa  Incorporated  Into  the  original  computer  i 

program  which  waa  written  to  solve  the  unateady  atate  aquatlona . Al-  i 

though  the  procedure  la  theoretically  and  numerically  conalatent , it  waa 
found  that  baoauae  of  the  large  numerical  valuaa  of  the  dlffualon  terma , | 

a vary  email  Integration  atep  alae  (At)  waa  required  In  order  to  obtain 
numerical  atabillty  and  generate  convergent  eolutlona . Thla  waa  not  I 

acceptabla , alnce  uae  of  an  extremely  email  At  reaulted  In  Intolerably  * 

long  computer  rune . To  avoid  thla  difficulty,  an  alternate  aolutlon  pro- 
cedure wae  devlaed  baeed  on  the  fact  that  the  accumulation  of  nltrobody  ^ 

In  the  acid  phaae  la  email  compared  to  that  in  the  organic  phaae.  Thla  , 

allowa  the  Integration  of  the  nltrobody  oomponenta  to  be  carried  out  over 
the  combined  total  amount  of  nltrobody  In  both  phaaea,  while  the  Integra-  | 

tlon  of  the  acid  componenta  In  the  acid  phaae  la  conducted  aa  before . The 
concentration  of  the  Individual  nltrobody  oomponenta  In  the  add  phaae  ^ 

can  than  be  obtained  by  aolvlng  the  appropriate  ateady  atate  equatlona 
after  the  Integration  haa  bean  carried  out.  In  effect,  It  la  aaaumed  that  | 

the  nltrobody  In  the  acid  phaae  la  In  equilibrium  with  that  In  the  organic  \ 

phaae.  Implementation  of  the  procedure  deacrlbed  above  la  carried  out 
by  flrat  adding  Equatlona  84  and  86  to  give:  | 


I Pj°  + Ai°  + R,  - R, 

I 

1 


(Pj  + Al) 


dP, 

dt 


A.J.  dt 
(91) 


(Note  that  the  troublaaome  maaa  tranafer  tarrna  have  dropped  out  of 
Equation  91 . ) 
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Equation  97  la  called  the  normal  or  forward  difference  form  of  the  dlecrete 
equation , where  the  term  in  brackets  la  evaluated  at  the  beginning  of  the 
Interval.  However,  the  efficiency  and  stability  of  the  numerical  Integra- 
tion procedure  is  enhanced  if  the  backward  difference  form  la  used , in 
which  case  the  bracketed  term  is  evaluated  at  the  end  of  the  Interval . In 
this  case: 


- T|*  + Atp  tPi°  + A»°  + Ri  - R,  - T, 


Solving  for  , we  got: 


T»*+At-  iPi®  ♦Ai^+Ri  -R,l* 


M a a _ 

Ti  ■ 


But  (from  Appendix  A) , 


R,  ■ Ai  (Q|6  Oil  Y ) ■ Ai  P, 
From  the  steady  stats  equation 

Ai®  + Rj  - R|  ♦ Di  M|  ■ A|  ■ 0 

we  can  determine  a value  for  At  as  follows: 


(100) 


(101) 


A.  •A.'tR.  -0,8A,  -0„yA,  * A„J,V  * A„^ ) 


NB  nb 


(1  0,a  + 0,,Y  + T*“  - At  * At®  + Ri  + Pi 


4 • 


Substitution  of  Equations  100  and  107  into  Equation  99  gives; 

t/  4 Atp  IA,°  + ♦ R|  - c,  (A,°  + n,  + C 


t>At 


1 + 


(108) 


RMrrtngamtnt  of  Equation  106  yialda: 

t 


- t+At  

■ I ♦ AtpY 


(109) 


wharo 


Pi°  + (Ri  + Ai°)  (1  - 


F|  +C. 


(110) 


Y * 


(ttVc^)  °x 


(111) 


t<«>At 


Th«  valuf  of  A|  can  now  b«  computad  from  Equation  107  uaing  tha 

vaiua  of  T datarminad  in  Equation  109.  A..  . M^.  . y . and  0 in  Equation  107 

t no  t no  ' 

ara  oomputad  uaing  th>  valuaa  of  A|  and  Pi  rat)iar  tlian  A|  and 

and  Pi^^^V  but  tliia  haa  not  poaad  any  problama  in  tha  axaoutlon  of  tha  pro- 
gram . Bquationa  aimllar  to  109  can  now  ba  writtan  to  detarmina 


through  with  A|**^Shrough  A^^^^  being  oalouiatad  from 


aquationa  analogous  to  107, 


For  the  acid  componanta  in  the  acid  phase,  forward  difference  inte- 
gration is  amployad . 


• a/  + IA,°  + P,°  - P,  - R,  - - R.  - R4  - R|  “ R| 


- 13  Rt  - la.P  R,  - 9 R,,  - 12.6  R,q  + (-33  •*  .67  HP)  R, 


-A,l  At^ 


(112) 
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AL*  ■ At^  + [A|f  - Ai  J + Ri  * R,  + R,  + R4  + R|  + R| 


+ 16  R,  + 16.5  R,  ♦ 12  R,,  16.5  R,q  - (.33  RP  + .87)  R, 


- ALl*  At 


(113) 


. A.t^At  * t+At  . t+At 

whart  AL  "An  ~ At4 


(114) 


If  AL  < 0 


All  • 0 
A|  4 * AL 


and  If  AL  > 0 


Alt  • AL 
A, 4 ■ 0 


u A,/  4 (A,,°  + (11  R,  M4.8  R,  4 8 R,,  ♦ 14.8  R»q) 


- A,|l* At 


(118) 


(A,,°  +At4‘'-(11RT^14.8R,  +9Ri, 
M4.8R,o)  (1-Fd)  ■ 
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Solution  Procodurt  for  Dynomtc  Simulation 


Tha  primary  saquenca  of  oparatlona  which  muat  ba  axaoutad  In 
order  to  obtain  the  unataady  atata  aolutlon  of  the  nitration  aaction  equa- 
tlcna  la  aa  followa; 

Stop  1;  Enter  the  Initial  ataady  atata  of  tha  prooaaa,  Thla 
inoludoi  tha  oompoaitlcin  and  flow  rataa  of  tha  idaally  aaparatad  aoid  and 
organic  phaaai  from  each  of  tha  14  prooaaa  vaaaala,  aa  wall  aa  tha  oor- 
raapondlng  faad  rataa , oompoaitiona , nltrator  tamparaturaa , and  notoh 
aottingi . 

Stop  2;  Prooaaa  tha  initial  oondltlona  through  tha  phyaioal 
aaparation  aquaUona  for  all  alx  aaparatora  to  gat  Initial  oondltlona  on  tha 
raoyola  atreamo. 

Stop  3i  Road  in  tha  parturbationa  to  tha  indapandant  varl- 
ablaa  whioh  oonatltuta  tha  davlationa  from  ataady  atata . 

Stop  4i  Bogin  axaoution  of  tha  voaaal  oquaUona  for  Vaatal  1 
(Nitrator  lA)  by  making  up  tha  tarma  required  to  oomputa  tha  generation 
axpraaalona  Dj,  M|^,  and  Rj. 

Stop  Si  Intagrata  ualng  Bquationa  107  to  117  to  gat  tha  valuaa 
of  Pi  and  Ai  at  tha  and  of  tha  flrat  time  atop. 

Stop  6i  Store  thoaa  valuaa  aa  tha  currant  outputa  from  Nltrator 
lA,  and  uaa  them  to  make  up  tha  inputa  to  Nltrator  IB, 

Stop  7;  Itapeni  Stapa  4 and  S for  Nltrator  IB. 

Stop  8:  Store  there  valuaa  aatha  currant  outputa  from  Nltrator 
IB  and  uaa  them  to  maka  up  tha  Inputa  to  Separator  1 , 

Stop  8:  Repeat  Stopa  4 and  0 for  Separator  1 , 

Stop  10:  Store  thaao  valuaa  aa  tha  currant  Idaally  aaparatad 
outputa  from  Separator  1 . 

Stop  11:  Prooaaa  tha  currant  valuaa  of  the  Idaally  aaparatad 
output  from  Separator  1 through  tha  phyaioal  reparation  aquatlona  to  give 
the  currant  oompoaltlon  and  flow  of  the  three  atraama  leaving  eeparator  1 , 


T 
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step  12:  Store  these  values  for  use  in  making  up  inputs  to 
other  vessels. 

Step  13:  Begin  execution  of  vessel  equations  for  Nitrator  2 and 
continue  the  procedure  for  all  remaining  vessels  through  Separator  6. 

At  this  point,  one  complete  pass  will  have  been  made  across  the  en- 
tire nitration  reactor  train  for  the  first  time  step  with  the  values  of  each 
vessel's  output  used  to  make  up  the  necessary  inputs  for  execution  during 
the  second  time  step , This  procedure  is  now  repeated  for  n time  steps  as 
desired.  The  program  currently  employs  a time  step  of  ,01  hours  (36 
seconds)  . Thus,  to  simulate  12  hours  of  process  operation,  1200  time 
steps  .lUSt  be  executed.  Typically,  the  execution  time  required  for  such 
a case  is  about  120  seconds  on  the  CDC  6500  computer  giving  an  average 
of  0.1  socond  of  machine  time  per  time  step. 

Program  Description 

The  nitration  section  dynamic  simulation  consists  of  a main  program , 
four  subroutine  subprograms,  and  one  function  subprogram,  Appendix  B 
presents  a complete  program  listing:  Appendix  C shows  the  program  flow 
chart:  Appendix  D gives  the  program  nomenclature:  and  Appendix  E pre- 
sents the  input  data  format. 

All  programs  are  written  in  FORTRAN  IV  for  compilation  on  the  CDC 
FORTRAN  extended  (FTN)  compiler,  version  3.0  and  require  65K  (octal) 
words  on  the  CDC  6500  computer  for  complete  execution.  Word  length  on 
the  CDC  6500  is  60  bits.  Since  there  are  no  machine-dependent  routines, 
the  program  should  execute  on  any  comparable  computer  with  only  minor 
modifications.  It  should  be  noted  that  a number  of  library  subroutines 
(NANCY,  NUCHAR,  and  NANCYM)  are  used  to  plot  desired  output. 

The  MAIN  program  reads  in  all  data  and  parameter  values , moni- 
tors and  controls  the  sequential  execution  of  the  vessel  and  physical  sep- 
aration equations,  and  keeps  track  of  all  vessel  inputs  and  outputs.  The 
program  also  computes  the  values  of  a number  of  process  performance  indi- 
cators from  the  current  values  the  vessel  outflows  at  various  times  during 
the  run  and  plots  them  as  a function  of  time.  These  performance  indica- 
tors constitute  the  primary  output  of  the  dynamic  simulation  and  include: 
concentration  of  DNT  in  the  crude  TNT  (wt  %),  flow  rate  of  crude  TNT 
(Ib/hr);  concentration  of  nitrobody  in  the  spent  acid  (wt  %) , the  fraction 
DNT  in  the  spent  ncld  nitrobody  (wt  ’h) , concentration  of  water  in  the  spent 
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acid  (vrt  %) , flow  rate  of  spent  acid  (Ib/hr) , total  carbon  and  nitrogen 
oxides  in  the  off-gas  (Ib-moles/hr) , heat  loads  in  each  nitrator  (BTU/hr) , 
concentration  of  nitric  acid  in  the  external  recycle  from  each  separator 
(wt  % actual  nitric) , ratio  of  nitric  to  sulfuric  acid  in  each  nitrator , the 
amount  of  nitration  occurring  in  each  nitrator , and  the  specific  gravities 
of  the  external  acid  recycle  and  organic  streams  leaving  each  separator, 

The  values  of  e^^  and  e^^  for  each  separator  given  by  Equations  66  and  67 

are  also  computed  in  the  main  program  for  later  use  in  the  physical  separa- 
tion equatlonB . 

Subroutine  INPUT  is  callod  by  the  main  program  to  read  in  the  start- 
ing values  of  the  independent  variables  including  raw  material  feed  rates 
and  compositions,  nitrator  temperatures,  and  internal  recycle  notch  setting, 
as  well  sfl  to  read  in  the  magnitude  and  time  of  the  desired  perturbations 
to  these  quantities.  Since  feed  rate  data  is  entered  in  Ib/hr,  the  subrou- 
tine makes  the  appropriate  conversions  to  molar  units.  Also,  the  exponen- 
tial tactors  in  the  rate  equations  (e.g. , 

K,.xp  , 

are  precomputed  in  INPUT  in  order  to  increase  the  efficiency  with  which 
the  vessel  equations  are  executed . 

The  heart  of  the  dynamic  simulation  is  subroutine  SUB , Here,  all 
of  the  vessel  and  physical  separation  equations  are  executed  and  integrated 
each  time  the  subroutine  is  called  by  the  main  program  as  it  steps  through 
the  reactor  train.  In  effect,  the  main  program  provides  the  necessary  input 
data  to  SUB  (l.e.  the  main  program  tells  SUB  whether  it  is  a nitrator  or 
s»^  irator  and  then  which  one  it  is  in  the  reactor  train;  it  provides  vessel 
cc  ic'tants  for  that  vessel  such  as  volume  and  notch  setting;  it  provides  the 
in)  it  and  output  phase  compositions  for  the  vessel  from  the  previous  time 
stfp).  SUB  then  generates  the  output  of  the  vessel  and  returns,  SUB 
also  computes  a temperature  rise  for  each  separator  based  on  the  adiabatic 
heat  generated  by  the  small  amount  of  reaction  which  takes  place  there , 

Subroutine  SUBl  computus  the  change  in  the  calculated  values  of 
the  phase  compositions  for  each  vessel  from  one  iteration  to  the  next.  An 
iteration  is  defined  as  the  execution  of  m time  steps  where  m is  specified 
at  the  beginning  of  the  run.  SUBl  also  flags  that  component  in  each  ves- 
sel which  shows  the  maximum  change  between  iterations . 
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Subroutine  SUB2  executes  the  physical  separation  equations  for 
each  separator  using  the  steady  state  phase  compositions  of  the  Input 
data.  SUB2  Is  called  by  the  main  program  prior  to  execution  of  the  first 
time  step  only  In  order  to  Initialize  the  recycle  flows . The  physical  separ- 
ation equations  are  thereasfter  executed  In  subroutine  SUB  as  the  run  pro- 
gresses. 

The  final  subprogram  In  the  dynamic  simulation  is  function  ENTH , 

This  function  Is  used  In  conjunction  with  the  nltrator  heat  balance  calcu- 
lations executed  by  the  main  program . Its  purpose  is  to  caloulste  the  en- 
thalpy of  any  specified  flow  stream  in  the  process  by  interpolating  enthalpy- 
composition  data  which  exits  as  a data  table  within  the  function.  The  en- 
thalpies of  the  input  and  output  streams  are  combined  with  the  heats  of 
reaction  end  heats  of  mixing  to  compute  each  nltrator 's  heat  duty . 

The  Steady  State  Simulation 

The  steady  state  simulation  is  a computer  program  which,  for  a given 
set  of  nitration  section  independent  variables,  will  execute  the  vessel  and 
physical  separation  equations  to  give  the  steady  state  values  of  the  compo- 
sition and  flow  rates  (moles/hr)  of  the  exiting  streams  from  each  of  the  14 
process  vessels.  In  addition,  the  values  of  a number  of  process  perform- 
ance functions  are  computed . 

At  steady  state,  the  accumulation  terms  In  Equations  A-57  to  A-B2 
go  to  zero  so  that  differential  equation  system  of  the  dynamic  case  becomes 
one  of  simultaneous  algebraic  equations . Because  of  the  implicit  nature  of 
\ I these  equations,  and  also  because  of  the  presence  of  numerous  recycle 

I I streams,  a multilevel  (nested)  iteration  procedure  is  used  along  with  a 

' ! simultaneous  equation  solving  technique  to  arrive  at  the  final  solution . j 

^ I ‘ 

1.  ' S 

1 Formulation  of  Equations 

! 1 i 

f » For  the  nltrobody  components  in  the  acid  phase,  the  following  sys- 

' tern  of  steady  state  equations  can  be  written;  | 

f.  ■;  ] 

For  component  1:  Ai°  - R3  + D,  + Mj  - A|  = 0 (IIB)  | 

i 

■I 

Ai  + R,  - Dj  - Ml  = Ai^'  (119) 


I 
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Substitution  of  Equations  IB,  20,  21,  and  38  from  Appendix  A into  Uqua 
tion  119  gives 


/Pi  Ai\ 

A,  + (0,8  + 0,sY)  A, -n(p— r^J  - 

N*^nh  nh' 


nb  ^nb'  ^nb 


Simplifying , 


A|  (1  + 0,8  + 0,g  Y)  - c7|L  - A,  - Cj^  Pi  - A 

' nb  nb  ^ 


where:  ■ q 


C . ■ for  M . < 0 

A nb 


®A  ■ “ «nb  i " 


Cj^  ■ 0 for  < 0 


Further  rearrangement  of  Equation  121  yields: 


(1  ♦ 0,8  + Q,g  8 + "°A^'^‘AP PjPi“A, 


Similarly,  for  component  2: 


A|  “ R4  + D|  + M|  “ A|  ■ 0 
Aj  ♦ R,  - D,  - M,  - A,° 


Appropriate  aubatituUon  frotn  Appendix  A ylolda: 


A|  + (G48  + G4S  Y ) A, 


Aid  + O46  + 04g  Y -C  ) ‘ )P,  -A,° 

nb  ^ ^ nb  « 

For  component  3: 

A|  + R|  “ R|  - R|  a"  ^iA  " ^tA  **'  “ A|  • 0 


A|  - R|  R,  -f  - R|A  + R|A  “ D|  - M|  ■ A|* 


A|  - (0|fi  ■♦•  0,g  Y ) Ai  + 0|YA|  + GtyA|  + 


VOif.A, 

A,  + Q,yA, 


+ 0,Y  +0,  Y + -Qi^  - + 0,  Y + 


A|  + (Gifl  + 0,g  Y)  A| 


For  component  4; 


A4  ^ R4  ■ R|  - Rtm  ” RiM  " + D4  + M4  - A4  « 0 


(1311 


A|  - R4  + R«  + RiM  ^'•M  - D4  - M4  » A4 


VO|f^A, 

A4  - (O48  041 Y)  A,  ♦ 0,yA4  + — qT^  A4  + 0,yA4 

''A 


/ V0,£^A,  C,  \ 

A|  ( 1 + OiY  + 0,Y  + Ql  ■ ♦ * AT  ■ ®A  ) 

\ '^A  nb  ' 

- ■*■  0 **«  " ^ ^ 


For  oomponant  S: 

A|®  + R|  - Rua  + D,  ♦ M,  - a,  ■ 0 

A|  - R|  ■♦■  Ri4A  - D|  - M|  ■ A|° 

VOnf.Ait  « Q 

A|  - OiYA|  + —77® A,  - P,  + at  A|  - 


A,  (1 


- Cj^Pi  ■ A| 


VOi.f  A„  c, 


“ A|  “ G|YAt 


*=rV- 


For  component  6; 


A|**  + 3 Rij  + D|  + Mg  - Ag  ® 0 

Ag  " 3 Rig  “ Dg  - M|  ■ Ag 

SVOigf.Ait  Cl  Cl 

A. gr (A.  *A,»  - _P,  *—  A,  -C^A. 

-0„».-A.“ 


(147 

(148 

(149 


For  the  nltrobody  epeolee  In  the  orginlo  phase,  the  steady  state 
veaael  aquations  can  be  written  as  follows: 

Pi°  + Ri  - Ml  - Di  - Pi  ■ 0 (15: 

P,  - Ri  + Ml  + Di  «Pi°  (15! 

C C 0 

Pi  - Dig  (1  - Oi)  + C.Ai  + CjjPi  + Pi  A,  ■ Pi  (1B; 


f 


■o  that , 


P,  -cIp,,  (1-0,)  +A,  +Pi(Cj^+|^)=  P,° 

nb  nb 


(C^-^  )A|  + (l4Cj^  + ^)  P,-P,‘ 


♦ C,  P„  (1  - 0|) 


For  oomponent  2: 


Pi  + R,  - M,  - D|  - P|  - 0 


(157) 


P,  - II,  M|  D,  ■ P,‘ 


P,  “ cl  P,,  0,  ♦ C^A,  + C^^P,  + 1^  P,  - A,  • P,°  (159) 

nb  nb 


(C^  - ^Jj^)  A.  (1  +Cj^  + |^)  P,  -P|°  + cl  P„  0, 


For  component  3: 


P,  - M,  - D,  - P,  ■ 0 


P,  + M,  + D,  • P, 


P,  + C.A,  4.  C-P,  + ^ p,  - A,  - P,° 
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(C^  - ^ ) A,  + (1  ♦ Cjj+  P,  - p,” 

nb  nb 

For  oomponont  4; 

p4°  - M4  - D4  - P4  - 0 
?4  ♦ M4  + D4  - P4® 

(C^  - ^ ) A4  ♦ (1  + Cj^  + ^ ) P4  . P4 

nb  nb 

Similarly,  for  oomponanta  S through  B: 

C|iA|  ^ C)|P|  « P|° 

C11A4  ♦ C||P|  • P|® 

C|,At  ^-C„P,  -P,® 

C|iA|  + C||p4  ■ P|® 


whara: 


o„.Ca 


^ and  0„  ■ (It  Cj,  ♦ 
nb  nb 


Pi?  " Dii  Pii  ■ 0 
Pi.  +01,  -P,S 
Pi.  (1  + ci)  -P»f 
Pi.  “Pi?/(1  + cl) 


Now,  for  toluana; 


For  the  Bold  ■pedes  In  the  add  phase,  the  steady  state  vessel  equa- 
tions are  written  as  follows; 

Af  P|  ~ P|  **  1^11  ~ P|  ~ P4  ~ P|  ” “ 13Ry  ~ 12.SR| 


- 9R||  - 12.5R|0  (.33  4 ,67  RP]  R|  - A,  ■ 0 (176) 


Making  the  appropriate  substitutions  from  Appendix  A gives; 

A|  + P|  4 cl  Pii  4 (0,6+0|sY)  A|  4 (Q48  + OisYlAi  4 G,yA| 


4 0«YA4  4 130,y(A|  4 A4)  4 12.6  ^ 


VO|f. 


[HNO| 


(A,  4A4I 


4 12.60,y(A4  4 A4)  4 9 (NOi  ]J  (A,  4 A,) 


- (.33  4 .67RP)PnC,  ■ A|°  4P,° 


' / VO,f. 

A,  4 p,  + C|P,4  4F,A,  - F4A|  4 f 0,Y  4 13O4Y  4 12.5-'—'-  [HNOjl 

4 12.6  0|yJ  A|  4 ^0,Y  4 130tY  + 12.5--^---  lHNO|]  4 12.60|yJ  A4 
VO,,f 

® ■’q^  A|  4 9 — — (Nol’]  A|  - (.334  .67RP)FjjQ| 


A,°4P,o 


Also, 


P| 


(178) 


■( 


1.61  Ai,P^N 

AjA||  / 


A| 


All  ♦ Au  “ C|p  “ AgQ^  * 0 

Ajq^  + (1  - PjjH  * "All  + Ai4 

All  + C|p  “ All  ■ J 
All  • C|p  « A|f 


(179) 

(180) 

(181) 

(182) 


AiJ  ” AiJ  + Ri  ♦ R|  + R|  ♦ R4  + R|  + J . ♦ IBRf  ♦ 16.8R|  ^ 12Rn 


♦ 18.8R,o  “ (-SSRp  f ,87)R|  • Aj^  ■ 0 


(183) 


Al  - Cl  Pii  - F,A,  - F4A1  - f 0,v  + 180i  18.8 


-(0, 


VOif, 


[HNO,  1 


f le.soiY 


0|Y  ISOfY  -*-18.6 


[HNO,] 


4-  16 


r')A.-(c 

.8G,y^  A«  - 12  [NO|  ] ^ - 12  (NOjl  A, 


♦ (.33RP  + .87)  FjjCi  - Ai?  + A14 


(184) 


Equations  122,  126,  130,  134,  138,  142,  148,  150,  156,  160,  164,  167,  168, 
188,  170,  171,  176,  177,  178,  178.  182,  and  184  are  linear  in  the  following 
unknowns:  Ai  through  A| , A^,^  .An,  A^  , and  Pi  through  Pn  . This 


constltutat  a system  of  22  equations  and  22  unknowns  which  Is  solvable  ^ 

provided  the  cueffiolents  of  all  the  unknown  terms  can  be  determined . 

Solution  Procedure  for  the  Steady  State  Simulation 

Execution  of  the  steady  state  simulation  consists  mainly  of  the  simul- 
taneous solution  of  the  22  vessel  equations  for  each  of  the  14  veosals  In  the 
nitration  train . Although  these  equations  could  be  solved  simultaneously 
via  Gauss  reduction , It  was  decided  that  because  of  the  sparclty  of  the  co- 
efficient matrix,  a straightforward  direct  substitution  approach  would  be 
more  efficient.  This  will  become  clear  as  the  step-by-step  solution  pro- 
cedure is  outlined  below.  The  physical  separation  equations  must  also 
be  employed  hare  as  in  the  dynamlo  case  to  convert  computed  pure  phase 
outputs  to  real  stream  outputs . As  mentioned  previously , a multilevel 
or  nested  Iterative  procedure  la  used  because  of  (1)  the  implicit  nature 
of  the  vessel  eqttatlona  (the  coefficients  of  the  linear  unknowns  are  func- 
tions of  these  unknowns) , and  (2)  the  numorous  recycle  streama  in  the 
procsas. 

The  primary  sequence  of  operations  which  must  be  executed  In 
order  to  obtain  the  steady  solution  is  as  follows; 

Step  1;  Enter  an  assumed  solution  for  the  phase  outputs 
from  each  of  the  14  process  vessels. 

Step  2:  Enter  the  values  of  the  process's  Independent 
variables  for  which  the  steady  state  solution  is  desired . 

Step  3:  Using  the  assumed  solution  for  the  vessel  phase 
outputs,  oomputfl  th<>  flow  rates  and  compositions  of  all  recycle  streams 
in  the  process  via  tl  j physical  separation  equations. 

Step  4:  Begin  the  Iterative  solution  of  the  vessel  equations 
for  the  first  vessel  (Nltrator  #1A) . Start  by  computing  the  values  of  all 
ooeffloients  of  the  22  x 22  vessel  equation  matrix  from  the  assumed  solu- 
tion for  that  vessel . Also,  make  up  the  correct  phase  Inputs  to  Nltrator  lA . 

Step  5:  Solve  Equation  175  for  Pk  . 

Step  6:  Solve  Equations  122  and  156  simultaneously  for 
values  of  Aj  and  Pt  • 
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step  7;  Solve  Equatloni  126  and  160  aimulteneoualy  for 
valuea  of  A|  and  P| . 


Step  8:  Solve  Bquationa  130  and  164  elmultaneousty  for 
valuee  of  A|  and  P| . 

Step  8:  Solve  Bquationa  134  and  167  elmultaneoualy  for 
valuaa  of  A4  and  P« . 

Step  10:  Solve  Bquationa  136  and  168  elmultaneoualy  for 
valuee  of  A|  and  P| . 

Step  11:  Solve  Bquationa  142  and  160  elmultaneoualy  for 
valuee  of  A|  end  P| . 

Step  12:  Solve  Bquationa  146  and  170  aimultaneoualy  for 
valuee  of  Af  and  Pf . 

Step  13:  Solve  Bquationa  180  aud  171  aimultaneoualy  for 
valuaa  of  A|  and  P| . 

Step  14:  Solve  Equations  176  and  178  aimultaneoualy  for 
valuee  of  A|  and  P| . 

Step  18;  Solve  Equation  182  for  Ah  . 

Step  16;  Solve  Equation  184  for  A^. 

Step  17;  Uae  the  logic  of  Bquationa  A-70  and  A-71  to  determ- 
ine valuaa  for  Ah  and  Ah  • 

Step  18:  If  computed  values  and  eoaumed  values  differ  by 
more  than  a predetermined  amount,  repeat  Steps  8 through  17  using  newly 
computed  values  for  the  Aj'a  and  Pj^'s. 

Step  19;  Repeat  Step  18  until  oonvergenoa  le  achieved  on 

Nltrstor  lA. 

Step  20;  Using  the  converged  solution  for  Nltrstor  lA,  and 
the  previously  assumed  recycle  stream  flows,  make  up  the  correct  phase 
Inputs  to  Nltrstor  IB  and  then  repeat  Step  4 for  Nltratnr  IB, 
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Step  21:  Repeat  Steps  5 through  17  for  Nltrator  IB.  J 

Step  22:  Repeat  Step  18  for  Nltrator  IB. 

Step  23;  Repeat  Step  22  until  convergence  is  achieved  on 

Nltrator  IB . 

Step  24;  Using  the  converged  solution  for  Nltrator  IB  as 
the  input  to  Separator  1 and  the  initially  assumed  phase  outputs,  -jcute 
Steps  4 through  19  for  Separator  1. 

Step  25:  Process  the  converged  outputs  for  Separsior  1 
through  the  physical  separation  equations  to  get  the  flow  rates  and  com- 
positions of  the  recycle  streams  leaving  that  separator. 

Step  26;  Compare  the  calculated  recycle  stream  between 
Separator  1 and  Nltrator  IB  (Internal  recycle)  to  that  assumed  when  the 
Inputs  to  Nltrator  IB  were  made  up.  If  they  differ  by  a preset  amount, 
use  the  newly  calculated  recycle  stream,  and  the  previously  converged 
values  of  the  outputs  from  Nitrators  lA  and  IB  to  repeat  Steps  4 through 
19  for  Nltrator  IB. 

Step  27:  Repeat  Steps  24  through  26  until  convergence  on 
the  Internal  recycle  stream  is  achieved. 

Step  26:  Using  the  converged  Separator  1 outputs  and  the 
originally  assumed  Nltrator  2 output  and  Separator  3 recycle  stream 
values,  execute  Steps  4 through  19  for  Nltrator  2. 

Step  29:  Repeat  Steps  24  through  27  for  Separator  2. 

Stop  30:  Repeat  Stops  28  and  29  for  Stager^  3 through  6 until 
convergence  on  all  internal  recycle  streams  is  achieved.  Note,  although  a 
converged  solution  for  each  stage  will  have  been  obtained  at  this  point, 
no  attemot  at  linking  the  six  stages  together  through  their  external  recycle 
streams  has  yet  been  made.  This  must  now  be  done  in  order  to  obtain  the 
complete  steady  state  solution  for  the  14-vessel  reactor  train. 

Step  31:  Begin  with  the  external  recycle  between  Separator  6 
and  Nltrator  5.  Using  the  previously  converged  outputs  from  Separators  4, 

5 , and  6 as  Nltrator  5 inputs  and  the  previously  converged  Nitrator  5 out- 
put, execute  Steps  4 through  19  for  Nltrator  5. 
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step  32;  Repeat  Steps  24  through  27  for  Separator  S. 


Step  33:  Execute  Steps  4 through  19  for  Nitrator  6 using  the 
newly  converged  Separator  S output  and  previously  converged  Separator 
6 and  Nitrator  6 outputs . 

Step  34;  Repeat  Steps  24  through  27  for  Separator  6 . 

Step  35:  Compare  the  newly  calculated  external  acid  recycle 
stream  from  Separator  6 to  that  employed  in  Step  31.  If  they  differ  by  a 
preset  amount,  use  the  newly  calculated  value  as  well  as  the  most  recently 
converged  values  for  the  other  required  flows  to  repeat  Steps  31  through  34. 

Step  36:  Repeat  Step  35  until  convergence  on  the  external 
recycle  from  Separator  6 is  achieved . 

Step  37:  Continue  with  this  stagewise  iteration  procedure 
until  the  remaining  external  recycle  streams  (and  thus  the  entire  nitration 
reactor  train)  have  b^-jn  converged, 

Typically , a steady  state  run  requires  approximately  5000  iterations 
of  the  vessel  equations  with  convergence  limits  set  at  .05  mole/hr  for  over* 
all  vessel  composition,  .02  moles/hr  for  components  in  the  Internal  recycle, 
and  .02  moles/hr  for  components  in  the  external  recycle.  Under  these  con- 
ditions, program  execution  time  is  approximately  30  seconds  on  the  CDC 
6600.  Obviously,  if  tighter  convergence  is  desired,  run  time  increases 
significantly . 

Program  Description 

The  nitration  section  steady  state  simulation  consists  of  a main  pro- 
gram, two  subroutine  subprograms,  and  one  function  subprogram.  A com- 
plete steady  state  run  Including  the  program  listing  is  presented  in  Appen- 
dix F.  Appendix  G shows  the  program  flow  chart;  Appendix  H gives  the 
program  nomenclature;  and  Appendix  I illustrates  the  input  data  format. 

All  programs  are  written  in  FORTRAN  IV  for  compilation  on  the  CDC 
FORTRAN  extended  (FTN)  compiler,  version  3.0,  and  will  compile  and 
execute  in  55, COO  (octal)  words  on  the  CDC  8500.  Since  there  are  no  ma- 
chine-dependent routines , the  program  should  execute  on  any  comparable 
computer  with  only  minor  modifications.  A brief  description  of  the  main 
program  and  the  threa  subprograms  will  now  be  given . 
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The  main  program  (TNTSIM)  reads  in  all  data  and  parameter  values 
which  form  the  fixed  data  base  for  the  steady  state  simulation . Included  in 
this  data  base  is  an  initiai  guess  at  the  steady  state  solution.  Also,  the 
program  is  responsible  for  reading  the  particular  set  of  values  for  the  in- 
dependent process  input  variables  for  which  the  steady  state  output  is  de- 
sired. All  feed  streams  which  are  entered  as  Ib/hr  are  converted  to  moles/ 
hr  and  the  Arrhenius  factors  are  precomputed  for  later  use  in  the  rate  ex- 
pressions as  are  values  for  ego  and  eQg  . However , the  primary  function 
of  the  main  program  is  to  monitor  and  control  the  execution  of  the  three- 
level  convergence  procedure  employed  to  arrive  at  the  steady  state  solution. 
The  program  is  actually  responsible  for  external  recycle  convergence  and 
calls  the  appropriate  subroutines  for  internal  recycle  or  vessel  convergence. 
The  program  also  assigns  and  updates  vessel  Identification  parameters 
so  that  it  can  correctly  sequence  through  the  reactor  system . After  overall 
convergence  is  achieved , the  main  program  prints  out  the  steady  state  values 
of  the  ideal  phase  and  actual  stream  compositions  and  flow  rates  from  each 
nitrator  and  separator  in  the  process  (including  the  off-gas  streams)  and 
then  proceeds  to  compute  and  print  out  the  values  of  a number  of  process 
performance  indicators  including  vessel  time  constants  (hours) , overall 
molar  yield,  flow  rate  of  crude  TNT  (Ib/hr) , DNT  in  the  crude  TNT  (wt  %) , 
raw  material  cost  ($/lb  TNT) , spent  acid  flow  (Ib/hr)  and  composition, 
extent  of  nitration  in  each  stage,  acid  composition  in  each  stage  (in  units 
similar  to  those  generated  by  plant  operators  during  control  checks) , ni- 
trobody  density  in  each  stage  (Ib/ft*) , and  heat  duty  in  each  nitrator. 

An  elemental  material  balance  is  also  run  to  check  the  consistency  of  the 
converged  solution.  The  complete  steady  state  printout  is  shown  in  Appen- 
dix F. 


Subroutine  NITSEP  is  employed  to  carry  out  the  Internal  recycle 
convergence  procedure.  H is  called  by  the  main  program  during  indi- 
vidual stage  convergence  and  then  again  during  total  process  conver- 
gence. The  subroutine  compares  the  current  and  previously  computed 
values  for  the  Internal  acid  recycle  leaving  a given  separator  and  does 
the  required  updating  if  the  difference  exceeds  a predetermined  value. 


The  heart  of  the  steady  state  simulation  program  is  :ubroutine 
VESSEL.  Here,  all  of  the  vessel  and  physical  separation  equations  are 
solved  to  give  the  output  of  each  nitrator  and  separator  during  execution 
of  the  nested  iteration  convergence  procedure.  The  subroutine  is  thus 
responsible  for  the  execution  and  convergence  of  the  primary  iteration 
loop,  i.e.  the  individual  vessel.  VESSEL  is  called  by  the  main  program 
and  by  subroutine  NITSEP  during  the  uncoupled  stage  convergence 
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procedure  and  then  again  by  subroutine  NITSEP  during  total  process  con- 
vergence. The  calling  program  provides  the  necessary  input  data  to  VESSEL, 
such  as  whether  it  is  to  simulate  a nitrator  or  a separator , the  values  of  all 
input  streams  and  the  currently  assumed  output  of  the  vessel.  VESSEL  will 
then  execute  Steps  4 through  19  of  the  solution  procedure  and , if  the  vessel 
has  been  flagged  as  a separator , process  the  converged  vessel  output  through 
the  physical  separation  equations . A modified  direct  substitution  procedure 
is  used , whereby  a percentage  of  the  difference  between  quantities  calcu- 
lated during  auccaasive  iterations  is  added  to  the  newly  calculated  quantity . 

The  final  aubprogram  in  the  steady  state  simulation  is  function  ENTH. 

It  serves  the  same  purpose  here  as  it  did  in  the  dynamic  aimulatlon. 

Parameter  Fitting 

Aa  has  been  illustrated,  tho  vessel  equations  ere  interspersed  with  a 
number  of  constants.  These  include  terms  associated  with  kinetic  phenom- 
ena such  as  frequency  factors,  activation  energies,  and  concentration  expon- 
ents, as  well  ss  mass  transfer  coefficients,  equilibrium  constants,  and  solu- 
bility fsotc^a.  In  most  oases,  the  appearance  of  these  constants  in  tho  vessel 
equations  is  thsoretloally  consistent  and  not  simply  a function  of  empirioal 
supposition.  An  attempt  was  made  to  use  publiahod,  exporlmsntally  determ- 
ined values  of  constants  when  such  values  were  available.  However , be- 
cause many  needed  values  were  simply  not  available,  and,  bsoause  most  of 
those  that  were  had  been  developed  under  laboratory  oonditlons  not  com- 
parable to  the  state  of  the  actual  nitration  process , it  was  necessary  to  fit. 
the  constarits  (parameters)  to  steady  state  process  operating  data  (steady 
state  snapshots) . 

The  parameter  fitting  procedure  consisted  of  (1)  assuming  an  initial 
value  for  each  of  the  27  parameters,  (2)  executing  the  steady  state  simula- 
tion with  these  pai  meter  values  and  input  variables  values  corresponding 
to  specific  data  snapshots,  (3)  comparing  computed  outputs  (molsr  com- 
ponent flows  of  all  process  streams)  with  observed  snapshot  values,  and 
(4)  making  appropriate  changes  in  parameter  values  to  give  better  agree- 
ment between  calculated  ond  observed  quantltios. 

A total  of  14  steady  state  snapshots  were  available  for  parameter 
fitting.  This  data  was  obtained  during  normal  operations  at  the  produc- 
tion facility  located  at  the  Hadford  Army  Ammunition  Plant . A typical  snap- 
shot is  shown  in  Table  I of  Appendix  J.  Of  the  14  snapshots,  seven  were 
used  for  fitting  with  the  remaining  sets  being  employed  to  chock  the  validity 
of  the  final  fit , 
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In  order  to  use  the  snapshot  data  (which  is  essenUally  stream  com-  O 

positions  in  weight  <1)  for  fitting  purposes,  it  was  first  necessary  to  convert 
from  units  of  composition  to  units  of  molar  flow  rate  for  eanh  component . 

This  was  accomplished  by  processing  the  snapshot  data  through  a material 
balance  program,  the  details  of  which  can  be  found  in  Roference  1.  With 
the  snapshot  data  in  correct  form  for  comparison  with  tht  Individual  stage 
outputs  calculated  by  the  simulation , it  was  simply  required  to  make  repeti- 
tive steady  state  runs  and  after  each  run,  to  evaluate  an  objective  function. 

This  function  was  in  effect  a measure  of  the  total  weight  deviation  from  ob- 
served conditions  in  all  nitration  stages.  Changes  were  thon  made  to  the 
parameter  set  which  drove  the  objective  function  toward  a minimum. 

Mlnlmliatlon  of  the  objective  function  was  accomplished  systematically 
by  using  a pattern  search  algorithm  which  was  atided  to  the  steady  state 
simulation  as  a subroutine  and  controlled  program  execution.  The  pattern 
search  technique  will  be  described  in  detail  in  the  next  section  of  the  report 
when  steady  state  optimization  la  discussed.  A list  of  parameters  with  their 
final  fit  values  is  included  in  Table  II  of  Appendix  J . 

APPLICATION  OF  THE  SIMULATION 

The  computer  simulation  of  the  nitration  section  has  boon  used  to 
carry  out  a variety  of  studies  on  the  nitration  process.  Most  significant 
among  these  has  been  the  use  of  the  steady  state  simulation  to  develop  op- 
timized operating  conditions  and  the  development  of  a multlvai  table  control 
procedure  for  vessel  acid  composition  using  the  dynamic  simulation  to  gen- 
erate the  process  transfer  function . 

It  is  not  Intended  that  this  report  serve  as  the  documentary  record 
of  the  exploitation  of  the  nitration  section  simulation.  Detailed  technical 
reports  on  specific  studies  using  the  mtratlon  section  simulation  have  been 
prepared  (see  Ref  15  to  17) . However,  It  is  desired  in  this  volume  to  give 
the  roeder  some  feeling  for  the  way  in  which  the  nitration  section  simula- 
tion can  be  used.  With  this  in  mind,  the  following  description  of  twomajor 
areas  of  exploitation  is  preaented. 

Steady  State  Applications 

Steady  state  optimization  was  carried  out  by  making  multiple  steady 
state  runs  under  the  control  of  the  same  pattern  search  algorithm  that  was 
used  during  parameter  fitting.  The  pattern  search  method  is  one  of  a num- 
ber of  direct  ollmlng  (Ref  13)  techniques  which  will  move  toward  an  extre- 
mum of  a multldlmenslona)  obJecMve  function  based  on  previous  information 
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about  the  response  surface.  The  particular  pattern  search  algorithm 
used  was  an  adaptation  of  a procedure  developed  by  Hooke  and  Jeevea 
(Ref  14}  which  employs  accelerated  cllmlng  and  has  rldge-followlng 
properties . 

For  each  steady  state  run  made  under  control  of  the  pattern  search 
subroutine,  an  objectlva  function,  defined  as  the  rsw  material  cost  per 
pound  of  TNT , was  evaluated . The  Independent  variables  ovsr  which  the 
search  was  made  included  feed  composition , fsed  flow  rate , nltrator  temper 
ature  and  Internal  recycle  rate.  The  feed  rate  of  toluene  to  the  process 
(and  therefore  the  production  rate)  was  not  varied  during  steady  state 
optimisation  so  that  what  was  sought  was  the  set  of  Independent  verlables 
which  gave  the  minimum  raw  material  cost  par  pound  of  a -TNT  produced 
at  a fixed  production  rate.  The  fixed  toluene  rate  employed  during  optim- 
ization was  2350  Ib/hr , which  oorresponds  to  a TNT  production  rate  of 
approximately  57  tons/day  based  on  an  87%  a -yield . The  raw  material 
cost  per  pound  of  TNT  was  calculated  from  the  aquation: 

(TxC«)  + (8xC„)  + (NxC„)  + (WxC,„) 


■ raw  material  coat  (i/lb  TNT) 

■ toluene  flow  rate  (lb/ hr) 

■ oleum  flow  rate  (Ib/hr) 

■ strong  nitric  acid  flow  rate  (Ib/hr) 

■ weak  nitric  add  flow  rate  (Ib/hr) 

■ cost  of  toluene  (I/lb  toluene) 

■ cost  of  oleum  (t/lb) 

> coat  of  strong  nitric  sold  (88.5%  HNO|)  (I/lb) 


B cost  of  weak  nitric  acid  (81.0%  HNO|)  (I/lb) 


yield  (lb  TNT/lb  toluene) 


Value!  of  coit  factor  employed  were  baaed  on  data  aupplled  by 
Radford  Army  Ammunition  Plant  and  thus  reflected  that  plant's  raw  ma- 
terial cost  an  of  July  1B72.  Specific  values  used  were:  ■>  .0293, 

C„  = .01192,  C,,,  = .01061,  C.,  = .01714. 
o W INI 

Rather  than  consider  all  Independent  variables  ilmultanaously.  a 
ssnsitlvlty  analysis  was  carried  out  which  enabled  ranking  these  vari- 
ables In  the  order  In  which  perturbations  in  them  affected  the  objective 
function.  Only  main  effect  sensitivities  wore  determined.  Although  the 
possibility  of  significant  higher  order  effects  could  not  be  discounted,  it 
was  felt  that  the  complexity  of  multiple-effect  analysis  (Including  increased 
machine  time  requirements)  was  not  juatlfled. 

The  Independent  variables  were  thus  divided  into  five  blocks  accord- 
ing to  their  effects  on  the  objective  function,  and  pattern  searches  were 
conducted  separately  over  each  block . Since  the  number  of  variables  per 
aearoh  was  reduced  by  this  procedure,  the  pattern  search  algorithm  func- 
tioned more  efficiently  and  saved  computer  time.  The  complete  Hat  of  in- 
dependent variables,  ranked  In  decreasing  order  of  their  effects  on  the  ob- 
jective function  and  broken  down  into  their  search  blocks,  is  shown  In 
Table  I of  Appendix  K , 

Because  of  physical  limitations  In  the  actual  prooesa,  a number  of 
constraints  on  the  Independent  variables  had  to  be  Included  in  the  objec- 
tive function  to  Insure  that  only  feasible  seto  of  process  operating  condi- 
tions were  considered , This  was  accompliohod  by  assigning  penalty 
values  to  the  objective  function  whenever  a constraint  value  of  one  or  more 
of  the  Independent  variables  was  under  consideration  of  the  search  rou- 
tine. This  "brute  torce"  method  worked  satisfactorily  in  keeping  th^ 
search  within  foaaible  oreos  of  the  response  surface.  The  constraints 
employed  were  based  on  Information  supplied  by  plant  operating  personnel 
at  Radford  and  are  listun  In  Table  II  of  Appendix  K, 

As  tho  optimisation  proceodod.  It  bneame  evident  that  a substantial 
reduction  In  oleum  was  being  made  by  the  search  program.  Because  of 
the  magnitude  of  the  prudlctfd  reduction  (nearly  50'kof  the  baseline  feed 
rate) , the  validity  of  the  model  In  a region  so  for  from  conditions  under 
which  parsmeturs  were  fit  came  into  question.  Also,  even  If  the  simula- 
tion's predictions  were  correct,  It  was  felt  that  plant  operating  personnel 
would  never  agree  to  such  a drastic  reduction  In  oleum  flow  , But  the  quali- 
tative result  was  clear-‘cut  back  on  oloum  flow. 
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Cnn.-cqucr.Uy , the  oleum  flow  rate  waa  treated  as  a parameter , with 
Cota  of  optimum  conditiona  being  developed  for  oleum  flow  rates  of  11,000 , 
10,000,  8,000,  8,000,  and  7,000  Ib/hr,  Of  course,  the  lowest  unit  cost  of 
TNT  is  predicted  for  the  lowest  oleum  flow  rate.  Predicted  optimum  condi- 
tions for  10,000  Ib/hr  oleum  are  shown  in  Table  K3  of  Appendix  K . For 
the  Interested  reader , a more  complete  deaorlption  of  certain  aspects  of 
the  steady  state  optimisation  studios  can  be  found  in  lieferenoe  3 . In  gen- 
eral , however , steady  state  optimization  studies  have  resulted  in  the  de- 
velopment of  operating  conditiona  which  predict  savings  on  the  order  of 
20%  in  the  raw  material  cost  per  pound  of  TNT. 

Dynamic  Applications 

The  primary  application  of  the  dynamic  simulation  has  been  its  use 
in  conjunction  with  the  design  of  a supervisory  control  strategy  for  acid 
composition  in  each  stage  of  nitration , The  need  to  control  the  composition 
of  the  acid  phase  in  each  nitration  stage  is  a result  of  the  critical  affect 
which  this  variable  has  on  the  extent  of  reaction.  The  object  of  the  control 
study  was  thus  to  develop  a prooedure  whereby  predicted  optimum  acid 
concentrations  could  be  maintained  throughout  the  process  in  view  of  un- 
controllable fluctuations  in  such  tilings  as  raw  material  acid  composition , 
internal  flows,  mixing  charaoteristlcs,  and  heat  transfer  ooeffioients , as 
wall  as  atop  changes  in  production  rate. 

It  was  determined  that  the  way  to  accomplish  this  was  to  control 
specific  indicators  of  acid  concentration  in  all  six  external  recycle  streams 
and  in  the  emulsion  stream  flowing  from  Nltrator  3A  to  Nltrator  3B  by  man- 
ipulating the  seven  primary  nitric  acid  feeds  and  the  yellow  waterfeeds  to 
the  process . A system  of  nine  manipulated  and  nine  controlled  variables 
results.  The  specific  indicator  of  acid  composition  is  optional  and  can  be 
either  total  acidity,  actual  nitric  acid,  or  nitric  to  sulfuric  acid  ratio.  Con- 
trol strategies  employing  various  combinations  of  these  controlled  variables 
have  been  tested  and  shown  to  be  feasible  and  the  moat  promising  approach 
appears  to  bo  control  of  tu  ul  acidity  and  actual  nitric  acid  in  Stages  1 and  2 
with  the  nitric  to  sulfuric  ratio  being  controlled  at  the  remaing  five  loca- 
tions. Controlled  variables  with  their  proposed  manipulated  variables  are 
listed  in  Figure  10. 

The  details  of  control  strategy  development  have  been  documented 
elsewhere  (Ref  15) , although  a number  of  sl|[*nlfiCBnt  points  should  be  noted 
here.  Because  of  the  many  recycle  streams  in  the  process,  it  was  necessary 
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to  develop  a decoupling  mechanism  in  order  to  counteract  the  disturbances 
created  in  adjacent  nitrators  when  a feed  change  is  made.  Secondly,  be- 
cause of  difficulties  in  obtaining  real  time  measurements  of  the  controlled 
variables,  it  was  also  necessary  to  develop  a computational  scheme  to  estim- 
ate the  values  of  tht  desired  controlled  variables  from  measurable  data  and 
certain  assumed  constant  values  generated  by  the  steady  state  program . 

A block  diagram  of  the  overall  scheme  for  acid  concentration  control 
is  shown  in  Figure  11,  Note  that  the  process  simulation  provides  the  neces- 
sary process  dynamics  in  the  same  way , as  an  assumed  transfer  function 
(e.g. , a first  order  lag  plus  dead  time)  would  have,  had  not  a detailed 
process  model  been  used , 

A number  of  runs  were  made  with  the  simulatod  control  scheme  in 
order  to  develop  reasonable  values  for  controller  tuning  parameters . The 
open  and  closed  loop  responses  of  the  system  to  a 10%  step  change  in  oleum 
flow  are  shown  in  Figures  12  c.nd  13.  respectively,  for  a typical  run. 

The  dynamic  simulation  was  also  used  in  conjunction  with  <n  simula- 
tion of  the  nitralor  temperature  cascade  control  loop  which  is  Included  in 
the  software  package  provided  by  the  Foxboro  Company  as  part  of  the 
Volunteer  AAP  direct  digital  computer  control  system , scheduled  to  begin 
operation  in  late  1074,  Performance  of  the  temperature  control  system 
was  tested  for  normal  continuous  operation  on  a single  nitrator  simulation 
and  also  for  start-up  uperatlons  on  the  complete  six-stage  nitration  section 
simulation.  The  details  end  significant  results  of  this  effort  can  be  found 
In  Reference  16. 


CONCLUSION 

The  modeling  and  simulation  study  of  the  continuous  TNT  process 
has  resulted  in  much  being  learned  about  the  chemical  and  physical  pheno- 
mena which  govern  the  process's  behavior.  Quantification  of  process  in- 
teractions hae  made  posslbla  tho  generation  of  optimum  operating  condi- 
tions, Improvements  in  tho  process  design,  and  has  allowed  the  develop- 
ment of  a supervisory  control  strategy . Implementation  of  the  results  of 
the  simulation  effort  is  anticipated  during  calendar  year  1974  when  nuw 
computer-controlled  TNT  lines  will  be  started  up  ut  Volunteer  Army 
Ammunition  Plant  in  Chattanooga,  Tennessee.  Only  after  verification  of 
model  predictions  on  the  actual  process  and  use  of  the  model -developod 
process  control  strategies,  will  the  full  benefit  of  the  simulation  under- 
taking be  realized.  IVhilo  simulation  studies  on  other  processes  are 
planned,  much  will  depend  on  the  successful  implementation  of  tho  TNT 
process  simulation  study's  results. 
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The  following  is  an  ordered  set  of  equations  which  make  up  the 
kinetics  and  mass  transfer  model  for  any  veeeel  in  the  nitration  section 
of  the  continuous  TNT  process.  The  order  In  which  they  are  listed  repre 
sents  the  logical  sequence  of  calculation  employed  In  SUBROUTINE  SUB 
of  the  dynamic  simulation  and  SUBROUTINE  VESSEL  of  the  steady  state 
simulation  program. 

Component  Index  Aselanments 

Component  Index  (11 


aMNT 

1 

mMNT 

2 

oDNT 

3 

mDNT 

4 

aTNT 

5 

mTNT 

6 

TNB 

7 

TNBX 

a 

HNOi 

8 

Toluene 

10 

H18O4 

11 

HNO8O4 

12 

H,0 

13 

80, 

14 
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IS 

Not  used 

16 
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Bouatloni 
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■ i.Vi 
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i*  *1 
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» E £l 
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■'p 

■ 1 - F. 

A 

(A-8) 

When  veiiel  la  a leparator,  • .758 

(A-8a) 

A. 

N-Tcir 

(A-9) 

(A-10) 

Where: 


P > 3.62 0.01023  T - 3.608  X|,  7.716  (A  ll) 

(Xi  X|)  - 8.452  (X,  + X,)  - 6.202  X, 


nA  • VPp  kma 


(A-12) 


I ' 


VF  (kma) 


^Tol  = 


i-  . _i 
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APPENDIX  C 

Dynamic  Simulation— Program  Flowchart 
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Preceding  page  blank 


I 


! 

) 


> 


Yes 


ift'ite  nitra- 
tor  ni'iiibar 
find  coupon 
fnt  flow 
leader. 


CALL  ‘^UB  1 


I 


3 


Q 


205 


m 


CAZX  8UB1 


write  mudmuffl  percent 
chnnce  in  acid  and 
organic  phase  oom- 
ponenta  for  this 
Iteration 


Compute  molar 
real  stream  com- 
ponent flows 
from  ideal  phase 
flows 


Write  stream  con^oo* 
ent  end  flow  data  for 
organic  stream  forward 
Internal  arecyole  und 
external  recycle  for 
each  stage  _ 


21 


208 


u 


m 


Read  tmcnltude.  tlnej 
hd  duration  of  pertur- 
bations to  initial  fltendyj 
state  of  the  proeeas 


•ak  Aold  Nitration? 


Datarmlna  whloh  atroh^ 

* 

mold  region  la  to  be 
anforoad  and  oomputa 
lonlo  apaoia  oompoai- 
tlona  for  aaoh  i^aaa 

- - 4 

Computa  nitro- 
nium  ion  oon- 
oantratlon  via 
waak  a61d 
aqullbrla 

Oonputo 

OANMA 


Computa  r«- 
aotlon  rata 
faotora 


Conputa  total  molar 
tarms  for  organic 
apaoiao  and  tha  oommon 
mas  tranafar  faotora 
CX  and  CZ  uaad  In  the 
Integration  aohama 


2 


Compute  naifl 
flow  rate  of  ni- 
tric plus  eq,ulv. 
Bidfurlc  In  acid 
phase  (LB^B) 


A«LBTO/tBS 


Conn)ute  adjusted  value 
of  WTA  for  use  In  lo- 
cating position  In 
the  enthalpy-onn^osi- 
tloD  data  table 


Compute  the  enthalp/ 
of  the  nltrobody  free 
sold  phase  by  inter- 
polntlon  within  the 
data  table 


Confute  the  en- 
thalpy of  organ- 
ic phase  plus 
remaining  aold 
phase  nltrobody 


Pet urn 


Main  Program 
A (I) 


molau/ht  of  component  i in  the  total  acid  phase 
leaving  the  vesael  cunently  being  calculated 


o 


AA(I.J) 

AFtI) 

Aid) 

AIR(I) 

AKMA 

AKIO 

AK12 

AK3S 

AK3W 

AK4S 

AK4W 

AK6S 

AK6S 

AK7 

AK8 


moles/hr  of  component  1 In  the  stream  phase  j 
of  the  separator  currently  being  calculated 

moles/hr  of  component  1 in  the  acid  phase  of 
ti<e  organic  atream  leaving  a separator 

moles/hr  of  component  1 in  the  total  acid  phase 
Input  to  the  vesael  currently  being  calculated 

molea/hr  of  component  1 in  the  acid  phase  of 
the  Internal  recycle  leaving  a separator 

acid  phase  mass  transfer  coefficient 

rate  constant  for  reaction  10 

ratio  of  rate  constants  for  reactions  1 and  2 

rate  constant  for  reaction  3 in  strong  acid 

rate  constant  for  reaction  3 In  weak  acid , NOT 
USED 

rata  constant  for  reaction  4 In  strong  acid 

rate  constant  for  reaction  4 In  weak  acid . NOT 
USED 

rate  constant  for  reaction  S In  strong  acid 
rate  constant  for  reaction  6 In  strong  acid 
rate  constant  for  reaction  7 
rate  constant  for  reaction  8 
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ANT(I.J) 


ARRAY  (I  .J) 

lS 


NOT  USED 


r 


percent  actual  nitric  acid  in  the  acid  phaae  of 
the  external  recycle  leaving  eeparator  j at  the 
and  of  the  Ifli  plot  Interval 

amount  of  o - nitration  taking  place  in  nitrator  j 
at  the  end  of  the  i**»  plot  interval 

running  aum  uaed  in  computation  of  ANA(I,J) 

amount  of  m-nltratlon  taking  place  In  nitrator  J 
at  the  end  of  the  plot  interval 

running  aum  uaed  In  computation  of  ANM  (I,J) 

moles /hr  component  i in  the  total  acid  phase 
leaving  nitrator  J 

moles/hr  of  component  i from  fresh  acid  feeds 
to  nitrator  j 

molea/hr  of  component  1 in  the  total  acid  phaae 
leaving  separator  J 

total  molea/hr  of  stream  phase  i leaving  the 
separator  currently  being  calculated 

total  Ib/hr  of  stream  phase  1 leaving  the 
eeparator  currently  being  calculated 

initial  value  of  moles/hr  component  J in  the 
total  acid  phase  leaving  vessel  1 

total  ou  ft/hr  of  stream  phase  1 leaving  the 
separator  currently  being  calculated 

moles/hr  of  component  1 in  the  acid  phase  ox 
the  external  recycle  leaving  a separator 

ffiolee/hr  of  component  i in  the  acid  phase 


AIM  (I) 


molei/hr  of  component  i in  the  fresh  acid  feeds 
to  nltrator  lA 


o 


■ 


AIB  (1) 

moles/hr  of  component  i 
nltrator  IB 

AlBl(l) 

moles/hr  of  component  1 
fresh  sold  feeds 

A18(I) 

moles/hr  of  component  1 
leevlng  separator  1 

A2(I) 

moles/hr  of  component  1 
nltrator  2 

A2K1) 

moles/hr  of  component  1 
acid  feeds 

A28  (I) 

moles/hr  of  component  1 
leaving  separator  2 

A3A(I) 

moles/hr  of  component  1 
nltrator  3A 

ASAKD 

moles/hr  of  component  1 
fresh  sold  feeds 

A3B  (1} 

moles/hr  of  component  1 
leaving  nltrator  3B 

A3BK1) 

moles/hr  of  component  1 
fresh  acid  feeds 

A3S  (I) 

moles/hr  of  component  1 
leaving  separator  3 

A4i:i) 

moles/hr  of  component  1 
nltrator  4 

A4S  (I) 

moles/hr  of  component  I 
leaving  separator  4 

A4I  (I) 

moles/hr  of  component  1 
fresh  acid  feeds 

in  the  add  phase  leaving 
to  nltrator  IB  from 
in  the  total  acid  phase 
in  the  acid  phase  leaving 
to  nltrator  2 from  fresh 
in  the  total  acid  phaae 
in  the  acid  phase  leaving 
to  nltrator  3A  from 
in  the  sold  phase 
to  nltrator  3B  from 
in  total  acid  phase 
in  acid  phase  leaving 
In  total  acid  phase 
to  nltrator  4 from 


u 


AS  (I) 


molei/hr  of  component  1 in  acid  phase  leaving 
nltrator  5 


I 


r 

r 


fi 


\ 

I 

! 


A8K1) 

rool«s/hr  of  component  1 to  nltrator  5 from 
froth  told  foods 

A58(l) 

moloo/hr  of  component  1 In  the  total  acid  phase 
leaving  separator  5 

A6(I) 

molet/hr  of  component  1 In  the  acid  phase  leaving 
nltrator  6 

A61(l) 

moles/hr  of  component  1 to  nltrator  6 from  freah 
acid  faede 

A68  (I) 

molee/hr  of  component  1 In  the  total  acid  phase 
from  oeparator  6 

BASIS 

NOT  USED 

CDNT 

weight  % DNT  In  crude  TNT  leaving  separator  6 

CDNT8A 

weight  % DNT  In  the  nltrobody  dissolved  In  the 
spent  sold 

CHAMXO 

NOT  USED 

CHPMX9 

NOT  USED 

CNBSA 

weight  % nltrobody  dissolved  In  the  acid  phase 
of  the  spent  acid  leaving  separator  1 

CP  (I) 

specific  heat  of  component  1 

CRUDE 

Ib/hr  crude  TNT  leaving  separator  6 (l.s. , the 
mass  flow  of  the  total  organic  stream  leaving 
separator  6) 

CWSA 

weight  % water  In  the  acid  phase  of  the  spent 
add  leaving  separator  1 

231 


DBN8  (I,J) 

Bpeclfic  gravity  of  the  organic  phase  of  any  of 
the  streams  leaving  the  Jfn  separator  at  the  end 
of  the  1th  plot  Interval 

DBNSAd.J) 

•peoiflo  gravity  of  the  acid  phase  In  any  of  the 
■treama  leaving  the  Jtti  separator  at  the  end  of 
the  1th  plot  Interval 

OFAMX0 

maximum  % change  In  an  acid  phase  component 
between  two  euoceeelve  Iteratlone  of  the  nitration 
aoctlon 

DFPMXa 

maximum  % change  In  an  organic  phase  component 
between  two  euocesalve  Iterations  of  the  nitration 
section 

DH(I) 

heat  of  reaction  for  the  1th  reaction  In  BTU/lb- 
mole  of  reactant 

DLT 

estimated  temperature  of  entrained  organic  phase 
that  would  be  returned  to  nltrator  IB  from  the 
after  separator 

DT 

Itegratlon  time  atep  (hrs) 

DTPR 

time  between  plotted  points  (mins) 

BAQ(I) 

molar  ratio  of  sold  phase  In  the  organic  atream 
to  total  organic  stream  leaving  the  1th  separator 

BA0P(I) 

molar  ratio  of  acid  to  organic  phase  in  the 
organic  stream  leaving  the  1th  separator 

BA0P1 

equivalent  to  BAOP  (1) 

EA0P2 

equivalent  to  EAOP(2) 

EA0P3 

equivalent  to  EAOPO) 

BA0P4 

equivalent  to  BAOP  (4) 

EA0P5 

equivalent  to  EAOP(5) 

equivalent  to  EAOP  (6) 


O EAQP6 

SAP 

BPB 

B9A(I) 

BOAPd) 

BQAPl 

BaAP2 

BOAPa 

B0AP4 

B9AP6 

B9AP6 

B9L 

B9F 

BWN 

BX 

BYW 
BIO 
62 1 


value  of  EAOP  (I)  for  the  separator  currently 
being  oaloulated 

enthalpy  of  mixing  for  fume  recovery  acid 

molar  ratio  of  organic  phaaa  in  the  external 
recycle  to  total  external  recycle  leaving  the 
aeparator  i 

molar  ratio  of  organic  to  acid  phase  In  either  of 
the  acid  recycle  itreame  leaving  the  ith 
separator 

equivalent  to  BOAP  (1) 

equivalent  to  BOAP  (2) 

equivalent  to  BOAP  (3) 

equivalent  to  BOAP  (4) 

equivalent  to  BOAP  (5) 

equivalent  to  BOAP  (6) 

enthalpy  of  mixing  for  40%  oleum 

value  of  BOAP  (I)  for  tliS  aeparator  currently 

being  calculated 

enthalpy  of  mixing  for  weak  nitric  ucld 

Impcaed  upper  bound  on  the  Arrhenius  co- 
efficient in  the  reaction  rate  expreesiona  for 
reactions  3 and  4 (MNT  to  DNT) 

enthalpy  of  mixing  for  yellow  water 

activation  energy  for  reaction  10 

ratio  of  activation  energies  for  reaotinns  2 and  1 
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E3S 


activation  energy  for  reaction  3 in  a atrong 
acid  medium 


B3W 

B4S 

B4W 

BS8 

B68 

B7 

B8 

B9 

PA  (I) 
FACIO(I) 

PACIDl 

FACID3 

FACID3 

7AC1D4 

FACID6 


NOT  USED 

activation  energy  for  reaction  4 in  a etrong 
acid  medium 

NOT  USED 

activation  energy  for  reaction  6 in  a atrong 
acid  medium 

activation  energy  for  reaction  6 in  a atrong 
acid  medium 

activation  energy  for  reaction  7 
activation  energy  for  reaction  8 
NOT  USED 


rptlo  of  acid  phaae  leaving  1^^  eeparator  that  li 
recycled  to  total  acid  phase  leaving  the  ith 
aaparator  (fraction  of  acid  phase  recycled) 

ratio  of  acid  phase  leaving  the  1^^  separator  in 
the  organic  stream  forward  to  total  acid  phase 
leaving  the  1^^  eeparator  (fraction  of  acid  phase 
sent  forward) 

equivalent  to  FACIO(l) 

equivalent  to  FACID(2) 

equivalent  to  FACIDO) 

equivalent  to  FACID  (4) 

equivalent  to  FACID  (8) 


FACID6 


equivalent  to  /ACID  (6) 
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I 


PAFR 


fraction  of  total  acid  phaae  leaving  aeparator  1 
that  appears  In  the  external  recycle  from  that 
aeparator 


PAR 

PAl 

PA2 

PAS 

PA4 

PAS 

FAS 

FD 

PlN(l.J) 

FMAX 

PMX 

POd) 

F0RO(I) 

Fanoi 


fraction  of  total  acid  phase  leaving  the  aeparator 
that  la  recycled  for  the  aeparator  currently 
being  oalouletad 


equivalent  to  PA  (1) 
equivalent  to  PA  (2) 
equivalent  to  PA  (3) 
equl'^alent  to  FA  (4) 
equivalent  to  PA  (8) 
equivalent  to  PA  (6) 


extent  of  reaction  for  nltroeylaulfurio 
daoompoaltlon 


Ib/hr  of  feed  t to  nitration  vessel  J 


maximum  allowable  voluroatrlo  flow  of  internal 
acid  recycle 


equivalent  to  FMAX 


ratio  of  organic  phase  leaving  the  ith  aeparator 
with  the  organic  atream  forward  to  total  organic 
phaae  leaving  the  lib  separator  (fraction  organic 
phaae  forward) 


i I 


ratio  of  organic  phase  leaving  the  1th  separator 
with  the  acid  recycle  atreama  to  total  organic 
phaae  leaving  the  itb  separator 


equivalent  to  P0RO  (1) 
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equivalent  to  PQR0(2) 
equivalent  to  F9R0(3) 
equivalent  to  P9RO(4) 
equivalent  to  F9RO(B) 
equivalent  to  PQRQ  (0) 
equivalent  to  P9(l) 
equivalent  to  PQ(2) 
equivalent  to  F90) 
equivalent  to  P9  (4) 
equivalent  to  P9  (S) 
equivalent  to  P9  (6) 

fraction  of  organic  phaee  leaving  a aeparator 
that  appeara  in  the  organic  atroam  forward  for 
the  aeparator  currently  being  calculated 

ratio  of  external  recycle  etream  to  total  recycle 
leaving  the  i^h  aoparator 

ratio  of  internal  recycle  etream  to  total  recycle 
leaving  the  aeparator  currently  being  calculated 

equivalent  to  PRC  (1) 

equivalent  to  PRC  (2) 

equivalent  to  PRC  (3) 

equivalent  to  PRC  (4) 

equivalent  to  PRC  (B) 


vj 


,1 


PRAC6 
PRC  (I) 

FRl 

PR2 

PR3 

FRl 

PRS 

PR6 

PV(I) 

PVKM 

FVKM8 

PVXMIA 

PVKMIB 

FVKM2 

PVKMSA 

FVKMSB 

FVKM4 

FVXMS 

FVKM0 

OMl) 


equivalent  to  FRC  (6)  f 

ratio  of  Internal  raoyole  atraam  to  total  raoyola 
atrean  laaving  tha  iUi  aeparator 

equivalent  to  FR  (1) 

equivalent  to  FR  (2) 

equivalent  to  FR  (S) 

equivalent  to  FR  (4) 

equivalent  to  FR  (B) 

equivalent  to  FR  (6) 

1 

volume  of  the  if^  nltrator  (ou  ft) 

volume  of  the  veaiel  currently  being  oaloulated  | 

volume  of  a aeparator  (all  aeparatora  have 
the  aame  volume) 

volume  of  nltrator  lA;  equivalent  to  FV(1) 

volume  of  nltrator  IB;  equivalent  to  FV  (2) 

volume  of  nltrator  2:  equivalent  to  FV  (3) 
volume  of  nltrator  3A;  equivalent  to  FV  (4) 

volume  of  nltrator  3B;  equivalent  to  FV(S) 

volume  of  nltrator  4;  equivalent  to  FV  (6) 

volume  of  nitretcr  8;  equivalent  to  FV  (7)  | 

volume  of  nltrator  6;  equivalent  to  FV  (8)  | 

ratio  of  Arrheniue  ooefflolenta  In  the  kinetlo 

rate  expreealone  lor  reastion  1 and  2 In  tha 

li*'  nltrator  i 
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010(1) 

03(1) 

038(1) 

048(1) 

048(1) 

05(1) 

05(1) 

07(1) 

05(1) 

09(1) 

1 

1C 

IGALL 

IC9UNT 

IN 

IPNCH 

IPRNT 


value  of  the  Arrhenlue  coefficient  In  the  kin  itlc 
rate  expreaalon  for  reaotion  10  In  the  1th  nitrator 

NOT  U8BD 

value  of  the  Arrhenius  ooefflolent  In  the  Itlnetlo 
rate  expression  for  reaction  3 In  the  ith  nitrator 
under  strong  sold  oondltlone 

NOT  U8BD 

same  as  038(1)  but  tor  reaction  5 

value  of  Arrhenius  ooefflolent  In  the  klnetlo 
rate  expression  for  reaction  5 In  the  1th  nitrator 

same  as  05(1)  but  for  reaction  B 

same  as  05  (1)  bur  for  reaction  7 

aama  as  05  (1)  bur  tor  raaotion  8 

NOT  U8BD 

00  LOOP  Index 

oounter  tor  plotting  Intervals 

flag  for  first  call  to  SUBROUTINE  INPT 

equivalent  to  1C 

DO  LOOP  Index  for  certain  loops  pertaining 
to  nltratore 

I/O  device  numbu  for  punched  cards 
print  control  flag: 

0 summary  printout  only 

1 debug  printout  of  Internal  calculations 
during  each  Integration  step 


i 

1 


33B 


number  of  stages  being  simulated 

DQ  LOOP  index  for  certain  loops  pertaining 
to  ■eparatora 

atage  number 

veaael  type  (lag: 

lTyPB*0  for  nitrator 
ITYFB«1  for  aaparitor 

nitration  vaaaal  number 


NO(l) 


NOTCH 


general  veaael  number 

DO  LOOP  counter  and  generally  used  integer 
operator 

counter  for  integration  time  itepa 

DO  LOOP  counter  and  generally  used  integer 
operator 

molecular  weight  of  component  1 

Integer  number  of  integration  paaaea  through  the 
prooeaa  (or  iterations)  required  for  a particular 
run . Alao  the  integer  number  of  Integration  time 
stops  required  for  the  run 

integer  number  of  iterations  between  summary 
printouts 

notch  setting  on  internal  rsoyola  gate  valve  for 
separator  1 

notch  setting  for  the  separator  currently  being 
oalculated 

Integer  number  of  iterations  between  computations 
of  process  performance  indicators  and  other  data 
for  plotting 


N1 

equivalent  to  N9(l) 

N2 

equivalent  to  N9(2) 

Na 

•qul valent  to  N9(3) 

N4 

equivalent  to  NQ(4) 

N6 

equivalent  to  NQ(B) 

N$ 

equivalent  to  NQ(6) 

P(I) 

molee/hr  of  component  1 In  the  total  organic 
phaee  leaving  the  vaaeel  currently  being 
oaloulated 

PA 

Ib/hr  of  total  aold  phoae  leaving  a aeparator 
for  which  aold  phaae  apeolflo  gravity  la  being 
oalaulated 

PF(I) 

molee/hr  oomponent  1 in  the  organic  phaae  of 
the  organic  atream  leaving  a aeparator 

PI  (I) 

molea/hr  oomponent  1 In  the  total  organic  phaae 
Input  to  the  veeael  currently  being  calculated 

Pin  (I) 

molea/hr  of  component  1 In  the  organic  phaae  of 
the  Internal  reoyole  leaving  a aeparator 

PNT  (I.J) 

nolea/hi  of  component  1 In  the  total  organic 
phaae  leaving  nitrator  J 

PP 

Ib/hr  of  total  orgenio  phaae  leaving  a aeparator 
for  which  orgenio  phete  apeolflo  gravity  le 
being  oaloulated 

PSd.J) 

molee/hr  of  oomponent  1 In  the  total  organic 
phaae  leaving  the  jih  aeparator 

PV  (I.J) 

Initial  value  of  molea/hr  component  J In  the 
total  organic  phof  'navlng  veaael  1 

o 


PXR  (I) 

moles/hr  of  component  1 in  the  organic  phase  of 
the  external  recycle  leaving  a separator 

PIA(I) 

moles/hr  of  component  i in  the  organic  phase 
leaving  nitrator  lA 

PIAI(I) 

molee/hr  of  component  i in  the  organic  feed 
stream  to  nitrator  lA  (moles/hr  toluene  to  lA) 

FIB  (T) 

moles/hr  of  component  i in  the  organic  phase 
leaving  nitrator  IB 

FIBI(I) 

moles/ hr  of  component  i in  the  organic  feed 
stream  to  nitrator  IB  (moles/hr  toluene  to  IB) 

PIS  (I) 

moles/hr  component  i in  the  total  organic  phase 
leaving  separator  1 

P2  (I) 

moles/hr  component  i in  the  organic  phase 
leaving  nitrator  2 

P2S  (I) 

moles/hr  component  i in  the  total  organic  phase 
leaving  separator  2 

F3A(I) 

moles/hr  component  1 in  the  organic  phase 
leaving  nitrator  3A 

P3B  (I) 

moles/hr  component  i in  the  organic  phase 
leaving  nitrator  3D 

F3S(I) 

moles/hr  component  i in  the  total  organic 
phase  leaving  separator  3 

P4(I) 

moles/hr  component  i in  the  organic  phase 
leaving  nitrator  4 

P4S  (I) 

moles/ hr  component  i in  the  total  organic 
phase  leaving  separator  4 

P5(I) 

moles/hr  component  1 In  the  organic  phase 
leaving  nitrator  5 
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SSirntm 


PSS  (1) 

moles/hr  component  1 In  the  total  organic 
phase  leaving  separator  5 

P6(I) 

molei/hr  component  1 in  the  organic  phase 
leaving  nitrator  6 

P6S(I) 

molea/hr  component  i in  the  total  organic 
phase  leaving  separator  6 

Q 

accumulation  term  used  in  calculation  of 
nitrator  heat  loads 

QA 

volumetric  flow  of  total  aold  phase  leaving  a 
separator  for  which  acid  phase  specific  gravity 
is  being  calculated 

QN(I.J) 

heat  load  in  nitrator  J calculated  at  the  and  of 
the  l^R  plotting  interval 

QP 

volument  flow  of  total  organic  phase  leaving  a 
separator  for  which  organic  phase  specific 
gravity  is  being  calculated 

QY(I) 

fraction  of  nitric  acid  converted  to  nitronlum 
ions  in  nitration  vessel  1 

R(I) 

NOT  USED 

FATES  (I,J) 

rate  of  reaction  j in  nitration  vessel  i 

(I) 

molar  density  of  component  i 

AHQl 

equivalent  to  RHO  (1) 

RH(?10 

equivalent  to  PHQ(IO) 

RHQtll 

equivalent  to  RHO(H) 

RH012 

equivalent  to  RH0(12) 

RH013 

equivalent  to  RHQX13) 

RH914 

equivalent  to  RH0(14) 
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-yO  ,V- V,‘ /vj'.k*;.  - • V 


RH02 


RH03 
RHQ4 
RHOS 
RHQ6 
RH07 
RHQ8 
RM09 
RN8  (I,J) 


RN88  (I.J) 


aqul valent  to  RHO  (2) 
equivalent  to  RHQ*  (3) 
equivalent  to  RHO  (4) 
equivalent  to  RH06 
equivalent  to  RHO  (6) 
equivalent  to  RHO  (7) 
equivalent  to  RHO  (8) 
equivalent  to  RHO  (9) 

nitric  to  aulfurlc  acid  ratio  In  nitration  veaael 
j calculated  at  the  end  of  the  ith  plotting  interval 
(for  acid  phase  only) 

nitric  to  auifurlo  acid  ratio  in  separator  j 
calculated  at  the  end  of  the  llh  plotting  Interval 
(for  acid  phase  only) 

coefficient  on  the  molecular  oxygen  term  In 
reaction  9 

Ideal  gas  law  constant 

adjustable  constant.  NOT  USED 

equivalent  to  organic  phase  mass  transfer 
coefficient  (KMAP) 

equivalent  to  AKEQ4  in  SUBROUTINE  SUB 
NOT  USED 

equivalent  to  AK9.  NOT  USED 


equivalent  to  E9 , NOT  USED 


SH 


equivalent  to  AKEQl  in  SUBROUTINE  SUB 
NOT  USED 


equivalent  to  AKBQ3  in  SUBROUTINE  SUB 

eooumulator  term  uaad  in  oaloulatlon  of  apent 
acid  maae  flow  rate 

apent  acid  maaa  flow  rata  at  the  and  of  the  1^^ 
plotting  interval  (flow  rate  of  aoid  phaae  of  the 
external  recycle  atream  leaving  aeparator  1) 

maaa  flow  rate  of  diaaolvad  organic  in  the  apent 
aoid 

accumulator  term  uaad  in  oaloulatlon  of  the  maaa 
flow  rate  of  crude  TNT  leaving  aaparator  8 

temperature  of  nitration  veaael  ourrantly  being 
calculated 

maaa  flow  rate  of  DNT  in  the  total  organic  atream 
leaving  aaparator  6 

cumulative  total  molaa/hr  oxidea  of  carbon  (gaa) 
leaving  all  veaaela  in  the  nitration  aaotlon  during 
a given  iteration 

temperature  of  the  i^  nitration  veaael 
NOT  USED 

almulated  elapaed  time  aince  beginning  of  run 

total  Ib/hr  of  aoid  phaae  leaving  a aeparator  on 
a nitrobody  free  baala.  Uaed  in  computation  of 
aoid  concentration  indioatora 


total  molaa/hr  of  atream  1 leaving  a aeparator 
currently  being  calculated 
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TNOX 

TR 

TSBPd) 

TTCOX(I) 

TTMB 

TTNOX(I) 

TVL(I) 

TW  (1) 

TIA 

TIB 

T2 

T3A 

T3B 

T4 

T6 

T6 

VKM  (I) 
VKMS 


total  iH'^laa/hr  oxidai  of  nitrogen  (gae)  leaving 
all  vess/rla  In  the  nitration  section  during  a given 
Iteration 


raferanoe  temperature 

oomputad  temperature  of  the  separator  currently 
being  calculated 

same  as  T€OX  except  that  value  is  particular 
value  at  the  end  of  the  1^^  plotting  Interval 

total  time  cl  process  operation  to  be  simulated 

same  as  TNOX  except  that  value  Is  particular 
value  at  the  end  of  the  1th  plotting  Interval 

volumetric  flow  (cu  ft/hr)  of  the  1*^  stream 
leaving  the  separator  currently  being  calculated 

mass  flow  (Ib/hr)  of  the  it^  stream  leaving  the 
separator  ourrently  being  calculated 


temperature  of  nltrator  lA  (equivalent  to  TBMP  (1) ) 
temperature  of  nltrator  IB  (equivalent  to  TBMP  (2) ) 
temperature  of  nltrator  2 (equivalent  to  TBMP  (3) ) 
temperature  of  nltrator  3A  (equivalent  to  TBMP  (4)) 
temperature  of  nltrator  3B  (equivalent  tc  TEMP  (6) ) 
temperature  of  nltrator  4 (equivalent  to  TEMP  (6) ) 
temperature  cf  nltrator  5 iequivalent  to  TEMP  (7) ) 
temperature  of  nltrator  8 (equivalent  to  TEMP  (8) ) 
<^olume  of  the  i**'  nltrator 


volume  oi  a separator  (same  for  all  separators) 
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W 


i 

I 


I 


X(l) 

elapsed  simulated  time  at  the  end  of  the  1^^ 
plotting  Interval  (nomenclature  Is  specific  to 
plotting  routine  for  plotting  one  dimensional 
variables  as  a function  of  time) 

XX(l.J) 

elapsed  simulated  time  at  the  end  of  the  1^^ 
plotting  interval  (this  form  of  elapsed  time  is 
used  for  plotting  variables  of  up  tu  j diroenelon) 

XYY(l) 

moles /hr  nitric  aoid  in  the  sold  phase  converted 
to  nitronlum  ions  in  nitration  vessel  1 

Subroutine  INPT 

The  definttione  of  the  following  varleblee  are  Identical  to  thoee  in  the  main 
program:  AF(I) , AIRd) . AKMA,  AXMATL,  AKIO.  AX12.  AX38.  AX3W, 

AX4a,  AK4W,  AK58.  AK68.  AK7.  AX8.  AK8,  ARRAY(I.J),  AXR(I) , A18  (I) , 

A2S(1),  AS8(I).  A48(l).  AS8(I),  A68(I).  CP(1) , DH(I),  DT.  BAOPl.  BA0P2.  BAOP 
BAOPS.  BAaP4,  BAQP6.  BA9P6,  BQAPl.  BaAP2.  SOAPS.  B9AP4,  SOAPS. 

BOAPS.  BX.  BIO.  B21.  B3S.  B3W,  £48.  B4W.  B58.  £68.  E7,  BB.  E9.  FACIDI. 
FAC1D2.  FAC1D3.  FAG1D4.  FAC1D5.  FAC1D6.  FAl . FA2,  FAS,  FA4.  PAS, 

FAB,  PIN(l.J),  FMX,  FOROl,  F0R02.  F0R03.  FOR04.  FQROS,  FOHQ6, 

FOl.  P02,  P03.  F04,  FOB.  FOB,  FRACl,  PRAC2.  FRAGS,  FRAC4,  FRAGS, 

FRAGB,  FRl,  PR2,  PR3,  PR4,  FR8,  FR6.  FVKMIA,  FVKM1B,FVKM,?. 

FVKM3A,  PVKM3B,  FVKM4,  FVKMB,  FVKMB.  01(1),  010(1),  03(1).  038(1), 

04(1),  048(1).  OB(l),  06(1).  07(1).  08(1).  08(1),  ICALL,  IPRNT,  JAGK, 

MW(1),  NM.  NN.  NP.  Nl,  N2,  N3.  N4.  NS,  N8,  PP(1),  PlR(l),  PXR(l), 

PlAI  (1) . PIBI  (I) , P18  (1) . P2S  (I) , P38  (1) , P48  (1) . PS8  (1) , PBS  (1) , 

RATES (1,J),  RHOl,  RHOlO,  RHOll.  RH012,  RH013,  RH014,  RHQ2.  RH03, 

RH04,  RHOS,  RHOB,  RH07,  HHOB,  RH08,  RP,  RR,  SA,  SB,  SC,  SD,  SE, 

SO,  8H,  81.  8J,  TGOX,  TEMP  (I).  TNOX,  TR.  TSBP(I),  TTME,  TIA,  TIB, 

T2,  T3,  T4,  TS.  T3, 

Remaining  variable!  are  defined  ai  follows: 

COMP(l.J)  weight  component  1 in  feedstock  typo  ] 

DDL  fractional  part  of  a change  In  feed  rate  or  nltrator 

temperature  which  la  applied  to  the  proceas 
during  a given  iteration 

DL  (I)  magnitude  of  the  ith  change  In  an  indepenent 

variable 
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ICH(J) 

Integer  code  number  for  the  independent 
variable  type  to  which  the  change  ia  being 

made 

IKT 

flag  to  indicate  whether  or  not  temperature 
chimgee  have  been  made  to  the  initial  oonditione 

IPRNT 

aifflulated  time  between  summary  printouts 

IVS(J) 

integer  code  number  for  the  vessel  to  which  the 
jth  change  in  one  of  the  independent  variable 
types  is  mads 

NCHNO 

number  of  given  change  which  is  to  be  made 
during  the  run 

NOTCH  (I) 

notch  setting  on  internal  recycle  gate  value  for 
the  ith  separator 

NTMB  (1) 

number  of  iterations  between  the  beginning  of 
a run  the  initiation  of  the  ith  change  in  an 
independent  variable 

NTVL  (1) 

number  of  iterations  over  which  the  magnitude 
of  the  ith  ohanga  is  to  be  equally  distributed 
(used  for  ramping  an  independent  variable) 

TME 

simulated  time  at  which  a given  change  is  to 
be  initiated 

TVL 

simulated  time  interval  over  which  a given 
change  is  to  be  made 

ZA 

temperature  coefficient  in  the  exponent  of  the 
Arrhenius  rate  expression 

Subroutine  SUB 

The  definitione  of  the  following  verliblee  ere  Identical  to  those  in  the 
main  program:  A (I),  AI(l).  AKMA.  AKMATL.  AKIO.  AX12,  AX38.  AK3W, 
AK48.  AX4W.  AK68.  AK68,  AK7,  AX3.  CP  (I).  DH(I),  DT,  BX.  BIO,  B21. 
BSW,  B88.  B4W,  B48,  BBS,  B7,  BB.  PAil,  FD.  ?MX,  FP9,  BBS,  FRAC,  FVKM, 
OKI),  010(1).  03(1),  038(1).  04(1),  048(1).  05(1),  03(1),  07(1).  OB  (I), 
08(1),  ICALL.  IPRNT.  I8T0B,  ITYPB,  IV.  JACK.  NM,  NN,  NOTCH.  NP, 

P(I),  PI(1).  QY(I),  RATBSd.J),  RHOl.  RHOH.  RH012,  RH013,  RH014, 
RH02,  RH03.  RH94,  RHOB,  RH96.  RH97,  RHOl.  RHOB,  RR.  P.P.  8A,  SB 
SC,  8D.  SB.  SO,  8H.  SI.  SJ.  T.  TCOX,  TI.  TNOX,  TR,  T8BP(I),  TTMB, 
XYY(l) 

The  definitions  of  the  following  variables  are  identical  to  other  variables 
dadnad  for  tha  main  program; 


Variable  in  SUB 


Variable  in  main  oroaram 


All  thru  A117 
A1  thru  A17 

BAOP 

B9AP 

Pll  thru  PI17 
PI  thru  P17 


AI  (1)  thru  AK17) 
A(l)  thru  A (17) 
BAP 
BOP 

Pl(l)  thru  PI (17) 
P(l)  thru  P (17) 


Remaining  variables  are  defined  aa  follows: 


AINB 


moles /hr  nltrobody  in  the  aggragate  sold  phase 
input  to  the  vessel  currently  being  oaloulated 


general  rate  factor  for  all  nitration  reactions . 
Also  if  Bennett's  strong  >aoid  rate  expression 
Is  used , the  uoaffiolent  uf  the  beaulfate  ion 
concentration 

coefficient  of  the  sulfuric  acid  concentratiii^n 
in  Bennett' K rate  expression 

coefficient  of  the  pyro  sulfate  ion  conoentration 
in  Bennett's  rate  expression 


SM  i"' 


r 


AKEQ 

equilibrium  constant  for  the  dissociation  of 
sulfuric  acid  in  an  aqueous  medium 

AKBQOV 

equilibrium  constant  for  the  dissociation  of  a 
mixture  of  aulfurlo  and  pyroaulfurio  acids  to 
ionio  species.  NOT  USBD 

AKBQl 

equilibrium  constant  for  the  dissociation  of 
nltrlo  sold  in  aqueous  sulfuric  acid 

AKBQ2 

equivalent  to  AXBQ 

AKBQ3 

equilibrium  constant  for  the  dissooiation  of 
nitric  sold  in  anhydrous  sulfuric  acid 

AXBQ4 

equilibrium  constant  for  the  recombination  of 
nltronium  and  blaulfata  ions  in  anhydrous  acid 
mlxtruea 

AXMAP 

organic  phase  mass  transfer  coefficient 

AKl 

computed  value  of  AKBQl  obtained  from  interval 
halving  procedure 

AK3 

computed  value  of  AKBQ3  obtained  from  interval 
halving  proosdura 

AK4 

computed  volua  of  AKBQ4  obtained  from  interval 
halving  procedure 

AL 

molas/hr  water  in  the  total  acid  phase  of  the 
vessel  current  being  calculated;  can  have 
negative  value  which  then  represents  molss/hr  SO| 

AMNB 

moles/hr  nitrobody  transferred  from  organic 
to  acid  phase  by  bulk  mass  transfer;  negative 
value  means  transfer  from  sold  to  organic  phase 

AMI 

tnoles/hr  uMNT  transferred  by  bulk  mass  transfer 

AM2 

same  as  AMI  but  for  noMNT 
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L 


AMS 

AM4 

AMS 

AMS 

■M7 

AMS 

ANB 

A804 

AT 

AUQA 

B 

BC 

C 

CB 

CF 


1 I 

i 

i 

I 

I 


lame  ai  AMI  but  for  aDNT 

■ama  aa  AMt  but  for  mDNT 

aana  aa  AMI  but  for  oTMT 

aana  aa  AMI  but  for  mTMT 

auna  aa  AMI  but  for  TNR 

aana  aa  AMI  but  for  TNBX 

total  molaa/hr  of  nltrubody  dlaaolvad  In  tha  add 
pbaaa  leaving  tha  vaaaal  ourrantly  being  oaloulatad 

ootnbinad  molaa/hr  of  80|  and  H|SO«  in  tha  add 
phaaa  leaving  tha  vaaaal  ourrantly  being  oaloulatad 

total  molaa/hr  of  add  phaaa  leaving  tha  vaaaal 
ourrantly  being  oaloulatad 

oonoantratlon  (molaa/ft* ) of  aulfurlo  add  in  tha 
aoid  phase  of  tha  vaaaal  ourrantly  being  oaloulatad 

initial  oonoantratlon  (tnolaa/ft*)  of  traa  aulfurlo 
aoid  in  tha  aoid  phaaa  of  tha  vaaaal  currently 
being  oaloulatad 

3 -t-  (3  X C) 

initial  oonoantratlon  (molaa/ft* ) of  pyroaulfurio 
add  in  tha  aoid  phaaa  of  tha  vaaaal  ourrantly 
being  oaloulatad 

(.5  X B)  ^ (1.8  X C) 


L.) 


I 


j 

)' 


C-P 

aqulllbrium  oonoantratlon  of  nitric  add  in  tha 
add  phaae  of  tha  vaaaal  ourrantly  balng  oaloulatad 


380 


CHN03 


o 


CHS04 


GH8207 

CH20 

CH2804 

CH30 

0N02 

CX 

02 

CS 

CPS 


equilibrium  oonoentratlon  nf  blsulfate  lone  in  the 
told  phase  of  the  vessel  ourrently  being  calculated 


equilibtltun  oonoentratlon  of  pyrosulfate  Ions  In 
the  told  phase  of  the  vessel  ourrently  being 
oaloulatod 

same  so  CHNOS  but  for  water 

equilibrium  oonoentratlon  of  sulfuric  sold  in  the 
sold  phase  of  the  vassal  currently  being  oaloulated 

equilibrium  cmosnfration  of  hydronlum 
ions  in  the  sold  phase  of  the  vessel  our*> 
rently  being  osloulatsd 

equlllbriuia  ocnoentration  of  nitronium  ions  in 
the  sold  phsQS  of  the  vessel  currently  being 
oaloulated 

oonaolid'^tton  *erm  fur  (2A  ^ YP) 

oonsolidation  term  for  (2B  - YA  4 i ♦ cX) 

moles/hr  nitrosylfulfurio  sold  produced  by  oxi- 
dation resotlons  In  the  vessel  ourrsnUy  being 
oaloulated 

moles/hr  nitrosylsulfurio  sold  which  ramsin 
undaoumposed  in  the  vessel  ourrently  being 
csloulated 


|1 


i 


i' 


1 

1 


D 

initial  sold  phase  concentration  of  combined 
sulfuric  sold  plus  pyrosulfurio  acid  m the 
vessel  ourrently  being  usloulsted 

DAL 

rate  of  ohsnge  of  AL  with  raspeot  tu  time  in  the 
vessel  ourrently  being  osloulstad 

DA804 

rate  of  ohsnge  of  AS04  with  rospeot  to  time  in  the 
vossol  ourrently  being  colculot^d 

DAI 

rate  of  change  of  A1  with  respeot  to  time  in  the 
vessel  ourrently  being  oaloulated 

2B1 

rftt«  of  ohango  of  A12  with  raopeot  to  time  In  the 
vaaaal  ourrantly  being  oaloulatod 

rata  of  change  of  A2  with  raapaot  to  time  in 
vaaaal  ourrantly  being  oaloulated 

rata  of  change  of  AS  with  raapaot  to  time  in  the 
vaaaal  ourrantly  being  oaloulated 

aama  aa  DAS  bin  for  A4 

aame  ar  DAS  but  for  AS 

aama  aa  DAS  bu'c  '^r  AS 

aama  aa  DAS  but  for  A7 

aame  aa  DAS  but  for  AS 

aama  aa  DAS  but  fot  \9 

aoid  concentration  v m in  tha  weak  acid  nitration 
rata  axpraaalon . Nut  U8BD 

rata  of  change  of  PI  with  raapaot  to  time  in  the 
vaaaal  ourrantly  being  oaloulated 

aama  aa  DPI  but  for  PIO 

aame  aa  DPI  but  for  PZ 

aamo  aa  DPI  but  for  P3 

aamo  aa  DPI  but  for  P4 

aama  aa  DPI  but  for  PS 

aame  aa  DPI  but  for  PO 

aama  aa  DPI  but  for  P7 


aama  aa  DPI  but  for  PB 


O DTA 
DTP 
D1 

DIO 

D2 

D3 

D4 

D5 

DO 

D7 

DO 

ETA 

ETAA 

ETAP 

ETATL 

F 

FA 


acid  phase  adjusted  integration  time  step 

organic  phase  adjusted  Integration  time  step 

moles/hr  of  o -MNT  that  diffuse  Into  the  sold 
phase  in  the  vessel  currently  being  calculated 

moles/hr  of  toluene  which  diffuse  from  the  bulk 
organic  phase  to  the  aoid-organlo  interfaos 

asms  as  D1  but  for 

same  as  D1  but  for 

same  as  D1  but  for 

same  as  D1  but  for 

same  as  D1  but  for 

same  as  D1  but  for 

same  as  D1  but  for 

overall  mass  transfer  coefficient  in  the  vessel 
currently  being  calculated 

acid  phase  mass  transfer  coefficient  in  the  vessel 
currently  being  calculated 

organic  phase  mass  transtar  coefficient  in  the 
vessel  currently  being  calculated 

mass  tiansfer  ocefficient  for  toluene 

initial  acid  phase  concentration  of  nitric  acid  in 
the  veseel  currently  being  calculated 

volumetric  fraction  of  aci ' phase  in  the  vessel 
currently  being  calculated 
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4 


m-MNT 

a -DNT 

m-DNT 

0 -TNT 

m-TNT 

TNB 

TNBX 


I'. 


PF(I) 


POV 

PP 


FlOA 


FiOM 

P3 


F4 

F8 

F6 

F7A 


F7M 

FBA 


F6M 

F9A 

F9M 


cross  sectionsl  ares  fraction  corresponding  to  the 
ith  notch  sotting  on  tho  gets  vslvs  in  ths  internal 
raoyols  lino  of  tho  sopsrstor  currently  being 
oslculstod 


Qonatont  sot  equal  to  one 


volumotrlo  fraction  of  organic  phase  in  the  vessel 
currently  being  osloulatsd 


raaotion  rate  factor  for  reaction  10  where  a -TNT 
is  the  reactant  for  the  vessel  currently  being 
calculated 


same  as  FlOA  but  for  m-TNT 


reaction  rate  factor  for  reaction  3 in  the  vessel 
currently  being  celoulated 


same  as  P3  but  for  reaction  4 
same  as  F3  but  for  reaction  5 
same  as  P3  but  for  reaction  6 


reaction  rate  factor  for  reaction  7 in  the  vessel 
currently  being  calculated  where  a -DNT  is 
the  reeotent 


same  as  F7A  but  for  m-DNT 


reaction  rate  factor  for  reaction  8 in  the  veseel 
currently  being  calculated  where  n -DNT  Is  the 
reactant 


same  as  F6A  but  for  m-DNT 

Zero 

Zero 


I 
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GAMMA 


acid  concentration  term  used  in  all  nitration 
rate  expressions 


i 


KLM  flag  to  Indicate  which  nitrstion  region  exists 

in  the  vessel  being  calculated 

FNB  total  flow  (moles/ hr]  of  organic  phase  nitrobody  jj 

components  leaving  the  vessel  currently  being 
calculated 

PP  exponent  of  nitrobody-in-acid-  solubility 

correlfadon 

Q fraction  nitric  acid  converted  to  nitronium  ions 

for  the  vesael  currently  being  calculated 

QA  volumetric  flow  of  the  acid  phase  leaving  the 

vessel  currently  being  calculated 

volumetric  flow  of  the  organic  phase  leaving 
the  vessel  currently  being  calculated 

total  volumetric  flow  of  combined  internal  and 
external  recycle  leaving  the  separator  currently 
being  calculated 

rate  of  the  ith  reaction  in  the  vessel  currently 
being  calculated 

rate  of  reaction  1 in  the  vessel  currently  being 
calculated;  equivalent  to  RATE  (1)  (moles/hr) 

same  as  R1  but  for  reaction  10;  equivalent  to 
RATE  (10);  also  RlOA  + RIOM 


QT 

RATE  (I) 

HI 

RIO 


RlOA 

rate  of  reaction  10  in  the  vessel  currently  being 
calculated  where  a -TNT  Is  the  reactant 

RIOM 

same  as  RlOA  but  for  m-TNT 

R2 

same  as  R1  but  for  reaction  2;  equivalent  to 
RATE  (2) 

I 
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I 


H3 

flame  as  R1  but  for  reaction  3;  equivalent  to 
RATE  (3) 

R4 

same  as  R1  but  f^r  reaction  4;  equivalent  to 
RATE  (4) 

R5 

lame  ae  R1  but  for  reaction  5;  equivalent  to 
RATE  (5) 

R6 

lame  ae  R1  but  for  reaction  6;  equivalent  to 
RATE  (6) 

R7 

same  as  R1  but  for  the  combined  reaction  7; 
equivalent  to  RATE  (7) ; also  R7A  R7M 

R7A 

rate  of  reaction  7 in  the  vessel  currently  being 
calculated  where  a -DNT  is  the  reactant 

R7M 

same  as  R7A  but  for  m-DNT 

R8 

same  ae  R1  but  for  the  combined  reaction  8; 
equivalent  to  RATE  (8);  also  R8A  ■*-  H8M 

RBA 

same  as  R7A  but  for  reaction  6 

R8M 

same  ae  R8A  but  for  m-DNT 

H9 

same  as  R1  but  for  reaction  9;  equivalent  to 
RATE  (9) 

R9A 

Zero 

R9G 

Zero 

R9M 

Zero 

TKEL 

temperature  of  the  vessel  currently  being 
calculated  C*K) 

T1 

A1  + PI 

T2 

A2  + P2 
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f 


T3 

A3  + P3 

T4 

A4  + P4 

T6 

AS  + P5 

T6 

A6  -t-PB 

T7 

A7  +P7 

TB 

AB  -t-PB 

W 

Initial  add  phase  concentration  of  water  In  the 
vessel  currently  being  calculated  (moler>/!t‘) 

X 

consolidation  term  used  in  the  Integration  sequence 
of  organic  species  in  the  acid  phase  of  the  vessel 
currently  being  calculated 

XEQA 

equilibrium  solubility  of  nltrobody  in  the  acid 
phase  of  the  vessel  currently  being  calculated 
(moles/hr) 

XL 

lower  limit  of  an  ionic  species  concentration  which 
la  applied  during  internal  halving  procedure 

XU 

upper  limit  of  an  Ionic  species  concentration 
which  Is  applied  during  Internal  halving  procedure 

XUP 


XI 

X13 

X14 


upper  limit  of  HS|Oif  concentration  (reaultlng 
from  dissociation  according  to  equation  19) 
which  is  applied  during  Internal  halving  procedure 


mole  fraction  a -MNT  In  the  acid  phase  of  the 
vessel  currently  being  calcul  \ted 


mole  fraction  water  In  acid  phase  of  the  vessel 
currently  being  calculated 


mole  fraction  sulfur  trioxide  In  the  acid  phase 
of  the  vessel  currently  being  calculated 


I 


t 

i 

i 


i 

) 

I 


I 


'i 

i . 
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mole  fraction  m-MNT  in  the  acid  phase  of  the  I 

vessel  currently  being  calculated 

mole  fraction  o -TNT  In  the  acid  phase  of  the  I 

vessel  currently  being  oaloulated  | 

mole  fraotlon  m-TNT  in  the  sold  phase  of  the  j 

vessel  currently  being  calculated  | 

i 

mole  fraotlon  nitric  acid  in  the  acid  phase  of  ' 

the  vessel  currently  being  oaloulated 

oonsolidation  term  used  in  the  integration 

sequence  for  organic  speolss  in  the  sold  phase  i 

of  the  vessel  currently  being  calculated;  also  j 

used  for  initial  conorntxation  (molos/ft* ) of  > 

free  sulfuric  in  the  aoid  phase  of  the  vessel 
currently  being  oaloulated 

ratio  of  AMNB  to  ANB  | 

ratio  of  AMNB  to  PNB 

initial  acid  phase  oonoentratlon  of  nitric  aoid  I 

which  is  used  under  strong  aoid  nitrating  | 

conditions  for  which  water  is  present;  also  < 

used  as  a oonsolidation  term  in  the  integration  i 

sequence  for  organic  species  in  the  acid  phase  I 

enthalpy  of  the  contents  of  a separator 

(BTU/hr  - •C)  , 

rate  of  heat  generated  in  a separator  due  to  I 

reaction  (BTU/Hr) 

AB-A17;  undissociated  acid  phase  nitric  acid 
remaining  at  equilibrium  in  the  vessel  cur- 
rently being  calculated  (moles/ hr) 
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Subroutine  SUB  2 


The  definitions  of  the  following  variables  are  identical  to  those  In  the  main 
program:  AKMA.  AKMATL.  AKIO.  AK12.  AK3S,  AK3W,  AK48.  AK4W, 
AKS8.  AX68,  AK7,  AK6.  CP  (I).  DH(I).  BIO.  E21,  B38.  E3W.  B48.  E4W. 
E5S,  BBS,  E7.  E8.  FD,  DMX,  QA,  QP.  RHO(I).  BP.  RR.  SA,  SB.  SC,  SD. 
SB.  SO,  SH.  SI,  SJ,  TR 


Remaining  variables  are  defined  as  follows: 


Add) 


EAQP 


EQAP 


moles/hr  oumponent  1 in  the  total  acid  phase 
leaving  the  separator  for  which  recycle  para- 
meters are  currently  being  calculated 

total  tnoles/hr  of  acid  phase  leaving  the  separator 
for  which  recycle  parameters  are  currently 
being  calculated 

molar  ratio  of  acid  to  organic  phase  in  the  organic 
stream  leaving  the  separator  for  which  recycle 
parameters  are  currently  being  calculated 

molar  ratio  of  organic  to  acid  phase  in  either  of 
the  recycle  streama  leaving  the  separator  for  which 
recycle  parameters  are  currently  being  calculated 

fracticn  cf  acid  phase  that  is  recycled  in  the  separ- 
ator for  which  recycle  parameters  sre  currently 
being  calculated 


FARl 


fraction  of  acid  phase  sent  forward  from  the  separ- 
ator for  which  recycle  parameters  are  currently 
being  calculated 


FF(I) 


identical  lo  definitions  given  in  subroutine  SUB 

fraction  of  organic  phase  sent  forward  from  the 
separator  for  which  recycle  parameters  are 
currently  being  calculated 


FPQl  fraction  of  organic  phase  which  is  recycled  in 

the  separator  for  which  recycle  parameters  are 
currently  being  calculated 

FRAC  ratio  of  Internal  recycle  stream  to  total  recycle 

streams  leaving  the  separator  for  which  recycle 
parametere  are  currently  being  calculated 

NOTCH  notch  setting  on  the  internal  recycle  gate  valve 

for  the  separator  for  which  recycle  parameters 
are  being  calculated 

PO  molee/hr  component  1 the  total  organic  phase 

leaving  the  separator  for  which  recycle  para- 
meters are  currently  being  calculated 

PT  moles /hr  total  organic  phase  leaving  the  separator 

for  which  recycle  parameters  are  currently  being 
calculated 

QT  volumetric  flow  (ft* /hr)  of  combined  recycle 

streams  leaving  the  separator  for  which  recycle 
parametere  are  currently  being  calculated 


Subroutine  8UB1 


The  definitions  of  the  following  variables  are  identical  to  those  in  the  main 
program;  AI(J) , AKMA,  AXMATL,  AKTO.  AK12.  AK38,  AK3W,  AK4S, 

AK4W,  AK5S,  AK6S,  AK7,  AK8.  AV(l.J).  CHAMX9,  CP  (I) , DFAMXO.  DH  (I)  , 
DT,  EX,  ElO,  E21,  E38,  E3W,  E4S.  E4W.  E5S.  EOS.  E7,  EB,  FD,  FMX,  01(1), 
OlO  (I) . G3  (I) . 038  (I) , 04  (I) , 048  (I) , 05  (I) , 06  (I) , 07  (I) , 08  (I) , 09  (I) , 
I^ALL,  IPRNT,  IV8L,  JACK,  NM,  NN,  NP,  PI  (I),  PV(I.J),  QY(I),  RATES  (I,  J), 
RhOl.  RH(910,  RHOill,  RH(212,  RHQ13,  RH914,  RH92,  RH93,  RH94,  RH(?5, 
RH96,  RH(?7,  RH9B,  RHOS,  RP.  RR.  SA,  SB.  SC,  SD,  SE,  SO,  SH,  SI,SJ, 
TCOX,  TNOX,  TR,  T8EP(I) , TTME,  XYY  (I) 

Remaining  variables  are  defined  as  follows; 

A(?  (I)  molee/hr  component  1 in  the  acid  phase  of  the 

vessel  for  which  component  % changes  are 
currently  being  calculated 
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DFA(I] 


OFAMX 


DFP(I) 


DFPMX 


POd) 


Ri-nio 


acid  phase  component  i % change  during  NN 
Iterations  In  the  vessel  currently  being  calculated 

maximum  % change  In  an  acid  phase  component 
In  the  vessel  currently  being  calculated 

same  as  DFA  (I)  but  for  organic  phase 

same  as  DPAMX  but  for  organic  phase 

moles/hr  component  1 In  the  organic  phase  of 
the  vessel  for  which  component  % changes  are 
currently  being  calculated 

NOT  USED 


Function  ENTH 


The  definitions  of  the  following  variables  are  Identical  to  those  In  the  main 
program:  AKMA.  AKMATL,  AKIO,  AK3S.  AK3W.  AK4S,  AK4W,  AKSS. 
AK6S,  AK7,  AKB.  CPd).  DNd).  ElO,  E12.  E3S,  E3W,  E4S.  E4W,  ESS, 
E6S,  E7,  E8,  FD.  PMAX.  RHOl.  RHOlO,  RH(?11.  RHQ12,  RH013,  RH014, 
RH92.  RH03.  RHq4,  RHOS,  RH06,  RHQ7,  RHOB.  RH9B.  RR.  SA.  SB,  SC. 
SD.  SE.  SO,  SH,  SI.  SJ,  TR 

Remaining  <^ariablss  arc  defined  as  follows: 


AX  d) 


moles/hr  component  1 in  the  acid  phase  of  the 
stream  for  which  the  enthalpy  Is  currently 
being  calculated 

moles /hr  component  1 In  the  total  acid  phase 
leaving  a vessel  for  which  the  enthalpy  of  a 
particular  leaving  stream  Is  currently  being 
calculated 


equivalent  to  CP  (11) 
equivalent  to  CP  (12) 
equivalent  to  CP  (13) 
equivalent  to  CP  (14) 


^ 9- 


DLT 

ENTH 

BNTHO 

ENTH2 

F 

LBS 

LBSE 

P(l) 

PX(I) 


T 


temperature  of  the  vessel  (relative  to  20**C) 
for  which  the  enthalpy  of  a particular  exiting 
stream  ie  currently  being  calculated 

enthalpy  of  the  stream  currently  being  calculated 
(kcal/hr) 

enthalpy  of  mixing  of  the  acid  components  in  the 
acid  phase  of  the  stream  currently  being  calculated 
based  on  the  zero  nitric  curve  of  the  McKinley  and 
Brown  data  (koal/lb) 

enthalpy  of  mixing  of  the  acid  components  in  the 
acid  phasa  of  the  stream  currently  being  calculated 
based  on  the  20%  nitric  curve  of  the  McKinley  and 
Brown  data  (kcal/lb) 

adjusted  weight  fraction  of  nitric  acid  plus  sulfuric 
acid  in  the  acid  phasa  of  the  stream  currently 
being  calculated  on  a nitrobody  free  basis 

Ib/hr  of  nitric  acid,  sulfuric  acid,  water  (or 
sulfur  trloxide)  in  the  acid  phase  of  the  stream 
currently  being  calculated 

Ib/hr  of  nitric  acid  and  equivalent  sulfuric  acid 
in  the  acid  phase  of  the  stream  currently  being 
calculated 

molea/hr  component  i in  the  organic  phase  of  the 
stream  for  which  the  enthalpy  is  currently  being 
calculated 

moles/hr  component  i in  the  total  organic  phase 
leaving  a vessel  for  which  the  enthalpy  of  a 
particular  leaving  stream  is  currently  being 
calculated 

moles/hr  combined  organic  and  acid  phase  nitro- 
body component  in  the  stream  for  which  enthalpy 
is  currently  being  calculated 


J 
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temperature  of  the  vesdel  for  which  the  enthalpy 
of  an  exiting  etream  Is  being  calculated  ('*0 

LB8E/LBS:  nitrobody  free  weight  ftactlon  of 
equivalent  eulfurlo  plua  nitric  acid  in  the  acid 
phaae  of  the  atream  for  which  the  enthalpy  la 
currently  being  calculated 

weight  fraction  of  nitric  acid  in  the  equivalent 
eulfurlo  plua  nitrlu  acid  of  the  acid  phaae  of  the 
etream  for  which  the  enthalpy  la  currently  being 
calculated 

the  1th  (iata  point  on  the  zero  % nitric  curve  from 
the  enthalpy  concentration  data  of  McKinley  and 
Brown 

FACID , FA  or  1 dependent  on  the  origination 
point  of  the  atream  for  whloh  the  enthalpy  la 
currently  being  calculated 

FR , FRAC  or  1 dependent  on  the  origination  point 
of  the  atream  for  which  the  enthalpy  la  currently 
being  otlculated 

F0 . F0RO  or  1 dependent  on  the  origination  point 
of  the  atream  for  which  the  enthalpy  la  currently 
being  calculated 

2 FR . FRAC  or  1 dependent  on  the  origination  point 

of  the  etream  for  which  the  enthalpy  la  currently 
being  calculated 

the  ifh  data  point  on  the  20%  nitric  curve  from  the 
enthalpy  concentration  data  of  McKinley  and  Brown 
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Data  Raad 

Format  j 

1 

MT12.  AK3V.  AKl^V.  AIBS.  AXl5;«  AK? 

, a»8,  aklo  Sno.o  ! 

2 

\K3fl,  AK*tT,  R33,  B‘»S 

BPIO.O 

3 

WSl,  K3V*  B*»W.  35«i  TS^n,  157,  bB, 

BIO  SPIO.O 

U 

HA»  8B.  3C.  SO.  n.  SO.  SII.  SI 

8P10.0 

5 

IS,  TO,  TS,  AIWA,  AlWATt,  TO,  RP 

8P10.0 

1 

7 

OH  (1),  >1.  10 

8F10.0 

B*  ; 

SH?I  fl^,  10,  !•' 

8P10.0 

10.  11 

CP  (I),  >1.  1** 

8P10.0 

r* 

■ 

Of 

r% 

AlA  (I),  I-l,  17 

5R16.7/5M(?.7Ah1'!.7,  3F10.0 

■A 

- 

V 

Ik  - l< 

PIA  (1),  T«r.  17 

5’iil'.:/3Fl'5.7Ani?.7,  3F10.0 

■ 

17  • 1; 

AIB  (I),  >1,  17 

5F1''^.7/5'^1''*7Ap1^*7,  3F10.0 

1 

,1 

20  • 

PIU  (T),  M . 17 

5Bi6.T/57i<>*7/''r:i^.7,  3FIO.0 

f 1 

23  • 25 

MS  (I),  T«l.  17 

5El6.7/5Kl6.7An^,7,  3F10.0 

' 1 

2^  • 2fl 

P18  (1),  •^M,  17 

5»!16.7/5E1<^*7A71^.7,  3F10.0 

1 

29  ■ 31 

A2  (I),  IM,  17 

5ni'i.7/3I!l^.7Al!1^.7,  3F10.0 

32  • 3*' 

P2  I-l,  17 

5E1B.7/5R16.7/''M'S.7,  3F10.0 

35  • 37 

A2n  (i[\,  >1,  17 

5B16.7/5E1<S.7A"51^-7,  3F10.0 

3B  - ItO 

P8T  (1),  I-l,  17 

5E16.7/5E16.7Aei6.7,  3F10.0 

i^K 

kl  - 1*3 

•3.‘.  (1),  I-l,  17 

5F16.7/5E18.7AeiB.7,  3F10.0 

lilt  . li(^ 

P3*'  '1).  I"l,  17 

5E16.7/5E1^.7AE1«.7,  2F10.0 

1*7  - h) 

A3B  (T),  I-l,  17 

5B16.7/5R1#5.7Aei^.7,  3F10.0  ) 

50  ■ 52 

P3B  (I),  1-1,  17 

5R1^.7/5El'^-7Arar'*7,  3F10.0  | * 

i 1 

53  - 55 

A3S  (T),  I-l,  17 

5E1'*’‘.7/5E1^.7A''1^*'^.  3F10.O  j ; 
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f 


Format 


Card 

- 58 
59  - 

fiZ  . 

- 70 

Tl  . 7? 

* * I# 

7'i  - 
77  - 79 

- ft2 
‘!3  - 85 

- «A 

A,’'  -.  'll 
o:  - ‘.'U 

'^5 

o< 

07 

90 

70 

lor 
101 
n2  - 
10?  - nB 
119 

320-(l20+fTCraW) 


Data  Raad 
P3S  (I),  1-1,  17 
AU  (I),  I-l,  17 
Fit  (I),  1-1,  17 
aUs  (I),  I-l,  17 
P»*B  (I),  1-1,  17 
A5  (I),  1-1,  17 
‘ P5  (1),  I-l,  17 
A59  (I),  I-l,  17 
P5S  (I),  1-1,  17 
,\6  (I),  I-l,  17 
(1).  1-1,  17 


5Elf:.7/5El6. 7/41516. 7.  3P10.0 
5l5l6.7/rEl6.7/4El6.7,  3P10.0 
5116. 7/5E16. 7/4816. 7,  3F10.0 
5B16.7/5E16. 7/4816. 7,  3F10.0 
5816.7/5816.7/4816.7,  3F10.0 
5816.7/5816.7/4816.7,  3F10.0 
5816.7/5816.7/4816.7,  3F10.0 
5816,7/5816.7/4816.7,  3F10.0 
5816.7/5816.7/4816.7,  3F10.0 
5816.7/5816. 7/4816. 7,  3F10.0 
5816. 7/5816. 7/4816.7.  3FIO.C 


a6o  (1),  >1,  17 

P6fl  fl),  I-l,  17 
Nl,  wa,  N3,  N4.  N5.  116 
(l),  1-1,6 
(I),  I-l,  6 
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APPENDIX  C 

StMdy  Stat«  Simulation— Program  PloMchart 


Praciilliii  Ml*  Mink 


Main  Program 


Read  values  of 
frequency  factors  and 
activation  energies 
for  use  In  kinetic 
rate  expressions 


^X^Read  values  of 
'pther  constants  used 
during  execution  of 
the  vessel,  equations 


ead  phase  oomposl 
Ions  for  each  vessel 
which  represent  triltla 
guess  at  a converged 
steady  state  solution 


Write  values  of 
kinetic  and  mass 
transfer  constants, 
and  component  physi- 
cal properties 


5 


O 


341 


10 


Sat  vaaaal 

parnmetars 

for 

nltrator  6 


Set  up  input 
tnolar  component 
flow  veotora 
for  nltrator  6 


i 


NC5<NC 


/Separator  5 ^ 
iarnal  raayola  atraa 
Vw  convargad?  ^ 


Has 

^,*<rfiparat,or  ax-\s, 
♦.«rnal  raoyola  atreM 
V.  •'onvarged? 


^^Haa 
^.^^parator  3 
mernal  raoycle 
sat ream  converged^ 


Sat  veaaal 
Parana tors 
for 

nltrator  4 


Stora  currant 
valuta  of  axtar- 
nal  raoycla  mo- 
lar oomponani 


348 


Store  current  values 
of  external  recycle 
component  molar 
flcrfs  leaving 
separator  ? 


‘’^et  up  Input 
component  molar 
flow  vectors 
for  nltrator  ? 

I 

CALL  NT'T'S'^F 

Set  up  Input 
component  molar 
flows  vectors 
for  nltrutor  1A 


CALL  VESSEL 


o»t  vessel 
parameters 
For 

nltratoi'  1B 


het  up  input 
component  molar  i 
flow  vectors 
for  nltra.tor 


CALL  NITSI'T 


u'hroutlnp  V’?SS"L 


Inltlallzfi  con 
vergence  para- 
meters! 


Sot  up  molar  component 
working  storage  ar- 
ray, V(l),  for  acl'i 
an'i  organic  phases  of 
the  vessel  currently 
being  computet 


Transfer  molar  compo- 
nent valurp  from 
working  s+onge  to 
the  variable  names 
uset  turlng  vessel 
execution 


irO”N'f-o 


A.-p 

oRn 

vessel 

temper- 

a 1 1 1 

ft.', 

verse  1 

vol  unir , 

the 

Arrhenius  rate 

f-ic 

tors 

f-'r  the  V 

sol 

curr'-n*  1 v 

be  1 lie 

Lzii 

evil -I 

t >>  i 

35 


6 


Separator?^ 


Set  'Uffusion 
factors  ' 


Compute  diffu- 
sion factors  re* 
qulrad  for  vas- 
sal aquation  so* 
lutlon  prooadur 


Compute  all  ooefflol- 
ents  required  for  sl- 
multaneoua  eolutlcn  of 
the  vessel  equations 
and  proceed  to  solve 
these  equations  via 
the  procedure  outlin- 
ed In  the  text 


III-TII  + 1 


,<^he  vessel  equa 
tlons  converged? 


yes 


Ai-i  separator  1 Intav- 
nal  recycle  molar  com- 
ponent vectors  anl  af- 
ter separator  recycl  e 
nltrohoiy  molar  compo- 
nent vector  to  nitra- 
tor  Input  molar 
component  flow 
vectors 


Compute 

wta-lmse/lbs 


Compute  ai  "usted  value 
of  VTA  for  use  In 
lon.'itln«  position  In 
tho  enthalpy-composi- 
tion data  table 


Compute  the  enthalpy 
of  the  nitroboiy  free 
acid  phase  by  Inter- 
polation within  the 
data  table 


Compute  the  enthalpy 
of  the  organic  phase 
plus  the  remain  In/? 
acl'l  phase  nitroboiy 
to  make  up  total 
ntream  enthalpy 


Peturn 


3 


s 

i 

J 
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APPENDIX  H 


StMdy  Statt  Simulation— Program  Nomenclature 


L' 

k 
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r 


Main  Program 
AAS  (I) 


AF(I.J) 


AFLBS(I.J) 


molea/hr  of  component  i in  the  acid  phase 
flowing  from  the  after  separator 

extent  of  nitration  (either  a or  m)  across  the 
stage  currently  being  calculated  (moles/hr 
NOi'*'  added  to  the  ring) 

molaa/hr  of  component  i in  the  entrained  acid 
phase  of  the  organic  stream  flowing  from 
separator  j 

same  as  AF  but  in  units  of  lbs/ hr 


AFl(l) 


AF2(1) 
AF3  (I) 


moles/hr  of  component  i in  the  entrained  acid 
phase  of  the  organic  stream  flowing  from 
separator  1 

same  as  API  but  for  separator  2 
same  as  AFl  but  for  separator  3 


AF4  (I) 
AFS(I) 
AF6(I) 


same  as  API  but  for  separator  4 

same  as  AFl  but  for  separator  5 

same  as  AFl  but  for  separator  6 

total  acid  phase  moles/ hr  of  component  1 
flowing  into  a nitrator  which  is  followed  by 
another  nitrator 


All  (I) 


AIR  (I,J) 


total  acid  phase  moles/hr  (less  acid  phase  of 
internal  recycle)  of  component  1 flowing  into 
a nitrator  which  is  followed  by  a separator 

moles/hr  of  component  1 in  the  acid  phase  of 
the  internal  recycle  stream  flowing  from 
separator  j 


374 


AIR  LBS(I.J) 
AIR  1(1) 

AIR2(I) 

AIRS  (I) 

AIH4  (I) 
AIRS(l) 

A1R6  (I) 

AKA 

AKB 
AKC 
AKEQ 
AKEQOV 

AKMA 
AKMATL 
AKIO 
AK12 

I 

! p 

I 

t 

I ' 

I 

I 

L 


same  as  AIR  but  in  units  of  Ibs/hr 

molea/hr  of  component  i in  the  acid  phaae  of 
the  internal  recycle  stream  flowing  from 
separator  1 

aame  as  AIRl  but  for  separator  2 

same  as  AIRl  but  for  separator  3 

same  as  AIRl  but  for  separator  4 

same  as  AIRl  but  for  separator  5 

same  as  AIRl  but  for  separator  6 

general  rate  factor  for  all  nitration  reactions . 
Also,  if  Bennett's  strong  acid  rate  expression 
is  used,  the  coefflolent  of  the  bisulfate  ion 
concentration 

coefficient  of  the  sulfuric  acid  concentration 
in  Bennett's  rate  expression 

uoefficient  of  the  pyrosulfuric  ion  concentration 
in  Bennett's  rate  expression 

equilibrium  constant  for  the  sulfuric  acid-water 
dissociation  reaction 

equilibrium  constant  for  the  dissociation  of  a 
mixture  of  sulfuric  and  pyrosulfuric  acids  to 
ionic  species.  NOT  USED 

acid  phase  mass  transfer  coefficient 

mass  transfer  coefficient  for  toluene 

rate  constant  for  reaction  10 

ratio  of  rate  constants  for  reactions  1 Q 2 


3 


AK3S 

AK3W 

AK4S 

AK4W 

AK5S 

AK6S 

AK7 

AK8 

AK9 

AN  (1) 

ANS  (I) 

ANSI  A 

\ 

ANS3A 
ARRAY  (I.J) 
' AS  (I) 


rate  constant  for  reaction  3 In  strong  acid 

rate  constant  for  reaction  3 in  weak  acid 
NOT  USED 

rate  constant  for  reaction  4 in  strong  acid 

rate  constant  for  reaction  4 in  weak  acid 
NOT  USED 

rate  constant  for  reaction  5 
rate  constant  for  reaction  6 
rate  constant  for  reaction  7 
rate  constant  for  reaction  8 
NOT  USED 

percent  actual  nitric  acid  in  tho  acid  phase 
of  the  external  recycle  stream  flowing  from 
the  ith  separator 

mole  fraction  nitrosylsulfuric  acid  in  the  acid 
phase  of  the  external  recycle  stream  Rowing 
from  the  ilh  separator 

mole  fraction  nitrosylsulfuric  acid  in  the  acid 
phase  flowing  from  nltrator  lA 

mole  fraction  nitrosylsulfuric  acid  in  the  acid 
phase  flowing  from  nltrator  3A 

molea/hr  of  component  1 in  the  raw  material 
feedstocks  flowing  to  nltrator  j 

percent  sulfuric  acid  and  SO|  (as  equivalent 
sulfuric)  in  the  acid  phase  of  the  external 
recycle  stream  flowing  from  the  ith  aeparator 
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o 


K 

1 


I 


AT 

AT^M  (1,J) 
AV(I.J) 

AVI(I.J) 

AXR(I,J) 

AXRLBS(I.J) 


total  moles/hr  of  acid  phase  flowing  from  the 
veaael  for  which  the  holdup  is  being  computed 

number  of  i-type  atoma  in  component  J 

molea/hr  of  component  i in  the  total  acid  phaae 
flowing  from  veaael  j 

firat  gueaa  at  converged  value  of  molea/hr  of 
component  i in  the  acid  phaae  flowing  from 
veaael j 

molea/hr  of  component  i in  the  acid  phaae  of 
the  external  recycle  stream  flowing  from 
separator  j 

same  as  AXR  but  in  units  of  Ibs/hr 


AXRl  (I) 

AXR2(I) 

AXR20(I) 

AXR3(I) 

AXR39(I) 

AXR4(1) 

AXR4Q(I) 

AXRS(I) 

AXRBO(I) 

AXR6  (I) 


molea/hr  of  component  1 in  the  acid  phase  of 
the  external  recycle  stream  flowing  from 
separator  1 

same  as  AXRl  but  for  separator  2 

previous  value  of  AXR2(I)  to  which  current 
value  is  compared  during  the  external  recycle 
convergence  procedure 

same  as  AXRl  but  for  separator  3 
same  as  AXR2(f  bit  fpr  AXR3 
same  as  AXRl  but  for  separator  4 
same  as  AXR20  but  for  AXR4 
same  as  AXRl  but  for  separator  S 
same  as  AXR20  but  for  AXR 5 
same  aa  AXRl  but  for  separator  6 
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AXR6Q  (I) 


same  as  AXRIQ  but  for  AXR6 


AlA(l) 

AlAI  (I) 

AIB  (I) 

AlBl(I) 

AIS(I) 


A2(I) 

A21(l) 

A2S(I) 

A3A  (I) 

A3A1  (I) 

A3B(I) 

A3BI(1] 

ASS  (1) 

A4(l) 

A4I(I) 

I A4SU) 

\ 

i A5(I) 

ASl(I) 
A5S(I) 


J 


moles/ hr  of  component  i in  the  total  add  phase 
flowing  from  nltrator  lA 

moles/hr  of  component  1 flowing  Into  nltrator  lA 
from  the  raw  material  feadetooks 

■ame  as  AlA  but  for  nltrator  IB 

aame  aa  AlAI  but  for  nltrator  IB 

moles /hr  of  component  1 In  the  total  acid  phaae 
flowing  ilrom  separator  1 

same  as  AlA  but  for  nltrator  2 

aame  as  AlAI  but  for  nltrator  2 

lame  ai  AIS  but  for  separator  2 

same  as  AlA  but  for  nltrator  3A 

aame  as  AlAI  but  for  nltrator  3A 

same  as  AlA  but  for  nltrator  3B 

same  as  AlAI  but  for  nltrator  3B 

same  as  AIS  but  for  separator  3 

same  AlA  but  for  nltrator  4 

same  as  AlAI  but  for  nltrator  4 

same  as  AIS  but  for  separator  4 

aame  as  AlA  but  for  nltr.itor  6 

same  aa  AlAI  bur  for  nltrator  5 

aame  aa  AIS  but  for  aeparator  R 


4 
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U A8(I) 

A6I  (1) 

A6S  (I) 

BTU  (I) 
BTUQ(I) 

CDNT 

COMP  (I.J) 

CP  (I) 
CRUDE 

CSL 

CSN 

CTF 

CTL 

CWN 

DEL 

DH(I) 

OIF 


same  as  AlA  but  for  nitrator  6 

same  as  AlAI  but  for  nitrator  6 

aame  as  AIS  but  for  separator  6 

total  heat  load  in  nitrator  1 , BTU/hr 

heat  evolved  In  nitrator  1 which  la  due  to 
chemical  reaction  only . BTU/hr 

weight  % DNT  in  the  nltrobody  stream  flowing 
from  separator  6 on  and  acid  free  basis 

weight  % component  1 In  raw  material  feedstock 

type  J 

heat  capacity  of  component  1 . BTU /lb-mole/% 

total  mass  flow  rate  (Ibs/hr)  of  the  organic 
stream  from  separator  6 

coat  of  oleum,  $/lb 

cost  of  strong  nitric  add.  I/lb 

cost  of  raw  materials  per  unit  of  TNT  produced , 
$/lb  a-TNT 

cost  of  toluene,  $/lb 

cost  of  weak  nitric  acid.  $/lb 

convergence  tolerance  for  components  in  acid 
recycle  streams 

heat  of  reaction  for  reaction  1,  BTU/lb-mole 
of  reactant 

storage  variable  for  the  difference  between 
atoms  of  a given  type  in  and  atoms  of  that 
type  out  during  material  balance  calculations 
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w 


r 

\ 


DLT 

temperator  of  separator  1 referenced  to  20**C 

HAO  (I) 

molar  ratio  of  add  phase  In  the  organic  atream 
to  total  organic  stream  flowing  from  the  1^^ 
separator 

EAQP(I) 

molar  ratio  of  acid  phase  to  organic  phase  in  the 
organic  atream  flowing  from  the  llh  separator 

EFR 

enthalpy  of  mixing  for  fume  recovery  acid 

EQA  (1} 

molar  ratio  of  organic  pha^a  in  the  external 
recycle  atream  to  the  total  external  recycle 
stream  flowing  from  the  separator 

EOAP  (I) 

molar  ratio  of  organic  phase  to  acid  phase  In  the 
external  recycle  stream  flowing  from  the  1^^ 
separator 

E0L 

enthalpy  of  mixing  for  oleum 

EWN 

enthalpy  of  mixing  for  weak  nitric  acid 

EYW 

enthalpy  of  mixing  for  yellow  water 

BIO 

activation  energy  for  reaction  10 

E21 

ratio  of  activation  energy  for  reaction  2 to 
that  of  reaction  1 

E3S 

activation  energy  for  reaction  3 

E3W 

NOT  USED 

E4S 

activation  energy  for  reaction  4 

E4W 

NOT  USED 

ESS 

activation  energy  for  reaction  5 

E6S 

acllvatlon  energy  for  reaction  6 

3BJ 


r? 

E8 

E9 

FA 

FACD(l) 

FACT2 

FACTS 

FACT4 

FACT5 

FACTS 

FD 

FDN 

FIN  (I,J) 
FMAX 

FO 

FOV 
GA  (I.J) 


activation  energy  for  reaction  7 
activation  energy  for  reaction  6 
NOT  USED 

volume  fraction  of  acid  phase  in  the  vessel  for 
which  holdup  is  being  calculated 

NOT  USED 

relaxation  factor  used  in  convergence  of  the 
external  recycle  from  separator  2 

same  as  FACT2  but  for  separator  3 

same  as  FACT 2 but  for  separator  4 

same  as  FACT2  but  for  separator  5 

same  as  FACT2  but  for  separator  6 

extent  of  reaction  for  the  decomposition  of 
nitiosylsulfuric  acid 

moles/hr  of  ring  associated  NO|  groups  flowing 
into  a nitration  stage  with  the  organic  stream 
from  the  previous  stage 

Ibs/hr  of  feedstock  i to  nitration  vessel  ] 

maximum  possible  volumetric  flow  rate  of 
material  through  any  internal  recycle  pipe 

volume  fraction  of  organic  phase  in  the  vessel 
for  which  holdup  is  being  calculated 

NOT  USED 

number  of  1-type  atoms  in  gaseous  product  J 


CUS  (I.J) 

moles/hr  of  gaseous  product  J evolved  in  vessel  i 

Gl(I) 

ratio  of  the  Arrhenius  coefficients  in  the  kinetic 
rate  expressions  for  reactions  1 and  2 in  the 
i^  nltrator 

010(1) 

Arrhenius  coefficient  for  reaction  10  In  nltrator  1 

015(1) 

moles/hr  of  COi.i  evolved  from  vessel  1 

016(1) 

moles/hr  of  NO^  evolved  from  vessel  1 

018(1) 

moles/hr  of  TNM  evolved  from  vessel  i 

019(1) 

moles/hr  of  CO|  evolved  from  vessel  i 

03(1) 

NOT  USED 

OSS  (I) 

Arrhenius  coefficient  for  reaction  3 in  nltrator  i 

04(1) 

NOT  USED 

04S  (I) 

same  as  OSS  but  for  reaction  4 

08(1) 

same  as  038  bur  for  reaction  6 

0G(I) 

same  as  OSS  but  for  reaction  8 

07(1) 

same  as  OSS  but  for  reaction  7 

08(1) 

same  as  03S  but  for  reaction  6 

09(1) 

Zero 

HDR  (I.J) 

array  used  to  store  alphanumeric  headings 
for  chemical  elements  in  material  balance 
printout 

I 

general  index  integer 

lAFT 

after  separator  switch . If  equal  to  1 then  after 
separator  is  considered  in  the  simulation  as  n 
dilution  vebbsl  to  which  1000  Ib/hr  of  yellow 
water  is  added 

382 


1 


o 


ICHECK 


ch(  ck  calculation  awltch . If  equal  to  1 then 
check  calculations  are  performed  on  the  vessel 
equations  in  SUBROUTINE  VESSEL 


1 


IE 

IFLG  (1) 

II 

III 
U 

ILIST 

IMBAL 

IN 

INIT(I) 

IPI 


entrainment  correlation  switch . If  equal  to  0 
then  entrainment  correlation  is  used  instead 
of  assuming  fixed  values  for  entrainments 

set  equal  to  1 whenever  Internal  recycle  loop 
for  stage  i does  not  converge 

index  used  in  calculation  of  extent  of  nitration 
in  each  stage 

number  of  passes  through  SUBROUTINE  VESSEL 
required  tr  reach  a converged  solution 

index  used  in  oaloulatlon  of  extent  of  nitration 
in  each  stage 

printout  switch . If  equal  to  1 printout  of  phase 
compositions , gas  flows , vessel  time  constants , 
volumetric  flowe , acid  phase  solubility,  vessel 
holdups  and  stream  compositions  is  obtained 

material  balance  switch.  If  equal  to  1 chemical 
element  balances  are  carried  out  on  the  entire 
nitration  section . each  individual  nitration  stage , 
and  each  individual  vessel 

nitrator  number 

vessel  index.  Has  value  of  1 if  vessel  is  a nitrator 
and  0 if  vessel  is  a separator 

intermediate  printout  switch . If  equal  to  1 a 
printout  of  external  recycle  molar  component 
flows  after  each  iteration  will  be  given 


1 

j 


I 


1 


I 

I 


1 
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IPNCH 


IPRNT 


ISTOE 

ISTOP 

ITYPE 
IV 
IV 1 


IX 

lY 

J 

K 

KEY 

L 


punch  switch.  If  equal  to  zero  no  punched  output. 
If  equal  to  1 final  steady  state  (i  ,e.  values  of 
component  molar  flows  for  both  phases  in  each 
vessel)  is  punched  out  for  use  as  initial  conditions 
In  dynamiu  simulation  or  new  first  guess  in 
steady  state  simulation 

SUBROUTINE  VESSEL  printout  switch . If  equal 
to  1 gives  a variety  of  debug  printout  which  can 
be  used  to  analyze  problems  associated  with 
the  vessel  equations 

stage  number 

equal  to  number  of  data  sets  plus  1 , Program 
execution  la  terminated  when  NSET  s ISTOP 

equal  to  0 for  a nitrator  and  1 for  a separator 

vessel  number 

vessel  number  used  during  Initial  determination 
of  recycle  streams  from  first  guess  values  of 
component  molar  phase  flows 

index  used  in  calculation  of  extent  of  nitration 
in  each  stage 

index  used  in  calculation  of  extent  of  nitration 
in  each  stage 

general  index  integer  used  for  a variety  of 
counting  purposes 

general  index  integer  used  for  a variety  of 
counting  purposes 

index  which  Indicates  whether  extent  of  alpha  or 
meta  nitration  is  to  be  determined  in  a given  stage 

vessel  index  used  in  holdup  calculation  printout 


I 


M 

MHOLDA 

MHOLO0 

MW  (I) 
NATM 

NC 

NCI 
Nf  2 
NC3 
NC4 
NG6 
NOS 

N9(I) 

NSET 


vessel  index  used  in  holdup  calculation  printout 

total  moles  of  add  phase  holdup  in  the  vessel 
currently  being  calculated 

total  moles  of  organic  phase  holdup  in  the  vessel 
currently  being  calculated 

molecular  weight  of  component  i 

index  of  atomic  species  used  in  elemental 
material  balance  calculations 

maximum  number  of  iterations  allowed  on  any 
given  recycle  (internal  or  external)  convergence 
loop 

running  sum  of  iterations  on  external  recycle 
loop  for  separator  2 

running  sum  of  iterations  on  external  recycle 
loop  for  separator  3 

running  sum  of  iterations  on  external  recycle 
loop  for  separator  4 

running  sum  of  iterations  on  external  recycle 
loop  for  separator  6 

running  sum  of  Iterations  on  external  recycle 
loop  for  separator  6 

number  of  sets  of  operating  conditions  for  which 
F steady  state  solution  is  desired  in  a given 
simulation  run 

notch  setting  on  the  Internal  recycle  gate  valve 
for  separator  i 

number  of  the  data  set  for  which  the  steady  state 
solution  is  currently  being  determined 


185 


NTOS  (I) 

NTQSIA 

NTCSS3A 

OUTN 

FACID  (I.J) 

PAS  (I) 

PCTAS 

PCTD 

PCTNB 

PCTNBD 

PF(I.J) 

PFLBS(I.J) 


molar  ratio  of  nitric  acid  to  the  sum  of  sulfuric 
acid  and  SO,  in  the  external  recycle  stream 
flowing  from  the  ith  separator 

molar  ratio  of  nitric  acid  to  the  sum  of  sulfuric 
acid  and  SO|  in  the  acid  phase  flowing  from 
nltrator  lA 

same  as  NTOSIA  but  for  nltrator  3A 

moles/hr  of  ring  associated  NO|  groups  flowing 
out  of  a given  nitration  stage 

pounds  of  component  i held  up  in  the  acid  phase 
of  vessel  ] 

moles/hr  of  component  i in  the  organic  phase 
flowing  from  the  after  separator 

weight  % dissolved  nitrobody  in  the  acid  phase 
flowing  from  the  after  separator  under  the  condi- 
tion of  1000  Ib/hr  of  yellow  water  added  to  the 
after  separator 

weight  % ONT  in  the  organic  components  dissolved 
in  the  acid  phase  flowing  from  separator  1 

weight  % total  nitrobody  contained  in  the  acid 
stream  flowing  from  separator  1 to  the  after 
separator 

weight  % dissolved  nitrobody  contained  in  the 
acid  stream  flowing  from  separator  1 to  the 
after  separator 

moles/hr  of  component  1 in  the  organic  phase 
of  the  organic  stream  flowing  from  separator  j 

same  as  PP  but  in  units  of  Ib/hr 
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K.) 

PFl(I) 

moles /hr  of  component  1 in  the  organic  phase  of 
the  organic  stream  flowing  from  seprtrator  1 

PF2  (I) 

same  as  PFl  but  for  separator  2 

PF3(I) 

same  as  FFl  but  for  separator  3 

PF4(1) 

same  as  FFl  but  for  separator  4 

PF6(I) 

same  as  PFl  but  for  separator  6 

FF6  (I) 

same  as  PFl  but  for  separator  6 

PI  (I) 

total  organic  phase  moles/ hr  of  component  1 
flowing  Into  a nltrator  which  Is  followed  by 
another  nltrator 

PII(I) 

total  organic  phase  molea/hr  of  component  1 
(less  organic  phase  of  internal  recycle)  flowing 
into  a nltrator  which  is  followed  by  another 
separator 

PIR(I,J) 

moles/hr  of  component  1 in  the  organic  phase  of 
the  Internal  recycle  stream  flowing  from 
separator  J 

PIRLB8(I,J) 

same  as  FIR  but  in  units  of  Ib/hr 

PIRl  (I) 

molea/hr  of  component  1 in  the  organic  phase 
of  the  internal  recycle  stream  flowing  from 
separator  1 

PIR2  (I) 

same  as  FIRl  but  for  separator  2 

PIR3  (I) 

same  as  FlRl  but  for  separator  3 

PIR4  (I) 

same  as  FIRl  but  for  separator  4 

PIRSd) 

same  as  FIRl  but  for  separator  5 

PIR6(I)  same  as  PlRl  but  for  separator  6 


POUNDS  (I, J) 


PV(I.J) 


PVI(I.J) 


PXR(I,J) 


weight  % nltrosylsulfurlc  acid  in  the  acid  phase 
of  the  external  recycle  stream  flowing  from  a 
given  separator 

pounds  of  component  1 held  up  in  the  organic 
phase  of  vessel  j 

total  moles/hr  of  organic  phase  flowing  from 
the  vessel  for  which  the  holdup  is  being  computed 

moles/hr  of  component  1 in  the  total  organic 
phase  flowing  from  vessel  J 

first  guess  at  converged  value  of  moles/hr  of 
component  1 in  the  organic  phase  flowing  from 
vessel  J 

molea/hr  of  component  1 in  the  organic  phase  of 
the  external  recycle  stream  flowing  from  separator  J 


PXRLBS  (I.J) 


same  as  PXR  but  in  units  of  Ib/hr 


PXRl(I) 


PXR  2 (I) 
PXR20(I) 


PXR3(I) 
PXR30(I) 
PXR4  (I) 
PXR49(I) 
PXR5  (I) 


PXR50(I) 


moles/hr  of  component  i in  the  organic  phase  of 
the  external  recycle  strr--'  flowing  from  separator  1 

same  as  PXRl  but  for  separator  2 

previous  value  of  PXR2(I)  to  which  current 
value  is  compared  during  the  external  recycle 
convergence  procedure 

same  as  PXRl  but  for  separator  3 

same  as  PXR2Q  but  for  separator  3 

same  as  PXRl  but  for  separator  4 

same  as  PXR20  but  for  separator  4 

same  as  PXRl  but  for  separator  5 

same  as  PXR 20  but  for  separator  S 
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PXR6  (I) 

PXR0O(I) 

PIA(I) 

PlAI 
PIB(I) 
PIBI 
PIS  (I) 

P2(l) 

P2S  (I) 
P3A(1) 
P3B  (I) 
P3S(I) 

P4  (I) 
P4S(I) 
PS  (I) 

P6S  (I) 

P8(I) 

P6S(I) 

Q 


same  as  PXRl  but  for  separator  6 

same  as  PXR20  but  for  separator  6 I 

moles/hr  of  component  i in  the  total  organic  | 

phase  flowing  from  nitrator  lA  H 

mole  a/ hr  toluene  to  nitrator  lA  | 

1 

I 

same  as  PIA  but  for  nitrator  IB 

7 

! 

molea/hr  toluene  to  nitrator  IB  | 

moles/hr  of  component  i in  the  total  organic  i 

phase  flowing  from  separator  2 

moles/hr  of  component  1 in  the  total  organic 
phase  flowing  ffom  nitrator  2 

same  as  PIS  but  for  separator  2 | 

same  as  F2  but  for  nitrator  3A 

same  as  P2  but  for  nitrator  3B  | 

same  as  PIS  but  for  separator  3 | 

i 

same  as  P2  but  for  nitrator  4 

1 

same  as  PIS  but  for  separator  4 

same  as  P2  but  for  nitrator  5 

same  as  PIS  but  for  separator  5 

same  as  P2  but  for  nitrator  6 

same  as  PIS  but  for  nitrator  6 

accumulation  term  used  in  calculation  of  nitrator 
heat  loads 


t 
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QA 

ft* /hr  of  total  acid  phase  flowing  from  the  vessel 
for  which  holdup  is  being  computed 

QP 

same  as  QA  but  foi  organic  phase 

QQ(l) 

fraction  nitric  acid  converted  to  nitronium  ion 
in  vessel i 

QT 

ft*  /hr  of  total  material  flowing  from  the  vessel 
for  which  holdup  is  being  computed 

R(l) 

rate  of  reaction  for  reaction  1 In  the  vessel 
currently  being  simulated  by  SUBROUTINE 
VESSEL 

RATES  (I,  J) 

rate  of  reaction  for  reaction  j in  nitrator  1 

RECN 

molas/hr  of  ring  associated  NO|  groups  flowing 
into  a nitration  stage  with  the  external  recycle 
from  the  next  higher  stage 

RHO  (0 

molar  density  of  component  1 

RHQA 

molar  density  of  the  acid  phase  in  the  vessel  for 
which  holdup  is  being  calculated 

RHQNBS  (1) 

mass  density  of  the  nitrobody  dissolved  in  the 
acid  phase  of  the  external  recycle  stream  flowing 
from  the  1^^  separator 

RHQtP 

molar  density  of  the  organic  phase  in  the  vessel 
for  which  holdup  is  being  calculated 

RHQIl 

equivalent  to  RHQ(l) 

RH(?10 

equivalent  to  RH0(1O} 

RHOll 

equivalent  to  RHQ(ll) 

RHQ12 

equivalent  to  RHQM12) 

RHQ13 

equivalent  to  RH0(13) 
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RH014 

equivalent  to  RHQ  (14) 

RHQ2 

equivalent  to  RH0(2) 

RHQ3 

equivalent  to  RHQ(3) 

RHQ4 

equivalent  to  RHQ(4) 

RHQS 

equivalent  to  RHO  (5) 

RHQ6 

equivalent  to  RHG  (6) 

RH07 

equivalent  to  RHQ(7) 

RH08 

equivalent  to  RHQ  (8) 

RH09 

equivalent  to  BHO  (9) 

RP 

ooeffiolent  of  the  molecular  oxygen  term  in 
reaction  9 

RPP 

precomputed  constant  based  on  RP  which  is  used 
in  material  balance  calculations 

RR 

ideal  gas  law  constant 

SA 

adjustable  constant.  NOT  USED 

SB 

equivalent  to  KMAP 

SC 

equivalent  to  AKEQ4  in  SUBROUTINE  VESSEL 

SD 

adjustaV'le  constant.  NOT  USED 

SE 

equivalent  to  AK9.  NOT  USED 

SO 

equivalent  to  E9.  NOT  USED 

SH 

equivalent  to  AKEQl  in  SUBROUTINE  VESSEL 

SI 

adjustable  constant.  NOT  USED 
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SJ 


equivalent  to  AKEQ3  in  SUBROUTINE  VESSEL 


SNA  (I)  net  moles/hr  of  ring  associated  alpha  NO,  groups 

flowing  through  the  iih  nitration  stage.  Extent 
of  alpha'-nltration  in  nitration  stage  i 

SNM(I)  same  as  SNA  but  for  meta-nitration 


i 

\ 

I" 

i: 


L 


SOL  (I) 

SPA 

SPACID 

SPP 

STN 

SULP 

SUM 

SUMA 

SUMI 

SUMO 


saturation  solubility  of  total  organic  con)ponents 
in  the  acid  phase  of  vessel  i 

mass  flow  rate  of  the  acid  phase  of  the  external 
recycle  stream  (spent  acid)  from  separator  1 

mass  flow  rate  of  the  total  external  recycle  stream 
from  separator  1 

mass  flow  rate  of  the  organic  phase  of  the 
external  recycle  stream  from  separator  1 

total  flow  rate  (Ib/hr)  of  strong  nitric  acid  into 
the  process  from  all  raw  material  feed  streams 

total  flow  rate  (Ib/hr)  of  oleum  into  the  process 
from  all  raw  material  feed  streams 

total  molar  flow  rate  of  gasses  evolved  in  the 
nitration  stage  for  which  material  balance  calcu- 
lations are  being  executed.  Also  the  moles/hr 
of  organics  flowing  from  separator  6 with  the 
organic  stream 

mass  flow  rate  acid  components  of  the  acid  phase 
from  the  after  separator  when  the  after  separator 
option  is  used 

accumulation  term  for  input  streams  used  during 
material  balance  calculations  (moles/hr) 

mass  flow  rate  of  organic  components  In  the  acid 
phase  from  the  after  separator  when  the  after 
separator  option  is  used 


i 

i 

k 

'll 

I 
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Il‘ 

?' 


..X.  .. . 


S15 

total  molea/hr  of  COi.s  evolved  in  nitration 
section 

S16 

total  moles/ hr  of  NO  evolved  in  the  nitration 
section 

S18 

total  moles/hr  of  TNM  evolved  in  the  nitration 
section 

w 


I 


! !: 


S19 

TANB 

TAS  (I) 
TATMI 

TATMQ 

TAU 

TC<3)X 

TC02 

TDNT 

TDNTA 

TDNTP 


total  moles/hr  of  C0|  evolved  in  the  nitration 
section 

total  moles/hr  of  dissolved  organics  in  the  acid 
phase  of  the  external  recycle  flowing  from 
separator  1 

total  acidity  of  the  acid  phase  of  the  external 
recycle  stream  flowing  from  separator  i 

total  moles/hr  of  a given  atomic  specie  entering 
that  section  of  the  process  over  which  an  elemental 
material  balance  is  being  computed 

total  rnoles/hr  of  a given  atomic  specie  exiting 
that  section  of  the  process  over  which  an  elemental 
material  balance  is  being  computed 

holdup  time  in  a given  vessel 

equivalent  to  S15 

equivalent  to  S19 

to^al  mass  flow  rate  of  DNT  in  the  organic 
stream  from  separator  6 

total  pounds  of  DNT  held  up  in  the  acid  phase 
of  all  process  vessels 

total  pounds  of  DNT  held  up  in  the  organic  phase 
of  all  process  vessels 
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TEMP  (I) 
TITLE  (I) 

TLBS 

TMQLES 
TN  (1) 

TNaX 

TNTI 

TOL 

TPNB 

TP(2UND 


TS(I} 

TSEP  (I) 
TSEPl 


temperature  of  the  nltrator 

alphanumeric  identifying  information  for  a given 
data  set 

total  lb /hr  of  acid  components  in  the  acid  phase 
of  the  external  recycle  stream  flowing  from  the 
separator  for  which  acid  composition  calcula- 
tions are  being  executed 

same  as  TLBS  but  in  units  of  moles/hr 

weight  % total  nitric  aoid  in  the  acid  phase  of  the 
external  recycle  stream  flowing  from  the  1th 
separator  on  an  organic  free  basis 

equivalent  to  S16 

moles/hr  of  TNT  entering  the  nitration  section 
with  feed  material 

Ib/hr  toluene  fed  to  the  nitration  section 

total  moles/hr  of  organic  components  in  the 
organic  phase  of  the  external  recycle  stream 
flowing  from  separator  1 

accumulation  term  for  total  pounds  of  organic 
material  held  up  in  either  the  acid  or  organic 
phase  of  the  vessel  for  which  holdup  calculations 
are  being  executed 

reference  temperature 

weight  % total  sulfuric  acid  plus  SOs  (as  equivalent 
sulfuric)  in  the  acid  phase  of  the  external  recycle 
stream  flowing  from  the  iih  separator  on  an  organic 
free  basis 

temperature  of  the  ith  separator 
equivalent  to  TSEP  (1) 


J 
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TSEP2 

TSEP3 

TSEP4 

TSEP5 

TSEP6 

TTDNT 

TTLBS  (I) 

TTLBSA  (I) 

TTNM 

TTNTA 

TTNTP 

TTTLBA 

TTTLBS 

TTTNT 

TTTTLB 

TVO»L 


equivalent  to  TSEP  (2) 
equivalent  to  TSEF  (3) 

equivalent  to  TSEP  (4) 
equivalent  to  TSEP  (6) 
equivalent  to  TSEP  (6) 

total  pounda  of  DNT  held  up  In  all  process  vessels 

total  pounds  of  organic  phase  held  up  In  the 
1th  vessel 

total  pounds  of  organic  components  held  up  in 
the  acid  phase  of  the  ith  vessel 

equivalent  to  S18 

total  pounds  of  TNT  held  up  In  the  acid  phase 
of  all  process  vessels 

total  pounds  of  TNT  held  up  In  the  organic  phase 
of  all  process  vessels 

total  pounds  of  organic  components  held  up  In 
the  acid  phase  of  all  process  vessels  between 
and  Including  nltrator  2 and  separator  6 

total  pounds  of  organic  phase  held  up  In  all 
process  vessels  between  and  including  nltrator  2 
and  separator  6 

total  pounds  of  TNT  held  up  In  all  process  vessels 

total  pounds  of  organic  components  (In  both  acid 
and  organic  phases)  held  up  In  all  process  vsssels 
between  and  including  nltrator  2 and  separator  6 

volumetric  flow  of  the  organic  compounds  In  the 
acid  phase  of  the  external  recycle  stream  flowing 
from  a separator  for  which  nitrobody  density  Is 
being  computed 
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TIA 

temperature  of  nitrator  lA 

TIB 

temperature  of  nltraior  IB 

T2 

temperature  of  nitrator  2 

T3A 

temperature  of  nitrator  3A 

T3B 

temperature  of  nitrator  SB 

T4 

temperature  of  nitrator  4 

T6 

temperature  of  nitrator  5 

T6 

temperature  of  nitrator  6 

V(I) 

atorage  vector  for  current  valuei  of  the  1 com- 
ponent flow  terma  generated  in  SUBROUTINE 
VESSEL 

VHQLDA 

volume  of  acid  phase  held  up  in  a given  veeeel 

VHOLOQ 

volume  of  organic  phase  held  up  in  a given  vessel 

VKM(I) 

working  volume  of  the  1^^  nitrator 

VKMS 

working  volume  of  any  separator 

VOL  (I) 

working  volume  of  the  llh  vessel 

WKN 

total  flow  rate  (Ib/hr)  of  weak  nitric  acid  into 
the  process  from  all  raw  material  feed  stream 

XP(I,J) 

mole  fraction  of  component  1 in  the  organic 
phase  within  vessel  j 

XPA  (I,J) 

mole  fraction  of  component  1 in  the  acid  phase 
within  vessel  J 

YELLOW  (I) 

moles/hr  of  component  1 in  the  yellow  water 
fed  to  the  after  separator  whenever  the  after 
separator  option  is  chosen . 
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YIELD  moles  a -TNT  In  the  organic  stream  flowing  from 

separator  6 per  mole  of  toluene  fed  to  the  process 
(expressed  as  %) 

ZA  temperature  coefficient  in  the  exponent  of  the 

Arrhenius  rate  law 

Subroutine  N1T8BP 

The  definitions  of  the  following  variables  are 
Identical  to  those  In  the  main  program:  AP  (I,  J) , API  (1) , AI  (1) , All  (I) , 
AIR  (I . J) , AIRl  (I) , AV  (I , J) , AXR  (I) , AXRl  (I) , DEL , EAQ  (I) , BAQIP  (D  . 
BOA  (I) , EQAP  (1) , FACD  (1) . 015  (I) . 016  (1) . 018  (I) , 018  (I) , I , lAFT . 
IGHECK,  IB,  IPLO(I) , IN,  IPRNT,  ISTQE,  ITYPB,  IV,  NOd) . PF(1,J) , 
PFl(I),  PI(I),  PII(I),  P1R(I,J),  PIRl(l).  PV(I.J),  PXR(I),  PXRKIl, 
RATES  (I,J) , SQL  (I) , TSEP  (I) , YELLOW  (I) , 


Remaining  variables  are  defined  as  follows: 


AIRQ  (I) 

previous  value  of  AIR(l.J)  to  which  current 
value  is  compared  during  the  internal  recycle 
convergence  procedure  for  stage  J 

ICOUNT 

counter  for  iterations  of  internal  reoyole 
convergence  loop 

IVA 

vessel  index  for  separator  1 used  whenever 
after  separator  is  included  in  the  simulation 

IVQ 

vessel  index  for  the  nitrator  which  preoeeds 
the  separator  for  which  internal  recycle 
convergence  is  being  executed 

PIRQ(I) 

previous  value  of  PIR(I,J)  to  which  current 
value  is  compared  during  the  internal  recycle 
convergence  procedure  for  stage  j 
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Subroutine  VESSEL 
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The  definitions  of  the  following  variables  are  identical  to  those  in  the 
main  program;  AF(I.J) , AI(I) , AIR(l.J) , AKA,  AKEQ,  AKEQOV,  AKMA, 
AKMAP.  AKMATL,  AKIO,  AK12,  AK3S.  AK3W,  AK4S,  AK4W.  AK48, 

AK6S.  AK7,  AK8.  AV(I.J).  AXR(I.J),  CP  (I).  DEL,  DH(I),  EA^d), 

EAaP(I),  BOA(l),  BQAP(I).  ElO,  E12,  E3S,  E3W,  B4S,  E4W,  E6S,  BBS, 

E7,  EB,  FACDd.J),  PP.  FMAX,  015 d).  OlBd).  018d),  018d).  I.  lAFT, 
ICHBCK,  IE.  Ill,  IN,  IPRNT,  ISTOB,  ITYPE,  IV,  J,  MNd) , NOd) . PFd.J) , 
Pld),  PIRd.J).  PVd.J).  PXRd.J),  QQd).  Rd).  RATES d.J) , RH91.  RMOIO. 
RH912,  RH013,  RHQ14,  RH92,  RH93.  RH94,  RH95,  RHQB,  RH97,  RH98. 
RH9B.  RP.  RR.  8A.  SB,  SC,  SD.  SB,  SO,  SH,  81,  SJ,  SOLd).  TEMPd), 

TR,  TSBPd) , V(I) , VKMd) , VKM8,  YELLOW  d) 

The  definitions  of  the  following  variables  are  identical  to  other  variables 
defined  in  the  main  program: 

Variable  in  VESSEL  Variable  in  main  program 

All  thru  AI17 
Fid)  thruFlOd) 

ISEP 

PIl  thru  PI17 
R1  thru  RIO 

Remaining  variables  are  defined  as  follows: 

Ad)  moles/hr  of  component  i in  the  acid  phase  flowing 

from  the  vessel  currently  being  simulated 

AINB  moles/hr  of  organic  components  in  the  combined 

acid  phase  flowing  to  the  vessel  currently  being 
simulated 

AKEQl  equilibrium  constant  for  the  dissociation  of 

nitric  acid  in  aqueous  sulfuric  acid 

AKEQ 2 equilibrium  constant  for  the  dissociation  of 

sullfuric  acid  in  an  aqueous  system 

AKSQ3  equiliblrum  constant  for  dissociation  of  nitric 

acid  in  anhydrous  sulfuric  acid 
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Aim  thru  AI(17) 
Old)  thru  OlOd) 
ISTOE 

PI(1)  thru  PI (17) 
R(l)  thru  R(IO) 





AKEQ4 


equilibrium  constant  for  the  recombination  of 
nitronium  and  bieulfate  ions  in  anydrous  acid 
mixtures 


AKl 

computed  value  of  AKBQl  obtained  from  Interval 
halving  procedure 

AK3 

computed  value  of  AKBQ3  obtained  from  Interval 
halving  procedure 

AK4 

computed  value  of  AKBQ4  obtained  from  Interval 
halving  procedure 

AL 

molea/hr  water  in  the  acid  phase  flowing  from 
the  vessel  currently  being  simulated;  negative 
values  represent  moles/hr  SO| 

AM(1) 

moles/hr  of  component  i transferred  from  the 
organic  phase  into  the  sold  phase  by  bulk  mass 
transfer  in  the  vessel  currently  being  simulated. 
Used  in  check  calculations  only 

AMNB 

total  moles/hr  of  organic  components  transferred 
from  the  organic  to  acid  phase  by  bulk  mass 
transfer  in  the  vessel  currently  being  simulated 

AMI 

moles/hr  o -MNT  transferred  from  the  organic 
phase  into  the  acid  phase  by  bulk  mass  transfer 
in  the  vessel  currently  being  simulated 

ANB 

moles/hr  organic  components  in  the  acid  phase 
flowing  from  the  vessel  currently  being  simulated 

AS04 

combined  moles/hr  of  H|SO«  and  SO|  in  the  acid 
phaae  flowing  from  the  vessel  currently  being 
simulated 

AT 

total  molos/hr  of  acid  phaso  flowing  from  the 
vessel  currently  being  simulated 

ATI 

reciprocal  of  AT 
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AXL 

same  aa  AL  but  uaed  only  in  check  calculationa 

AXl 

1.  + ClA-CA 

AXl  thru  AX8 

■ante  aa  A1  thru  AS  but  uaed  only  in  check 
calculationa 

AX12 

aame  aa  A12  but  uaed  only  in  check  calculationa 

A1 

molea/hr  a MNT  in  the  acid  phaae  flowing  from 
the  vaaaal  currently  being  aimulated 

A2 

aame  aa  Al  but  for  m MNT 

A3 

aame  aa  Al  but  for  a DNT 

A4 

aame  aa  Al  but  for  m MNT 

A5 

aame  aa  Al  but  for  a TNT 

A6 

aame  aa  Al  but  for  m TNT 

A7 

aame  aa  Al  but  for  TNBX 

A8 

aame  aa  Al  but  for  TNB 

A9 

molea/hr  nitric  acid  in  the  acid  phaae  flowing 
from  the  veaael  currently  being  aimulated 

All 

aame  aa  AS  but  for  H|SO( 

A12 

aame  aa  AS  but  for  nltroaylaulfuric  acid 

A13 

aame  aa  AS  but  for  water 

A14 

aame  aa  AS  but  for  SO| 

A17 

molea/hr  nitronium  ion  in  the  acid  phaae  flowing 
from  the  veaael  currently  being  aimulated 

B 

free  aulfuric  acid  concentration  in  oleum  in  the 
veaael  currently  being  aimulated 

( 


[■ 


BC 


B (3  X C) 


B1 

B2 

C 

CA 

CB 

CF 

CHNela 

CHSfif4 

CHS2gi7 

CH2Q 

CH2804 

CH30 

CN02 

CR 


term  representing  moles/hr  of  an  organic 
component  flowing  Into  a vessel  with  the  combined 
acid  phase  plus  the  amount  of  that  component 
formed  by  reaction 

moles/hr  of  an  organic  component  flowing  Into 
a vessel  with  the  combined  organic  phase 

concentration  of  pyrosulfurlo  acid  In  the  vessel 
currently  being  simulated 

ratio  of  AMNB  to  ANB  (constant  used  In  determ- 
ining bulk  mass  transfer  during  slmulataneous 
solution  of  vessel  equations) 

(.5  X B)  ^ (1.6  X C) 

C - F 

equilibrium  concentration  of  nitric  acid  in  the 
acid  phase  of  the  vessel  currently  being  simulated 
(moles/of*) 

same  as  CHN03  but  for  blsulfate  Ion 

same  ae  CHN03  but  for  pyrosulfate  Ion 

same  as  CHN(2f3  but  for  water 

same  as  CHN93  but  for  sulfuric  acid 

same  as  CHN03  but  for  hydronlum  Ion 

same  as  CHN03  but  for  nltronlum  Ion 

ratio  of  AMNB  to  PNB.  Constant  used  In  determ- 
ination of  bulk  mass  transfer  during  simultaneous 
solution  of  vessel  equations 


CV0L 


computed  volumetric  flow  of  an  internal  recycle 
stream  based  on  entrainment  correlation 


cxx 


covrelation  used  to  determine  fraction  organic 
entrained  in  the  internal  recycle  stream 


ratio  of  fraction  Internal  recycle  to  fraction 
external  recycle  used  in  entrainment  correlation 

equal  to  overall  mass  transfer  coefficient  (ETA) 
during  elmultaneous  solution  of  vessel  equations 

coefficient  equal  to  Cl/ANB  used  during  simul- 
taneous solution  of  vessel  equations 

coefficient  equal  toCl/PNB  used  during  simul- 
taneous solution  of  vessel  equations 

coefficient  related  to  toluene  diffusion  which  is 
used  during  simultaneous  solution  of  vessel 
equations 

coefficient  equal  to  acid  phase  mass  transfer 
plus  the  sum  of  the  reaction  coeffioients  used 
during  simultaneous  solution  of  the  vessel 
equations 

coefficient  equal  to  (-  C10  - CR)  used  during 
elmultaneous  solution  of  the  vessel  equations 

coefficient  equal  to  (CA  - CIA)  used  during 
simultaneous  solution  of  the  vessel  equations 

coefficient  equal  to  ( 1 ciQI  CR)  used  during 
simultaneous  solution  of  the  vessel  equations 

equal  to  the  sum  of  reaction  rates  R1  thru  R6. 
Represents  the  moles/hr  water  produced  by 
reactions  1 thru  6 

moles/hr  nltrosylsulfuric  acid  produced  by 
oxidation  reactions 


net  moles/hr  nltrosylsulfuric  acid  produced 
after  decomposition  takes  place 
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C6 

organic  phase  nitric  acid  solubility  coefficient 

DD(I) 

rate  of  diffusion  of  organic  component  1 . Used 
only  in  oheok  calculations 

DLA 

oonvargenoe  tolerance  for  components  in  both 
acid  and  organic  phases  flowing  from  the  vessel 
currently  being  simulated 

DRM 

difference  in  molar  densities  between  the  acid 
and  organic  phases  used  in  the  entrainment  nor- 
relation 

01 

molea/hr  of  a MNT  which  diffuse  from  the  organic 
into  the  acid  phase  in  the  vessel  currently  being 
simulated.  Negative  values  means  diffusion  is 
from  acid  to  organic  phase 

OlO 

molee/hr  of  toluene  which  diffuse  from  the  bulk 
organic  phase  to  the  acid-organic  interface 

02 

same  as  D1  but  for  mMNT 

03 

earns  as  D1  but  for  o DNT 

04 

earns  as  D1  but  for  mDNT 

OS 

same  as  Dl  but  for  a TNT 

06 

sime  as  Dl  but  for  mTNT 

07 

same  as  Dl  but  for  TNBX 

08 

same  ss  Dl  but  for  TNB 

BA94 

molar  ratio  of  acid  phase  in  the  organic  stream 
to  total  organic  stream , based  on  the  entrainment 
correlation 

BTA 

overall  mass  transfer  coefficient  in  the  veeeel 
currently  being  simulated 

I 
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ETAA 

ETAP 

BTATL 

BXX 

P 

FA 

FACT 

PAR 

FARl 

FC 

FF(I) 


] 

acid  phase  mass  transfer  coefficient  in  the  ' ' O \ ' 

currently  being  simulated 

organic  phase  mass  transfer  coefficient  in  the 
vassal  currently  being  simulated 

mass  transfer  coefficient  for  toluene  In  the 
organic  phase 

ratio  acid  phase  to  organic  phase  In  the  organic 
atream , baaed  on  the  entrainment  correlation 

Initial  oonosntratlon  of  nitric  acid  prior  to 
aatabllahment  of  equilibrium  oonditiona  In  the 
vessel  currently  being  simulated 

volume  fraction  acid  phase  In  the  «;easel 
currently  being  simulated 

relaxation  factor  used  in  convergence  procedure 
for  composition  of  phases  flowing  from  the  vassal 
currently  being  simulated 

fraction  of  the  total  acid  phase  flowing  from  a 
separator  that  Is  raoyolad  for  the  separator 
currently  being  simulated 

fraction  of  the  total  acid  phase  flowing  from  a 
separator  that  Is  entrsined  In  the  organic  stream 
which  flows  to  the  next  nitrator  for  the  separator 
currently  being  simulated 

constant  equal  to  (2  x F)  - C used  in  computation 
of  nltronium  Ion  equilibrium  concentration  in 
region  2 for  the  vessel  currently  being  simulated 

cross  sectional  area  fraction  corresponding  to  the 
1th  notch  setting  on  the  gbte  valve  in  the  internal 
recycle  lino  of  the  separator  currently  being 
simulated 
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FOQC 


FP91 


FRAC 


FVKM 


GAMMA 


upper  boundary  value  of  Internal  recycle 
fraction  computed  uaing  the  entrainment 
coefficient 

internal  recycle  fraction  computed  uaing  the 
entrainment  correlation  for  the  aeparator 
currently  being  aimulated 

aame  aa  FHQ  except  lower  boundary  value 

molar  ratio  of  organic  phaae  in  the  acid  recycle 
atreama  to  total  organic  phaae  flowing  from  the 
aeparator  currently  being  aimulated  (baaed  on 
entrainment  correlation) 

current  value  of  F(?G 

volume  fraction  of  organic  phaae  in  the  veaael 
currently  being  aimulated 

molar  fraction  of  the  total  organic  phaae  flowing 
from  a aeparator  which  appeara  In  the  organic 
atraam  for  the  aeparator  currently  being  aimulated 

molar  fraction  of  the  total  organic  phaae  flowing 
from  a aeparator  which  appeara  in  the  recycled 
acid  atreama  for  the  aeparator  currently  being 
aimulated 

volume  fraction  of  the  combined  recycled  acid 
atreama  which  la  external  recycle  for  the  aeparator 
currently  being  simulated 

volume  fraction  of  the  combined  rncycled  acid 
atreama  which  la  Internal  recycle  for  the  aeparator 
currently  being  simulated 

volume  of  the  vessel  currently  being  simulated 

acid  concentration  factor  for  all  nitration  reactions 
in  the  vessel  currently  being  aimulated 


G1 


OlO 

OlOP 

015P 

03 
03D 

03S 

04 
04D 

04S 

05 
050 

06 
060 


equivalent  to  FI  (I)  where  I Is  the  Index  of  the 
nltrator  currently  being  simulated  or  the  index 
of  the  nitrator  preceding  the  separator  currently 
being  simulated 

same  as  01  but  for  FIO  (I) 

rata  of  reaction  10  per  mole  of  TNT  in  the  acid 
phase  flowing  from  the  vessel  currently  being 
simulated 

equivalent  to  015  (1)  where  I is  the  index  of 
the  vessel  currently  being  simulated 

same  as  01  but  for  F3(I) 

rate  of  reaction  3 per  mole  of  a MNT  in  the  acid 
phase  flowing  from  the  vessel  currently  being 
simulated 

same  as  01  but  for  F3S  (I) 

same  as  01  but  for  F4(I] 

rate  of  reaction  4 per  mole  of  mMNT  in  the  acid 
phase  flowing  from  the  vessel  currently  being 
simulated 

same  as  01  but  for  F4S  (I) 

same  as  01  but  for  F5(I] 

rate  of  reaction  5 per  mole  of  a DNT  in  the  acid 
phase  flowing  from  the  vessel  currently  being 
simulated 

same  as  01  but  for  FB(I) 

rate  of  reaction  6 per  mole  of  mDNT  in  the  acid 
phase  flowing  from  the  vessel  currently  being 
simulated 
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O G7 

G7G 

08 

G8P 

09 

G9G 

ICOUNT 

ITH 

MWAVGA 

MWAVGP 

NC 

NOTCH 

P(I) 


same  as  G1  but  for  F7  (I) 

rate  of  reaction  7 per  mole  of  DNT  in  the  acid 
phase  flowing  from  the  vessel  currently  being 
simulated 

same  as  01  but  for  F8(I) 

rate  of  reaction  B per  mole  of  DNT  in  the  acid 
phase  flowing  from  the  vessel  currently  being 
simulated 

same  as  01  but  for  FB  (I) 

:(ero 

number  of  passes  made  through  suoroutlne 
VESSEL  during  a particular  convergence  iteration 
for  the  vessel  currently  being  simulated 

number  of  convergence  iterations  executed  for 
the  vessel  currently  bein^j  simulated.  Five 
Iterations,  each  consisting  of  NC  passes  through 
VESSEL  are  allowed 

average  molecular  weight  of  the  acid  phase 
flowing  from  the  separator  which  is  being  simu- 
lated and  for  which  entrainment  is  being  determ- 
ined via  the  correlation 

same  as  MWAVGA  but  for  the  organic  phase 

maximum  number  of  passes  through  VESSEL 
allowed  per  iteration  during  vessel  equation 
convergence  procedure 

value  of  NO  vl)  where  I is  the  index  of  the  separator 
currently  being  simulated 

moles/hr  of  component  i in  the  organic  phase 
flowing  from  the  vescel  currently  being  simulated 
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PNB 


moles /hi  of  combined  organic  components  in 
the  organic  phase  flowing  from  the  vessel 
currently  being  simulated 


PP 

PT 

PTI 

PX(I) 

PXIO 

PX9 

PI 

PIO 

P2 

P3 

P4 

P5 

P6 

P7 

PS 


exponent  of  the  nltrobody-in-acid  solubility 
correlation 

moles/hr  of  combined  organic  components  in 
both  the  acid  and  organic  phases  flowing  from 
the  vessel  currently  being  simulated 

the  reciprocal  of  PT . Used  for  computation  of 
entrainment  via  the  correlation 

same  as  P (1)  for  components  1 thru  8 but  used 
only  in  check  calculations 

same  as  P (10)  but  used  only  in  check  calculations 

same  as  P(9)  but  used  only  in  check  calculations 

molns/hr  o MNT  in  the  organic  phase  flowing 
from  the  vessel  currently  being  simulated 

same  as  PI  but  for  toluene 

same  as  FI  but  for  mMNT 

same  as  PI  but  for  a DNT 

same  as  FI  but  for  mDNT 

same  as  PI  but  for  o TNT 

same  as  PI  but  for  mTNT 

same  as  PI  but  for  TNBX 

same  as  FI  but  for  TNB 

same  as  PI  but  for  nitric  acid 


rmrn 


O Q 

QA 

QP 

QT 

RHL9 

RH09 

RRl 

RRIO 

RR3 

RR7 

RR8 

RR9 

RR9G 

RTL 

RT9 


•'■'i 


fraction  of  molecular  nitric  acid  converted  to 
nitronium  ion  in  the  veaael  currently  being 
aimulated 

ft*  /hr  of  acid  phaae  flowing  from  the  veaael 
currently  being  aimulated 

ft*  /hr  of  organic  phaae  flowing  from  the  veesel 
currently  being  aimulated 

volumetric  flow  (ft*  /hr)  of  combined  internal 
plui  externa*  recycle  atreama  flowing  from  the 
aeparator  currently  being  aimulated 

maaa  denalty  (Ib/ft* ) of  the  organic  phaae 
flowing  from  the  aeparator  currently  being 
aimulated  and  for  which  entrainment  la  being 
determined  via  the  correlation 

same  as  RHL9  but  molar  denalty  (molea/ft*) 

same  aa  R1  but  used  only  for  check  calculation 

same  aa  RIO  but  used  only  for  check  calculation 

same  aa  RRl  but  for  R3 

same  as  RRl  but  for  R7 

same  ae  RRl  but  for  R8 

same  aa  RRl  but  for  R9 

same  aa  RRl  but  for  R90 

maaa  ratio  of  organic  to  acid  phaae  in  the  aeparator 
currently  being  simulated  and  for  which  entrain- 
ment la  being  calculated  via  the  correlation 

same  as  RTL  but  molar  ratio 
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Ir 


RVQL 


RlOA 


RIOM 


R9QA 


R8QM 


8ANB 


SANBC 


SPNB 


maxlaum  volumetrio  flow  of  Internal  recycle 
Btream  poaalble  for  a given  notch  setting  In  the 
separator  currently  being  simulated  and  for  which 
entrainment  is  being  calculated  via  the  correlation 

rate  of  reaction  10  for  o TNT  In  the  vessel 
currently  being  simulated 

rate  of  reaction  10  for  mTNT  In  the  vessel  cur- 
rently being  simulated 

rate  of  reaction  7 for  mDNT  In  the  vessel  cur- 
rently being  simulated 

rata  of  reaction  8 for  o DNT  In  the  vessel  cur- 
rently being  simulated 

rate  of  raatlon  8 for  mDNT  In  the  vessel  currently 
being  simulated 


Initial  concentration  (moles/ft*)  of  total  sulfuric 
acid  In  the  add  phase  of  the  vessel  currently 
being  simulated 

same  as  ANB  but  used  cnly  In  check  calculations 

same  as  8ANB  but  based  on  computed  solubility 

same  as  PNB  but  used  only  In  check  calculations 

sum  of  the  Individual  organic  component  diffus* 
slons  used  only  In  the  check  calculations 

moles/hr  of  entrained  add  phase  [based  on  the 
correlation)  contained  In  the  organic  stream 
flowing  from  the  separator  currently  being  simulated 


n>olai/hr  of  acid  phase  (baaed  on  the  correlation) 
contained  in  the  internal  recycle  stream  flowing 
from  the  separator  currently  being  simulated 

moloa/hr  of  acid  phase  (baaed  on  the  correlation) 
contained  in  the  combined  internal  and  external 
recycle  atreams  flowing  from  the  separator 
currently  being  eimulated 

molea/hr  of  sold  phase  (baaed  on  the  correlation) 
contained  in  the  external  recycle  stream  flowing 
from  the  aeparator  currently  being  simulated 

temperature  (in  Degrees  Kelvin)  of  the  vessel 
currently  being  simulated 

total  mass  flow  (Ib/hr)  through  the  separator 
currently  being  simulated  and  for  which  entrainment 
is  being  computed  via  the  correlation 

molea/hr  of  organic  phase  (based  on  the  correla- 
tion) contained  in  the  organic  stream  flowing 
from  the  separator  currently  being  simulated 

moles/hr  of  entrained  organic  phase  (based  on 
the  correlation)  contained  in  the  internal  recycle 
stream  flowing  from  the  separator  currently 
being  simulated 

moles/hr  of  entrained  organic  phase  (based  on 
the  correlation)  contained  in  the  combined  recycle 
streams  flowing  from  the  separator  currently 
being  simulated 

moles/hr  of  entrained  organic  phase  (based  on 
the  correlation)  contained  in  the  external  recycle 
stream  flowing  from  the  separator  currently 
being  simulated 


VO  (I) 

storage  vector  for  the  values  of  the  i component 
molar  flows  generator  on  the  previous  pass  through 
VESSEL  for  the  vessel  currently  being  simulated 

w 

concentration  (moles/ft*)  of  water  in  the  vessel 
currently  being  simulated 

X 

computed  equilibrium  concentration  of  nitronium 
ion  in  the  acid  phaae  under  aqueous  conditions 
and  under  conditions  of  oleum  region  3.  Also 
equilibrium  ooncentratlon  at  molsoular  nltrlo  sold 
in  dleum  region  2 

XA 

NOT  USED 

XB 

NOT  USED 

XC 

NOT  USED 

AD 

NOT  USED 

XE 

NOT  USED 

XBQA 

equilibrium  solubility  of  organic  material 
(nltrobody)  in  the  acid  phase  of  the  vessel 
currently  being  simulated 

XF 

NOT  USED 

XL 

^ower  boundary  value  used  In  an  Interval  halving 
convergence  procedure 

XU 

upper  boundary  valued  used  In  an  Interval  halving 
convergence  procedure 

XI 

mole  fraction  a MNT  in  the  aold  phase  o!  tno 
vessel  currently  being  simulated 

Xll 

same  as  XI  but  for  sulfuric  acid 

X12 

same  as  XI  but  for  nltrosylsulfurlc  acid 
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uma  at  XI  but  for  S0| 

■ame  aa  XI  but  for  mMNT 

aama  aa  XI  but  for  a DNT 

aama  aa  XI  but  for  mDNT 

aama  aa  XI  but  for  o TNT 

aama  aa  XI  but  for  mTNT 

aama  aa  XI  but  for  TNBX 

aama  aa  XI  but  for  TNB 

aama  at  XI  but  for  nitric  acid 

aqulllbrlum  oonoantratlon  of  tha  pyroaulfate  ion  in 
tha  add  phaaa  of  tha  vaaaal  baing  aimulatad 
undar  oonditona  of  oleum  region  3 

lower  boundary  value  used  in  an  Interval  halving 
oonvarganca  procedure 

upper  boundary  value  uaed  in  an  Interval  halving 
oonvarganca  prooadurn 

mole  fraction  of  o MNT  in  the  organic  phaae  of 
the  veaael  currently  being  aimulated 

aame  aa  Y1  but  for  toluene 

aama  aa  Y1  but  for  mMNT 

aame  aa  Y1  but  for  a DNT 

aame  aa  Y1  but  for  mDNT 

eame  aa  Y1  but  for  a TNT 
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lame  ai  Y1  but  for  MTNT 

same  oa  Y1  but  for  TNBX 

aama  aa  Y1  but  for  TNB 

lama  aa  Y1  but  for  nitric  acid 

dlffuilon  ooefficient  for  organic  components 
in  the  organic  phase 

diffusion  ooefficient  for  organic  components 
in  the  acid  phase 


Funotlon  BNTH 


The  definitions  of  the  following  varlabla  are  identical  to  those  in  the 
main  program;  AKA,  AKBQ.  AKMA.  AKMATL.  AKIO,  AK12,  AK3S,  AK3W, 
AK4S,  AK4W,  AK6S,  AK6S.  AK7.  AKB,  CP  (I),  OH(l).  ElO.  E21,  E3S. 

B3W,  E48,  E4W,  ESS.  BBS.  E7.  BB,  FD.  FMAX,  RHai,  RHOlO,  RHOll. 
RHai2.  RH013.  RH014,  RH02.  RH93.  RH94.  RHaS,  RHQB,  RHa?,  RH08, 
RH9B,  HP,  HR,  SA,  SB.  SC.  8D.  SB.  SO.  SH.  SI.  SJ.  TR 

Remaining  variables  art  defined  as  follows; 

A(l)  moles/hr  of  component  1 in  the  acid  phase  of 

the  stream  for  which  the  enthalpy  Is  currently 
being  calculated 

CPU  equivalent  to  CP  (11) 


CP13 


CP14 


equivalent  to  CP  (12) 
equivalent  to  CP  (13) 
equivalent  to  CP  (14) 
equivalent  to  CP(B) 

temperature  of  the  vessel  (relative  to  20°C) 
for  which  the  enthalpy  of  a particular  exiting 
stream  is  currently  being  calculated 
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BNTH 

BNTH9 

BNTH2 

F 

I 

J 

LBS 

LB8B 

P(I) 

T 

TP 

WTA 


enthalpy  of  the  stream  currently  being  calculated 
(kcal/hr) 

enthalpy  of  mixing  of  the  acid  components  in 
the  acid  phase  of  the  stream  currently  being 
oaloulated  based  on  the  sero  nitric  curve  of  the 
McKinley  and  Brown  data  (kcal/lb) 

same  as  ENTH9  but  baaad  on  the  20%  nitric  curve 

adjusted  weight  fraction  (nltrobody  free  baala) 
of  nitric  plus  equivalent  sulfuric  acid  in  the  acid 
phase  of  the  stream  currently  being  calculated 

general  index  integer 

general  index  integer 

Ib/hr  of  nitric  acid , sulfuric  acid , water  and 
80|  flowing  with  the  acid  phase  of  the  straatn 
fox  which  the  enthalpy  la  currently  being  calculated 

Ib/hr  of  nitric  acid  plus  equivalent  sulfuric  acid 
in  the  acid  phase  of  the  stream  for  which  the 
enthalpy  is  currently  being  calculated 

moles/hr  of  component  1 in  the  organic  phase 
of  the  stream  for  which  the  enthalpy  Is  currently 
being  calculated 

moles/hr  of  a given  organic  component  in  the 
combined  acid  and  organic  phases  (i.e. , in  the 
total  stream)  of  the  stream  for  which  the  enthalpy 
is  currently  being  calculated 

temperature  of  the  vessel  for  which  the  enthalpy 
of  a particular  exiting  stream  is  being  calculated 

LB8E/LB8;  nltrobody  free  weight  fraction  of  nitric 
plus  equivalent  aulfurio  acid  in  the  acid  phase  of 
the  stream  for  which  the  enthalpy  is  currently 
being  oaloulated 


u 


weight  fraction  of  nitric  add  In  the  nitric  plua 
equivalent  aulfuric  acid  portion  of  the  acid  phase 
of  the  stream  for  which  the  enthalpy  la  currently 
being  calculated 

the  1th  data  point  on  tha  taro  % nitric  curve  from 
the  anthalpy-oonoantratlon  data  of  McKinley  and 
Brown 

tha  1th  data  point  on  tha  20%  nitric  curva  form  tha 
anthalpy-conoantratlon  data  of  McKinley  and  Brown 


CARD 


DATA  READ 


FORMAT 


o 


1 

AK12,  AK3W,  AK4W,  AK88.  AXaS, 
AK8.  AKIO 

AK7, 

8F10.0 

2 

AK38,  AK48.  B38,  E48 

8F10.0 

3 

B21,  B3W.  B4W,  BBS,  B88,  B7.  BB, 

BIO 

6F10.0 

4 

8A.  SB,  8G.  SB.  80.  8H.  81 

BPIO.O 

5 

SJ.  RR.  TR.  AKMA,  AKMATL,  PMAX 
FD,  RP 

8F10.0 

6 

AKA,  AK8,  AKC,  AKBQ,  AXBQOV, 

> 

O 

BPIO.O 

7,  8 

AVI  (1.1).  1-1,14 

8P10.0/8P10.0 

8.10 

PVKl.l).  1-1.14 

6P10.0/8P10.0 

11.12 

AVI  (1.2).  1-1,14 

8F10.0/6F10.0 

13.14 

PV1(1,2),  1-1,14 

8P10.0/ BPIO.O 

18.18 

AVI  (1.3),  1-1,14 

BP10.0/8P10.0 

17,18 

PVKI.3),  1-1,14 

BPIO.O/OPIO.O 

18,20 

AVI  (1.4).  1-1,14 

BPio.o/e* 

21,22 

PVI(1,4),  1-1,14 

8P10.0/6FiL , : 

23,24 

AVI  (1,8),  1-1,14 

BPIO.O/BPIO.O 

26,28 

PVI(1,6),  2-1,14 

BPIO.O/OPIO.O 

27,28 

AVI  (1,6),  1-1,14 

BP10.0/6P10.0 

28,30 

PVK1.6),  1-1,14 

8P10.0/6F10.0 

31,32 

AVI  (1,7),  1-1,14 

8F10.0/6P10.0 

33,34 

PVI(I,7),  1-1,14 

8P10.0/6P10.0 

38,36 

AVI  (1,8),  1-1,14 

BF10.0/6F10.0 

37,38 

PVI(1,8),  1-1,14 

BPIO.O/BPIO.O 

38,40 

AVI  (1, 8).  1-1,14 

BPIO.O/BPIO.O 

41,42 

PVKl.B),  1-1,14 

BPIO.O/BPIO.O 

43,44 

AVI  (1,10),  1-1,14 

BFIO.O/BPIO.O 
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) 


•.V  - •. 


CARD 

DATA  READ 

FORMAT 

45,46 

PVKI.IO),  I»l,14 

8F10.0/6F10.0 

47.48 

AVl(I.ll),  I«l,14 

8F10.0/6F10.0 

48,50 

PVI(I.ll),  1-1,14 

8P10.0/6F10.0 

81,62 

AVI  (1,12),  1-1,14 

8F10.0/6F10.0 

83,84 

PVI(I.l2),  1-1,14 

8F10.0/6F10.0 

68,68 

AVI  (1. 13),  1-1,14 

8P10.0/6F10.0 

67,88 

PV1(1,13).  1-1,14 

8F10.0/6F10.0 

88,60 

AVI  (1. 14),  1-1,14 

8F10.0/6F10.0 

61,62 

PV1(1,14),  1-1,14 

8F10.0/6F10.0 

63 

PACT  8.  FACT  8,  FACT  4,  FACT  3, 

5F10.0 

FACT  2 

64 

NDA 

12 

68 

TITLE 

20A4 

66 

IPRNT,  El,  1L18T,  IMBAL,  IPI,  lAFT, 

911 

IPNCH,  ICHBCK 

67 

FlN(l.J),  J-1,8 

8F10.0 

68 

PIN(2,J),  J-1,6 

8F10.0 

68 

FIN(3,J),  J-1,8 

8F10.0 

70 

PIN(4,J),  J-1,8 

8F10.0 

71 

FIN(5,J),  J-1,8 

BFIO.O 

72 

F1N(6,J),  J-1,8 

8F10.0 

73 

FIN(7,J),  J-1,8 

BFIO.O 

74-83 

((COMP(J,I),  1-1,14),  J-1,8) 

BF10.0/6F10.0 

84 

N0(1),  1-1,6 

61 

86 

TEMP  (I),  1-1,6 

BFIO.O 

86 

VKM(I),  1-1,8 

BFIO.O 

87 

VKM8 

FIO.O 

86 

BAQ(I),  1-1,6 

BFIO.O 

68 

E(9A(1),  1-1,6 

BFIO.O 

? 

t. 
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Pnceilliii  pigi  Mink 


Table  J1 

Nitration  snapshot 
B-Line  - 8/18/71  - 1:30 
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e c iH  I 


eo  cn  in  e <B 

N in  ^ in  iH 

m o in  N ca 

e 09  <H  T 


(0  C4  r%  m eo 

eo  cn  09  tn  C4 

■ • ■ 4 • 

e CO  cn  « cn 

e eo  H I 


lo  cn  CO  eo  09 

M 09  eo  m M 

eo  <o  cn  (o  e 

09  eo  rH  I 


09  in  o <H  e* 

w CO  (0  H eo 

eo  eo  o w 09 

eo  t«.  »H 


eo  09  (0  « 
CN  CO  eo  09  in 


C4  CO  in  CO 
tH  CO  04 


CO  C4  cn  o iH 

CO  C4  CO  e er 


o*  r*  cn  P4  m 
CO  >H  C-4 


llil 


H H C 


Nitrobociy  0.38  0.49  7.12  23. IB  16.12  19.04  20.27 

Acid  Density  1.619  1.716  1.786  1.798  1.790 


(Cont'd) 


Trinitrobenzeae 


Table  J1  (Cont'd) 


Table  J1  f 


2,3,4-TNI 
TOTAL  Tf 


Tabl«J2 

Pinal  valuta  of  paramatara 


Paranatar 

Codad  Nama 

Value 

Unlta 

’*12 

AK12 

B8.8 

dlmanalonlaaa 

^3a 

AK3S 

1.8x10* 

hr“^ 

‘'aa 

AK48 

1.0x10* 

hr“^ 

k_. 

AK8B 

124.8 

hr“^ 

HBa 

’'ea 

AK68 

508.8 

hr"^ 

^7 

AK7 

.OS 

hr"^ 

AK8 

2.1 

hr"^ 

**10 

AKIO 

.001 

hr"^ 

AKA 

.2 

(ft*)  (Ib-molaa)"^ 

^ql 

8H 

.01 

(lb-molaa)(ft*)"^ 

K 

AKBQ 

50. 

dlmanalonlaaa 

aq2 

K 

8J 

6.0x10* 

(Ib-molaa)*  (ft*)  ^ 

aq3 

K 

SC 

0. 

(ft*)  (Ib-molaa)  ^ 

eq4 

«21 

B21 

88. B 

(Kcal)  (Ib-mola)  ^ 

®3a 

E38 

2.2x10* 

(Koal)  (lb-mola)~^ 

»4a 

«6a 

E48 

1.0x10* 

(Koal)  (Ib-mola)  ^ 

E58 

9.3x10* 

(Koal)  (Ib-mola)"^ 

^6a 

E6S 

8.5x10* 

(Kcal)  (Ib-mola)"^ 

«7 

E7 

5.0x4x10* 

(Kcal)  (lb~mole)  ^ 

! 

f 


-« 
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Parameter 


Coded  Name 


Value 


Unite 


*8 

E8 

7.0x10* 

(Koal)  (Ib-mola)  ^ 

*10 

BIO 

1.0 

(Koal)  (lb-mola)~^ 

k^a 

in 

AKMA 

2.7SX10* 

[ (Ib-molea)  (hr)"’,^  X 

|(ft*)(ft»)"^] 

«V>p 

AKMAP 

1.58x10* 

Idb-moleel  (hr)"Nft*)"^]  *< 
[(ft»)(ft*)"^] 

AKMATL 

8.0x10* 

t (lb-mole)  (hr)“^(ft*)“^J  x 
t(ft*)(ft»)“‘] 

1 

SI 

1.51 

dlmenalonlesa 

PD 

.68 

dlmeneionleei 

RP 

RP 

.6 

dlmeneionleee 
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APPENDIX  K 

StMdy  Stat«  OptimiMtIon 


42B 


PriNdliii  pan  tlmk 


T«bl«  K1 

lnd«p«nd<int  variablas 


Variabla 

Waak  Nitric  Acid  to  lA  6 

Waak  Nitric  Acid  to  IB  20 

Waak  Nitric  Acid  to  2 6 

Fuma  Raoovary  Nitric  Acid  to  lA  IB 

Puma  Raoovary  Nitric  Acid  to  IB  IB 

Yallow  Watar  to  lA  22 

Yaliovr  Water  to  IB  28 

Yellow  Watar  to  2 17 

Praah  Watar  to  lA  28 

Fraah  Watar  to  2 27 

Strong  Nitric  Aold  to  8A  4 

Strong  Nitric  Acid  to  3B  12 

Strong  Nitric  Acid  to  4 3 

Strong  Nitric  Acid  to  8 11 

Strong  Nitric  Acid  to  6 2 

Oleum  to  6 1 

Concentration  of  Weak  Nitric  Acid  13 

Concentration  of  Fume  Recovery  Nitric  Acid  10 

Concentration  of  Strong  Nitric  Acid  D 

Sulfuric  Content  of  Oleum  8 

Tamparature  ofNitrator  lA  7 

Temperature  oi  Nitrator  IB  26 

Temperature  of  Nitrator  2 21 

Temperature  of  Nitrator  3A  14 

Temperature  of  Nitrator  3B  16 


J 


■i 


I 


I 

t 
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Tabu  K1  (Cont'd) 


Variable 

Temperature  of  Nltrator  4 
Temperature  of  Nltrator  5 
Temperature  of  Nltrator  6 
Reoyele  Valve  Notch  Setting  In  Separator  I 
Reoyole  Valve  Notoh  Setting  In  Separator  2 
Reeyole  Valve  Notch  Setting  In  Separator  3 
Reoyole  Valve  Notoh  Setting  In  Separator  4 
Reoyole  Valve  Notoh  Setting  In  Separator  5 
Reoyole  Valve  Notoh  Setting  In  Separator  6 


3 


1 . I DNT  in  Crude  TNT  mutt  be  lese  than  .<i 

2.  I DIaaolved  NItrobody  in  Spent  Acid  < .5 

3 . % DNT  In  DIaaolved  NItrobody  < S3 . 

4 . Weak  Nitric  Concentration  < 61% 

6.  Puma  Recovery  Concentration  < 86% 

6 . Strong  Nitric  Concentration  < 99% 

7 . % H18O4  in  Oleum  > 87% 

6.  Wt  % Actual  Nitric  (AN)  in  atagea  3-6  < 16  | 

9.  AN(S)>7  j 

10.  AN(4)  >6 

11.  AN(8.6)  > 10  j 

12.  AN(1)  <6  ' 

13.  AN(2)  <11 

14.  AN(1)>2.8  I 

16.  AN(8)>4.3  I 

16.  20.  < (Temperature  of  Nitrator  lA  (TIA)  1 

17.  20.  < TIB  < 60.  "C  j 

18.  80<  T2  < 70.  "C 

19.  e0.<  T8A<  100.  ’*0 

20.  80.<  T3B  < 100.  “C 

21.  70.<  T4  < no.  «C 

22.  80.<  TB  < 118.  "C 

23.  90.<  T8  < 120.  ’C 

24.  Heat  Load  in  Nitrator  lA  <2,8x10*  BTU/hr 
28 . 66 . < (%  Total  Acidity  of  Spent  Acid)  < 72 , 
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Tabia  K3 


o 


Optimum  oparating  conditloni 

Bails:  59.2  tons/ day  production  rats 
10000  Ib/hr  olsuin  faad  rats 


Nltrator 

U 

IB  2 

SA  SB  4 

8 

8 

Walk  Nitric 

2000 

1800 

Puma  Rsoovsry  Nltrlo 

800 

Ysllow  Water 

1000 

8trcn|  Nltrlo 

800  800  400 

too 

1100 

Olsuin 

10000 

Toluene 

1783 

887 

Temperature 

81 

S3  84 

83  88  90 

08 

no 

SUBS 

Actual  Nltrlo 

Total  Aoldltv 

1 

3.1 

70.8 

2 

8.1 

88.0 

3 

11.1 

93.8 

4 

13.8 

101.2 

6 

13.8 

103.8 

8 

13.9 

104.2 

Ylsld  « 80.2% 

DNT  In  oruds  TNT  « .17% 

Wt  % dlsaolvad  nltrobody  in  spsnt  sold  ■ .3 
Wt  % of  dissolved  nltrobody  which  ii  DNT  ■ 19,7 
Raw  Material  Cost  ■ 6.48thb  TNT 
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Fig  4 Nitrator-separator  combination 


Fig  7 Chefrical  components  considered  in  the  model 
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Int«rfhc(i 


Bulk  Organio 

Phatt 


Organic  Fhaa« 
Pwalstano* 
Layer 


Acid  Phaae 

Reslatanoe 

Layer 


Bulk  Acid 
Phase 


- net  rate  of  dlffusslon  of  species  1 Into  the  acid  phase 

- BOle  fraction  species  1 in  bulk  organic  phase 

- mole  fraction  species  1 in  bulk  acid  phase 

- eqjilllbrluiB  solubility  of  all  organic  species  In  the  acid  phase 
or  mole  fraction  organic  species  In  add  phase  at  saturation 

- mass  transfer  coefficient  In  acid  phase 

- mass  transfer  coefficient  In  organic  phase 

- moles  species  1 In  acid  phase 

- moles  npecles  1 In  organic  phase 

- total  moles  In  sold  phase 

- total  moles  In  organic  phase 

- total  moles  of  nltrobody  In  add  phase 

- ‘*^otal  moles  of  nltrobody  In  organic  phase 

- superscript  denoting  equilibrium  at  the  Interface 


Fig  8 Diffusion  in  nitration  process 
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OF  VARIABLES 


Manipulated  Variable 

1 . Water  to  Stage  1 

2.  Weak  Nitric  to  Stage  1 

3.  Water  to  Stage  2 

4.  Weak  Nitric  to  SUge  2 

5 . Strong  Nitric  to  Stage  3A 

6 . Strong  Nitric  to  Stage  3B 

7 . Strong  Nitric  to  Stage  4 
6.  Strong  Nitric  to  Stage  8 

9 . Strong  Nitric  to  Stage  6 


Controlled  Variable 


Total  Acidity . Stage  1 
Nitric  Concentration . Stage  1 
Total  Acidity , Stage  2 
Nitric  Concentration , Stage  2 
Nitric  to  Sulfuric  Ratio , Nitrator  3A 
Nitric  to  Sulftirio  Ratio.  Stage  3B 

Nitric  to  SulfUrio  Ratio.  Stage  4 
Nitric  to  SulfUrio  Ratio.  Stage  8 

Nitric  to  Sulfuric  Ratio.  Stage  6 


Fig  10  Acid  concentration  control— varlablaa  operated  on 
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